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‘^tot  i s  minimized  by  taking  partial  derivatives  and  setting  them  equal 
to  zero.  Since  only  two  equations  are  needed,  we  can  select  which  partial 
derivatives  to  use.  We  used  the  partial  derivatives  with  respect  to  vg  ancj 
f  given  by: 


Mtot/3ve  =  0 

Mtot/3f  =  0  (24) 

These  equations  are  canplex  and  numerical  methods  were  used  to  solve  them. 

2.3.3  Summary 

The  process  of  selecting,  sizing,  and  designing  an  integrated  electro¬ 
magnetic  propulsion  system  can  be  summarized  as  follows: 

a.  The  mission  must  be  defined  (AV,  M  yf  tm,  a). 

b.  The  thruster  and  primary  power  types  must  be  chosen. 

c.  The  outputs  of  the  propulsion  system  (mp,  ve,  f,  and  t^) 

are  selected  to  minimize  the  propulsion  system  and  to  satisfy  the 
mission  requirements  (  i.  e.  ,  solving  equations  20,  21,  and  24 

simultaneously)  . 

To  accomplish  this  se  lection  me  thodo  1  ogy  we  nee  d  to  dev  el  op: 

a.  Thruster  performance  models, 

b.  Primary  power  source  mass  functions, 

c.  Propellant  system  mass  functions,  and 

d.  ES/PC  mass  functions. 

All  mass  and  performance  models  must  be  developed  as  functions  of  the 
propulsion  system  outputs. 
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Aa.1 


“pay  +  «PP  *  Mespc  *  Vp  sys  =  Mprp/ ( e xp ( AV/ vg )  “  1}  (20) 

In  order  to  solve  this  eqwtioL,  the  mass  of  the  propellant  and  the  dry 
engine  mass  must  be  related  to  the  thruster  ontput  parameters  as  described 
by  Equations  17  to  19. 

Hie  second  equation  derived  from  the  m.ision  requirements  is  that  the 
total  thrust  time,  t  is  related  to  the  mission  time,  tgj*  by: 

td  =  atB  (21) 


where  a  =  the  thruster  duty  cycle. 

The  thruster  duty  cycle  is  defined  by  the  mission,  and  Equation  21  uniquely 
defines  the  total  thrust  time. 

No  further  independent  relations  among  the  thruster  output  parameters 
can  be  derived  from  the  mission.  The  two  remaining  equations  must  be  de¬ 
rived  from  propulsion  system  optimization. 

2 J .2  Propulsion  System  Optimization 

Optimized  design  of  electromagnetic  propulsion  systems  involves  the 
minimization  of  deployed  mass.  The  criterion  of  mass  minimization  applied 
to  system  selection  and  sizing  is  illustrated  as  follows. 

Mathematically,  minimization  of  the  total  propulsion  system  mass  may  be 
thought  of  as  minimizing  a  function  (the  propulsion  system  mass)  of  several 
independent  variables  (the  propulsion  system  outputs)  subject  to  certain 
auxiliary  restrictions  (the  mission  requirements)  The  propulsion  system 
mass  is  the  sum  of  the  subsystem  masses'- 

Mtot  “  Mpp  +  Mespc  ^prp  sys  ^22^ 

which  can  be  expressed  as  a  function  of  the  propulsion  system  outputs  since 
all  the  subsystem  masses  are  uniquely  related  to  these  variables  (Equations 
17  to  19).  Then  Mtot  can  be  written  as: 

Mtot  =  F(Bp've’ f' td)  (23) 


13 


If  these  relations  were  known,  inspection  of  Equations  S,  8  and  12  shows 
that  all  the  major  propulsion  subsystem  masses  (and  other  performance  re¬ 
quirements)  are  uniquely  related  to  the  propulsion  system  outputs  as: 


MPP  =  Mpp(Bp'Ve'f'td) 


(17) 


Mespc  “  Mespc(“p«ve'f'td) 


(18) 


prp  sys 


Mprp  sys(Bp'f'td) 


(19) 


Therefore,  the  performance  characteristics  of  the  major  subsystems 
must  be  known  before  an  integrated  propulsion  system  can  be  selected,  sized, 
or  designed  in  an  optimal  manner.  In  order  to  illustrate  this  more  clearly 
we  next  examine  the  generic  problem  of  propulsion  system  selection  and 
siz  mg. 


2.3  PROPULSION  SYSTEM  SELECTION  AND  SIZING 

The  selection  and  sizing  of  an  optimal  propulsion  system  for  a  specific 
mission  is  a  complex  task.  The  basic  principles  are  simple.  Propulsion 
system  selection  and  sizing  may  be  thought  of  as  the  process  of  choosing  a 
combination  of  nip,  we,  f  and  t  ^  which  meet  the  mission  objectives  and 
satisfy  the  optimization  criteria.  Also,  the  process  of  choosing  the  pro¬ 
pulsion  system  outputs  results  in  a  unique  definition  of  the  size  and  per¬ 
formance  of  each  of  the  propulsion  subsystems. 

Mathematically,  the  process  may  be  thought  of  as  simultaneously  solving 
four  equations  with  four  unknowns  (i.e.,  the  propulsion  system  outputs).  In 
order  to  find  unique  solutions  for  these  output  parameters,  four  independent 
equations  which  relate  the  parameters  to  one  another  must  be  developed.  The 
process  by  which  this  is  accomplished  is  illustrated  in  the  following 
paragraphs. 

2.3.1  The  Mission 

The  mission  provides  us  with  two  of  the  four  required  equations.  The 
first  equation  is  derived  from  the  momentum  change  which  must  be  imparted  to 
the  payload  mass.  Momentum  must  also  be  imparted  to  the  dry  engine  mass  and 
the  propellant.  The  complete  equation  is: 


(11) 


Mtank  ~  klMprp 

and  the  propellant  system  mass  can  be  written  as: 

^prp  sys  =  ^2mp^td  (12) 

where  kj  >  1.  As  can  be  seen  from  equation  12,  thruster  characteristics  do 
not  directly  affect  propellant  mass  requirements. 

2 .2 .4  The  Thruster 

The  thruster  converts  the  electrical  energy  delivered  to  it  from  the 
ES/PC  into  kinetic  energy  of  the  pulse  mass.  The  thruster  is  required  to  do 
this  at  a  frequency  and  for  a  total  time  consistent  with  the  mission 
requirements.  The  manner  in  which  the  thruster  accomplishes  this  task  has  a 
major  impact  on  the  other  subsystems  in  the  propulsion  system.  Thruster 
performance  has  a  first  order  effect  on  the  ES/PC  and  primary  power  systems. 
The  efficiency  of  the  thruster  will  determine  the  energy  storage  and  average 
power  requirements  of  the  system.  The  physics  of  thruster  operation  will 
determine  the  current  and  voltage  levels  and  pulse  duration  required  of  the 
ES/PC. 

Under  'steady'  average  thrusting  conditions  the  thruster  produces  iden¬ 
tical  pulses  and  its  performance  should  be  largely  independent  of  the  pulse 
frequency  or  the  total  operating  time.  This  implies  that  thruster  per¬ 
formance  is  only  dependent  on  pulse  mass  and  exhaust  velocity.  The  modeling 
relationships  which  are  required  may  then  be  expressed  as 


*th* 

Jth 


’Uh(Vve) 

(13) 

Jth(Bp've) 

(14) 

Vth(V  v 

(15) 

tp(Bp've) 

(16) 
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(8) 


Mespc  =  Me spc^Ep' J’ V# f * £p' ld^ 

where  J  =  ES/PC  output  current 
V  =  ES/PC  output  voltage 
f  »  operating  frequency 
tp  »  pulse  duration 
td  »  thrusting  time 

The  efficiency  of  the  ES/PC  system  will  also  depend  on  these  factors  and  may 
be  expressed  as: 

^e  spc  =  ’'espc^p'  J,V’  f'  tp'  cd}  (9) 

The  relationships  described  by  Equations  8  and  9  have  not  been  extensively 
investigated.  The  determination  of  these  relationships  was  one  of  the 
primary  goals  of  this  program. 

Examination  of  Equations  8  and  9  irdicates  that  four  of  the  six  param¬ 
eters  on  which  ES/PC  performance  depends  are  determined  by  thruster  perfor¬ 
mance.  These  four  parmaeters  are  Ep,  J,  V  and  tp.  Thruster  modeling  must, 
therefore,  not  only  address  pulse  energy  requirements  (or  efficiency)  but 
must  also  treat  current,  voltage,  and  pulse  duration. 


The  Propellant  System 


The  primary  functions  of  the  propellant  system  are  to  store  the  propel¬ 
lant  and  to  provide  the  reaction  mass  for  the  propulsion  system.  The 
propellant  mass  is  a  significant  mass  component  of  the  total  propulsion 
system  and  in  some  missions  will  dominate  propulsion  system  mass.  The 
propellant  mass  is  related  to  propulsion  system  outputs  by: 


Mprp  *  mpftd  (10) 

The  propellant  storage  container  mass  may  be  a  considerable  fraction 
of  the  propellant  system  mass,  depending  on  the  propellant  type  (gas  or 
solid)  and  on  the  way  it  is  stored  (i.e.,  gas  propellant  may  be  stored 
cryogeni  cal  ly  or  in  high  pressure  tanks).  In  general,  the  container  mass  is 
expressed  as  a  function  of  the  propellant  mass  as: 
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Inspection  of  Equations  5,  6,  and  7  indicates  that  the  mass  of  the 
primary  power  system  and  the  oatput  power  requirements  are  directly  related 
to  the  propulsion  system  output  parameters,  mpj  Ve,  f  and  td,  plus  the 
thruster  and  ES/PC  efficiencies.  The  ES/PC  efficiency  is  usually  high 
080%)  and  has  a  small  effect  on  the  primary  power.  The  thruster  efficiency 
is  not  as  high  (<60%)  and  has  a  major  impact  on  system  performance.  Thrust¬ 
er  efficiency  is  an  important  aspect  of  thruster  performance  which  must  be 
analyzed,  understood,  and  modeled. 

2  .2  .2  The  ES/PC 

A  generic  ES/PC  is  presented  in  Figure  3.  The  primary  function  of  the 
ES/PC  is  to  store  energy  for  each  thruster  pulse.  In  addition,  the  ES/PC 
system  must  condition  the  output  to  deliver  an  optimum  pulse  to  the  thruster 
(i.e.,  the  correct  current  and  voltage  profiles).  Input  conditioning  may 
also  be  required  of  the  ES/PC  to  permit  efficient  charging  of  the  energy 
store  from  the  primary  power  source. 


(694) 


Figure  3.  Generic  ES/PC. 

The  mass  of  the  ES/PC  is  dominated  by  the  mass  of  the  energy  storage 
component.  The  mass  of  energy  storage  devices  is  directly  related  to  the 
energy  stored.  Therefore,  it  is  expected  that  the  mass  of  the  ES/PC  will 
primarily  depend  upon  the  pulse  energy  requirements.  The  masses  of  the 
input  and  output  conditioning  systems  depend  on  current  and  voltage  levels 
and  to  some  extent  on  frequency,  pulse  duration,  and  total  propulsion  system 
operating  time  (i.e.,  component  life).  Total  ES/PC  mass  may  then  be 
expressed  as: 


9 


Each  of  the  major  subsystems  plays  an  important  part  in  propulsion 
system  performance.  The  performance  characteristics  of  each  subsystem, 
particularly  subsystem  mass,  must  be  investigated  and  understood  in  order  to 
optimally  select,  size,  and  design  integrated  propulsion  systems.  The  ge¬ 
neric  characteristics  of  each  subsystem  and  the  implications  these  charac¬ 
teristics  have  on  the  propulsion  system  performance  modeling  are  examined  in 
the  following  paragraphs. 

2.2.1  The  Primary  Power  System 

The  primary  power  system  is  the  basic  source  of  energy  for  the  propul¬ 
sion  system.  The  output  of  the  source  is  electrical  power  at  some  charac¬ 
teristic  voltage  and  current  level.  The  power  is  delivered  at  a  constant 
level  to  maximize  the  efficiency  of  the  primary  power  source  and  minimize 
its  mass. 

The  mass  of  the  primary  power  source  is  primarily  a  function  of  the 
average  power  level  required.  Primary  power  sources  for  space  application 
have  been  investigated  extensively,  and  their  mass  can  be  expressed  as: 

»PP  •  *,,<»»>  (5’ 

where  Mpp  =  primary  power  mass 
Ppp  =  average  power  level 

Ppp  is  related  to  the  pulse  energy,  E^,  (i.e.,  the  energy  required  by  the 

thruster  for  each  pulse)  by: 

PPP  *  Epf/%spc  (6> 

where  f  *  propulsion  system  operating  frequency 

’lespc  =  the  efficiency  of  the  ES/PC. 

Ep  is  determined  by  the  kinetic  energy  output  and  the  efficiency  of  the 
thruster,  and  is: 


E 


P 


V2e/2" 


th 


where  Titb  =  thruster  efficiency. 


(  7) 
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The  range  of  mission  parameters  nsed  in  this  stndy  is  smnmarized  in 
Table  3. 

TABLE  3 

SUMMARY  OF  MISSION  CHARACTERISTICS 

Payload  Mass  AV  Mission  Duty 

(kg)  (m/s)  Time  Cycle 

Stationkeeping  2000-25000  50-185/yr.  1-10  yrs  0.1-1 

Orbit  Transfer  2000-25000  5800  30-300  days  1 

2.2  INTEGRATED  ELECTROMAGNETIC  PROPULSION  SYSTEM 

A  generic  electromagnetic  propulsion  system  is  illustrated  in  Figure  2. 
It  consists  of  four  major  subsystems  or  mass  components:  a  primary  power 

source,  an  ES/PC,  a  propellant  system  (propellant  and  its  storage  con¬ 
tainers)  and  the  thruster  itself.  All  of  these  subsystems  are  required  to 
produce  the  thruster  output,  a  pulse  of  propellant  mass,  mpp  at  an  average 
exhaust  velocity,  v e-  The  propulsion  system  produces  these  mass  pulses 
repetitively  at  frequency,  f,  for  a  total  time,  td. 

The  parameters  mp,  ve,  f  and  td  are  the  propulsion  system  outputs. 
Therefore,  system  performance  must  be  modeled  with  these  parameters  as  the 
independent  variables. 


(594) 

Figure  2.  Generic  Electromagnetic  Propulsion  System. 
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Figure  1.  Velocity  Increaent  Per  Year  Required  To  Correct  Eccentricity 
Changes  Caused  By  Solar  Radiation. 


possible  (about  30  days)  to  300  days.  The  payload  mass  of  interest  is 
similar  to  stationkeeping  payload  mass  (2000  to  25000  kg). 

TABLE  2 

AV  REQUIRED  FOR  LEO  TO  GEO  TRANSFER2 


— 

Analytic 

SECKSP0T  Results 

Transfer 

Par  ame  ter 

Predi  ction 

w/o  Shadow 

w  /  Shadow 

No  Inclination 

An  (bn/ sec) 

4.674 

4.674 

4.811 

Change 

tf  (days) 

40.0 

40.0 

51.5 

Incl  ina  tion 

Au  (ka/sec) 

5.757 

5.802 

5.814 

Change 

tf  (days) 

40  .0 

40.3 

49.6 

^ Bauer,  David  P. ,  Barber,  Jobn  P.,  Swift,  Hal  1  ock  F. ,  Vsblberg,  C.  Julian, 
'Electric  Rail  Gun  Propulsion  Study  (Advanced  Electric  Propulsion  Technol¬ 
ogy — High  Thrust),'  AFRPL-TR-81-02 ,  1981. 
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A/M  for  *  generic  class  of  space  structures  which  can  be  deployed  in  a 
single  shuttle  launch  is  presented  in  Table  1  and  ranges  from  0.08  to  0.S98 
m^/kg.  The  AV  requirements  calculated  for  this  generic  class  of  space 
structures  are  presented  in  Figure  1. 


TABLE  1 

CHARACTERISTICS  OF  A  GENERIC  CLASS  OF  SPACE  STRUCTURES1 


SINGLE  SHUTTLE 


LAUNCH 

MASS 

AREA/  MASS 

EFFECTIVE 

GENERIC  CLASS 

SIZE 

(kg) 

( m^/kg ) 

A/M(m2/kg; 

I  PLATE 

MEDIUM 

(100 

m) 

1618 

4.014 

0.201 

STRUCTURE 

LARGE 

(250 

m) 

3672 

11.055 

0.553 

W/O  BLANKET 

II  PLATE 

SMALL 

(  30 

m) 

1334 

0  .43  o 

0.438 

STRUCTURE 

MEDIUM 

(100 

m) 

11350 

0.572 

0.572 

W/ BLANKET 

LARGE 

(150 

m) 

24420 

0.598 

0.598 

III  MODULAR 

SMALL 

(  15 

m) 

23  00 

0.165 

0.091 

ANTENNA 

MEDIUM 

(  60 

m) 

8375 

0.433 

0.113 

LARGE 

(200 

m) 

18017 

1.980 

0.236 

IV  SERIES  OF 

SMALL 

(2) 

7500 

0.826 

0.085 

ANTENNAS 

JEDIUM 

(3) 

11250 

0.802 

0.085 

LARGE 

(4) 

15000 

0.764 

0.084 

The  total  AV  requireme 

nts 

for 

all  stationkeeping  functions  range 

SO  to  184  n/s/yr  mission  times  ranging  from  1  to  10  years. 


2.1.2  Orbit  Transfer 

The  AV  for  orbit  transfer  is  5800  m/  s.  The  AV  required  for  LEO  to  GEO 
transfer  was  investigated  in  a  previous  program  and  the  results  are  pre¬ 
sented  in  Table  2.  The  mission  times  of  interest  range  from  the  minimum 


1Smith,  W.  W.  and  Clark,  J.  P. ,  'Study  of  Electrical  and  Chemical  Propulsion 
Systems  for  Auxiliary  Propulsion  of  Large  Space  Systems,'  NASA  CR-165502, 
1981 . 
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The  AV  required  to  correct  the  payload  inclination  caused  by  solar  and 


lunar  perturbations  is  given  by: 


A V/ year  =  A0AVq£q 


(2) 


where  A0 
A0 
A 

0 

AVGEO 


required  payload  orbital  plane  inclination  change  (rad/yr) 

A  sin  0  cos  0 

constant  whose  value  depends  on  the  body  creating  the  effect 
( 0.74  for  solar  and  1.613  for  Inner) 

inclination  of  equatorial  plane  with  respect  to  the  ecliptic 
or  loner  orbital  planes 

tangential  velocity  of  GEO  payload  (3074.7  m/s) 


Substituting  back  into  Equation  2  we  get: 


AV/year  =  -( 3074 . 7/ ( 360/ 2n ) ) [ ( 0 .74  sin  0  cos  0)solar 

+  (1.615  sin  0  cos  0)lunar^ 

The  AV  ranges  from  41.2  to  51.8  m/s/yr. 

The  AV  required  to  correct  eccentricity  changes  caused  by  solar  radia¬ 
tion  effects  is  given  by: 


AV/year  =  (  3  s  kn/ 21.)  (an/2  sin  (an/2) )  (0/sin-1  0)  (4) 


wh«*re  s 


solar  constant  at  1  AO  *  4.5xl0""<’kg/(m-s^) 


k  -  (1-w)  (A/M) 


a  *  average  reflectivity  of  payload  -  0.3 
A/M  *  area-to-mass  ratio  (m^/kg) 

X  =  mean  angular  velocity  of  earth  motion  around  sun 
=  1 .99xl0-^rad/ s 


a 

0 


■  duty  cycle 


■  eccentricity  ratio 


eP' 


e* 

e  _ 


maximum  allowable  eccentricity 
peak  eccentricity  that  would  occur 
was  circular  and  no  stationkeeping 
appl ied. 


if  initial  orbit 
corrections  were 


(3  defines  the  stationkeeping  accuracy  required  and  is  set  to  zero  to  calcu¬ 
late  the  maximum  AV.  AV  is  a  function  of  a  and  A/M%  a  may  be  as  high  as  1. 
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SECTION  2 


ES/PC  ANALYSES  APPROACH 

The  ES/PC  will  be  examined  within  the  context  of  an  integrated  propul¬ 
sion  system.  The  process  of  propulsion  system  selection  begins  with  identi¬ 
fication  of  the  mission.  The  propulsion  system  is  then  sized  to  accomplish 
that  mission  in  the  most  economical  manner  (the  propulsion  system  is 
optimized  for  the  mission).  Optimized  design  can,  in  principle,  always  be 
accomplished  if  criteria  for  optimization  are  established.  We  assumed  that 
the  optimization  criterion  is  minimization  of  the  integrated  propulsion 
system  mass. 

2.1  MISSION  CHARACTERISTICS 

All  space  propulsion  systems  are  required  to  accomplish  the  same  gen¬ 
eric  task:  imparting  a  velocity  increment,  AV,  to  a  payload  mass,  Mp#yj  in 
some  mission  time,  tD.  These  mission  parameters  are  the  basic,  independent 
variables  which  constrain  propulsion  system  selection  and  sizing,  and  to 
which  propulsion  system  performance  must  be  related.  The  range  of  mission 
parameters  of  interest  must  be  identified. 

2.1.1  Stationkeeping 

On-orbit  stationkeeping  of  geosynchronous  payloads  is  required  to  re¬ 
move  accumulated  disturbances  caused  by:  longitude  drift  caused  by  triaxi- 
ality  of  the  earth,  long  term  increase  in  inclination  of  the  payload  orbital 
plane  due  to  lunar  and  solar  perturbations,  and  change  in  eccentricity 
caused  by  solar  radiation  effects.  Stationkeeping  is  accomplished  by  sup¬ 
plying  external  torques  which  are  provided  by  the  propulsion  system.  The 
AV  imposed  on  the  propulsion  system  by  these  disturbances  were  calculated. 

The  AV  due  to  longitudinal  drift  is  expressed  as: 

AV/year  =  (1.72  sin  (2A+15))  (1) 

where  X  is  the  longitude  in  degrees  East  of  Greenwich.  The  AV  range  is  from 
0  to  1.8  m/s/yr. 
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SECTION  1 


INTRODUCTION 

An  energy  storage /power  conditioning  system  (ES/PC)  is  the  interface 
between  an  electromagnetic  thruster  and  a  primary  power  source  in  a  space 
propulsion  system.  An  ES/PC  accepts  the  relatively  low,  continuous  power 
output  from  a  primary  source,  converts  it  to  a  series  of  high  power  pulses 
and  delivers  it  to  the  thrusters.  Electromagnetic  thrusters  require  peak 
power  levels  well  above  the  average  power  level  required  to  accomplish 
mission  objectives.  The  ES/PC  is  required  to  minimize  the  primary  power 
source  requirements  while  maximizing  thruster  performance. 

Development  of  electromagnetic  propulsion  systems  has  continued  for 
over  twenty  years.  Most  of  the  work  has  concentrated  on  the  thrusters.  As 
thruster  technology  matures,  more  emphasis  must  be  placed  on  the  ES/PC  and 
on  system  integration.  This  report  describes  the  results  of  a  program 
conducted  at  LAP  Research,  Inc.  to  establish  guidelines  by  which  program 
planners  can  optimally  select  and  evaluate  the  performance  of  integrated 
propulsion  systems.  The  guidelines  allow  the  program  planner  to  make 
tradeoffs  among  the  propulsion  system  components  and  compare  alternative 
concepts. 

The  propulsion  systems  considered  in  this  program  covered  a  range  of 
missions,  thrusters,  and  primary  power  options.  The  missions  were:  ten 

year  stationkeeping  and  orbit  transfer  from  low  earth  orbit  (LEO)  to  geosyn¬ 
chronous  orbit  (GEO)  of  large  payloads  (up  to  25  metric  ton).  The  thrusters 
were:  the  Teflon  pulsed  plasma  (TPP),  the  magne  topi  asmadynamic  (NPD),  the 

pulsed  inductive  (PIT),  and  the  electric  rail  gun  (ERG).  The  primary  power 
options  were:  solar  cell  arrays,  solar  thermophotovol  taic,  nuclear  ther¬ 
mionic,  nuclear  thermoelectric,  nucl  ear  Rankine  cycle,  nucl  ear  Bray  ton  cycle 
and  nuclear  magne tohydrodynam ic  (MHD). 

Selection,  sizing,  and  optimization  of  an  ES/PC  can  be  achieved  if  an 
optimization  criterion  is  established  and  if  the  ES/PC  is  examined  within 
the  context  of  an  integrated  propulsion  system.  To  develop  the  selection 
guidelines,  we  identified  the  requirements  imposed  on  the  ES/PC  and  paramet¬ 
rically  related  the  ES/PC  performance  to  the  propulsion  system  outputs. 


1 


SECTION  3 


INTEGRATED  PROPULSION  SYSTEM  PERFORMANCE  AND  MODELING 

In  this  section  we  describe  the  characteristics  of  the  subsystems  in 
the  integrated  propulsion  system  and  present  the  performance  and  mass  models 
of  these  subsystems. 

3.1  PRIMARY  POWER  CHARACTERISTICS 

The  primary  power  source  was  assumed  to  contain  all  the  components 
necessary  to  deliver  electrical  power,  at  some  characteristic  voltage  and 
current  level,  to  the  ES/PC.  The  performance  characteristics  of  several 
power  sources  were  studied  and  mass  models  were  developed.  The  characteris¬ 
tics  and  mass  models  are  described  in  detail  in  Appendix  A.  A  summary  is 
presented  in  this  section. 

The  characteristics  of  interest  are  the  output  constraints  (current  and 
voltage),  the  minimum  and  maximum  power  levels  available,  and  the  mass.  The 
characteristics  were  determined  by  conducting  a  literature  survey  and  con¬ 
tacting  technical  experts  to  confirm  the  results  and  make  refinements  where 
needed.  The  literature  contained  considerable  point  design  data  but  very 
little  parametric  or  limit  data.  The  primary  power  mass  is  best  described 
as : 


MPP  “  “sp^pp 

where  me „  i*  the  specific  mass.  The  specific  mass  was  modeled  by  curve- 

sp 

fitting  the  data  gathered  in  the  literature  search  for  power  levels  ranging 
from  2  kW  to  7  MW.  The  results  are  presented  in  Table  4.  All  power  sources 
investigated  in  this  study  can  be  theoretically  designed  for  any  power 
level.  If  a  power  level  lower  than  the  minimum  specified  in  Table  4  is 
required,  the  mass  is  set  to  the  minimum.  For  high  power  levels  0500  kW) 
multiple  shuttle  launches  may  be  required.  These  effects  are  shown  in 
Figure  4. 

The  output  voltage  range  of  these  systems  is  presented  in  Figure  5. 
There  is  a  tradeoff  between  ohmic  losses  in  the  busses  and  insulation  degra¬ 
dation  due  to  the  high  voltage.  The  optimum  output  voltage  cannot  be  deter¬ 
mined  until  the  mission  is  defined. 
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TABLE  4 


SPECIFIC  MASS  OF  THE  PRIMARY  POWER  SYSTEMS 


POWER  SOURCE 

SPECIFIC  MASS  (m  )  MODEL*  MINIHJM  MASS 

(kg) 

_ _ -  -  -  ■ _ 1 _ 

Solar  Cell  Arrays 

0.013  +  40 .61/Ppp,  Ppp  >1000  W 

54 

S  ol  ar 

Thermophot  ovol  taic 

0.005  +  4.6x10-7P  ,  1000  W  <  <  10.000  W 

pp  P  P 

6 

Nncl  ear 

Thermoelectric 

0.021  +  370 . 3/P__,  Pn„  >  1000  W 
pp  PP 

3  92 

Nucl  ear 

Them  ioni  c 

0.018  +  514/P  ,  Pon  >  10,000  W 

pp 

700 

Nuclear  Brayton 
Cycle 

0.019  +  190.6/P  .  P__  >  1500  W 

pp  PP 

220 

Nuclear  Raukine 
Cycle 

0.013  +  178.6/P  ,  P__  >  1000  W 

pp  PP 

191 

Nu.iear  MED 

10-4  +  1990 .1/Ppp,  Ppp  >  100,000  W 

2000 

6 


< 

X 


X 

o 

a. 


Figure  5.  Power  Output  Voltage  Range  of  Primary  Power  Systems. 


3.2  TH ROSTER  CHARACTERISTICS 

The  performance  characteristics  of  the  four  electromagnetic  thrusters 
were  identified  and  their  inputs  were  modeled.  The  characteristics  and 
models  are  described  in  detail  in  Appendix  B.  A  immary  of  the  characteris~ 
tics  and  a  brief  description  of  the  operation  of  each  thruster  is  presented 
in  the  following  paragraphs. 


3.2.1  Pulsed  Inductive  Thruster  (PIT) 

The  PIT  consists  of  a  flat,  spiral  coil  through  which  an  energy  supply 
is  discharged.  A  cross-sectional  view  is  shown  in  Figure  6.  A  thin  layer 
of  propellant  gas  is  injected  over  the  face  of  the  coil.  When  this  propel¬ 
lant  layer  has  achieved  the  desired  density,  the  energy  supply  is  discharged 
through  the  coil. 

A  magnetic  field  is  produced  from  the  current  flowing  in  the  coil  and 
induces  a  sheet  of  current  to  form  in  the  gas.  Current  in  this  sheet  flows 
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Figure  6.  Cross- sectional  View  of  Pulsed  Inductive  Thruster  (PIT). 

in  the  opposite  direction  to  current  flow  in  the  coil,  and  the  two  repel 
each  other.  The  sheet  accelerates  away  from  the  coil  in  the  axial  direction 
and  produces  the  desired  inpul  se.  The  acceleration  of  the  propellant  ceases 
at  sane  distance  from  the  primary  coil,  when  the  magnitude  of  the  primary 
magnetic  field  is  insufficient  to  induce  current  flow  in  the  propellant. 

3.2.2  Teflon  Pulsed  Plasma  Thruster  (TPP) 

The  TPP  thruster  is  illustrated  in  Figure  7,  and  consists  of  a  set  of 
parallel  electrodes  separated  by  several  Teflon  fuel  bars.  The  Teflon 
propellant  is  held  against  the  electrodes  by  a  constant  force  spring.  A 
solid  state  ignitor  plug  is  located  in  the  cathode  of  the  thruster,  between 
the  Teflon  bars,  and  an  energy  supply  is  connected  directly  to  the  thruster 
electrodes. 

In  operation,  the  energy  supply  is  charged  and  supplies  a  potential  to 
the  thruster  electrodes.  This  applied  voltage  is  sustained  by  the  thruster 
until  a  pulse  is  supplied  to  the  ignitor  plug  initiating  a  discharge  between 
the  electrodes.  Ihe  mini-discharge  from  the  ignitor  plug  produces  enough 
electrical  conductivity  between  the  thruster  electrodes  for  the  main  energy 
supply  to  discharge.  An  arc  forms  between  the  electrodes,  ablates  and 
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Figure  7.  Teflon  Pulsed  Plasma  Thruster  (TPP) 


ionizes  material  frcn  the  face  of  the  Teflon  fuel  bars.  The  ionized  propel¬ 
lant  is  then  accelerated  by  gas  dynamic  and  electromagnetic  forces  produced 
from  the  interaction  of  the  arc  current  and  its  self-generated  magnetic 
field.  After  the  energy  supply  has  discharged,  the  cycle  can  be  repeated. 


3.2.3  Magnetoplasmadynamic  Thruster  (MPD) 

The  MPD  thruster  has  a  coaxial  cylindrical  geometry  as  shown  in 
Figure  8.  It  consists  of  an  outer  cylindrical  anode  and  an  axial  cathode. 
The  back  wall  of  the  thruster  chamber  is  an  insulating  backplate  which 
contains  propellant  injection  holes.  A  power  supply  is  connected  across  the 
electrodes.  Propellant  is  injected  into  the  chmber  through  the  backplate 
at  a  constant  flow  rate.  The  presence  of  the  propellant  gas  between  the 
electrodes  causes  a  breakdown,  forming  a  plasma.  The  high  radial  plasma 
current  (104A)  induces  a  toroidal  self-generated  magnetic  field  which  accel¬ 
erates  the  ionized  gas.  After  a  short  time  into  the  discharge  (10  us),  the 
current  distribution  in  the  plasma  becomes  constant.  The  stable  current 
then  continues  to  ionize  and  accelerate  propellant  gas  as  it  enters  the 
chamber.  Very  high  power  (>1  MV)  is  required  to  operate  this  thruster 
efficiently,  and  only  pulsed  (or  repetitive)  operation  has  been  considered 
here.  The  term  'qua si- steady1  is  used  to  describe  this  thruster  operation. 
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Figure  8.  Magne  topi  a  sin  adynamic  Thruster  (MPD) 


3.2.4  Electric  Rail  Gun  Thruster  (ERG) 


The  acceleration  process  of  an  electric  rail  gun  is  very  similar  to 
that  of  a  Teflon  pulsed  plasma  thruster.  It  is  a  parallel  rail  device  as 
shown  in  Figure  9.  The  mass  to  be  accelerated,  a  'pellet'  which  contains  an 
electrically  conducting  armature  at  its  rear,  is  placed  in  the  bore  of  the 
gun.  The  presence  of  the  pellet  is  the  major  difference  between  the  ERG 
and  TPP.  The  electric  rail  gun  does  not  rely  on  mass  ablation  for  arc  (or 
armature)  formation  or  for  reaction  mass. 

In  operation,  an  electric  potential  is  applied  across  the  electrodes  at 
the  breech  end  of  the  gun.  Current  flows  down  one  rail,  across  the  con¬ 
ducting  armature,  and  back  down  the  other  rail.  The  interaction  of  the 
armature  current  with  the  self-generated  magnetic  field  caused  by  current 
flowing  in  the  rails,  accelerates  the  armature  and  projectile  down  the  bore 
and  out  of  the  gun  producing  the  desired  impulse.  The  cycle  is  repeated  as 
required. 


3.2.5  Thruster  Performance  Modeling 

For  each  combination  of  pulse  mass  and  exhaust  velocity,  there  is  only 
one  thruster  geometry  which  has  maximum  efficiency.  It  is  this  thruster 
which  is  of  interest  here.  Thruster  performance  modeling  could  be  readily 


Projectile 


Figure  9.  Electric  Rail  Gan  (ERG). 


accampl ished  if  a  large  body  of  data  existed  (either  experimental  or  analyt¬ 
ical)  for  aaxiniai  efficiency  thrusters  and  if  this  data  contained  values 
for  the  input  requirements  over  a  vide  range  of  mp  and  ve>  A  simple  anl  ti- 
variate  curve  fitting  technique  could  then  be  used  to  develop  the  required 
relationships  for  E  ,  J,  V  and  tp.  This  was  the  approach  need  for  modeling 
the  electric  rail  gun.3  This  simple  date-based  approach  could  not  be  ap¬ 
plied  to  the  other  thrusters,  because  not  enough  data  was  available  over  a 
wide  range  of  mp  and  v^.  An  alternate  approach  was  developed  based  on 
methods  of  extrapolating  thruster  performance  from  existing  performance  data 
nsing  physically  reasonable  assumptions  concerning  thruster  operation. 
Using  these  methods  we  were  able  to  develop  the  required  relationships 
between  the  thrusters  inpnt  requirements  and  the  propulsion  system  outputs. 
The  results  of  thruster  performance  modeling  are  presented  in  paragrsphs 
3. 2. 5.1  through  3. 2. 5. 4,  along  with  the  assumptions  made  to  develop  the 
model s. 


3 .2.5.1  Pulse  Energy 

We  asstased  that  the  pulse  energy  required  (for  any  of  the  thrusters)  is 
the  sum  of  three  components:  (1)  kinetic  energy.  Egf'  ^  tonizstion 


’’Bauer,  D.  P.,  Barber,  J.  P.  and  Vahlberg,  C.J.,  The  Ei  ear  B»i.  Gan  for 
Space  Propulsion,'  NASA-CR-165312,  Final  Report.  Eebraarv  1  v  (j ; . 
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energy',  Eion,  end  (3)  dissipation  energy,  Edis.  Then  the  pulse  energy  can 
be  expressed  as: 


EKE+  Eion 


E 


di  s 


(26) 


The  kinetic  energy  is  the  useful  energy  imparted  to  the  pulse  mass 
during  the  pulse  and  is  given  by: 


EKE  “  Ve/2  (27) 

The  ionization  energy  includes  all  energy  that  is  'frozen'  into  the 
propellant.  This  term  is  not  strictly  ionization  but  includes  depolymeriza¬ 
tion,  ionization,  and  heating  of  the  gas.  This  term  can  be  written  as: 


Eion 


=  Vi 


(28) 


where  e^  is  a  'specific  ionization  energy1  characteristic  of  the  state  of 
the  plasma  at  the  thruster  exit  and  has  units  of  energy/mass. 

The  dissipative  energy  loss  term  includes  both  electrode  losses  and 
radiative  losses  from  the  plasma.  This  term  includes  only  those  'steady' 
losses  that  occur  during  gas  acceleration.  The  dissipative  energy  is  re¬ 
quired  to  maintain  the  plasma  at  temperature  and  does  not  include  ionization 
and  other  one  time  penalties.  The  energy  loss  mechanism  was  characterized 
as  a  resistance  for  both  the  TPP  and  PIT  and  as  a  resistance  and  'electrode 
fall'  for  the  MPD.  The  dissipative  loss  term  can  be  written  as: 


Edis  *  /RJ2dt  +  /Vf Jdt 


(29) 


where  Vf  is  the  electrode  fall  or  voltage  penalty  which  must  be  paid  to 
remove  current  from  the  electrodes.  Vf  was  assmed  to  be  zero  for  the  TPP 
and  PIT.  Using  the  definition  of  the  dissipation  energy  and  assuming  that: 


(a)  The  resistive  dissipation  in  the  thruster  is  primarily 
due  to  plasma  resistance, 

(b)  The  apparent  dissipative  resistance  remains  constant  during 
accelera  tion, 

(c)  The  acceleration  of  the  plasma  is  primarily  magnetic,  and 
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(d)  The  pi  asm*  body  scales  geometrically  with  pulse  mass  (i.e.,  ratios 
of  dimensions  remain  fixed  as  the  thrnster  size  is  varied). 

the  dissipation  energy  for  the  PIT  and  TPP  can  be  written  as 


E 


dis 


(30) 


and  for  the  MPD  as: 


Edis  =  V$/2vetp  *  VVe)1/2tp  (31) 

See  Appendix  B  for  the  complete  derivation  of  Equations  30  and  31.  The 
terms  k,  k^#  and  are  constants  wuich  were  determined  from  available 
performance  data  and  are  shown  in  Table  5. 


TABLE  S 

PULSE  ENEBGY  COEFFICIENTS 


ei 

(MJ/kg) 

CONSTANTS 

ist 

IONIZATION 

(MJ  /  kg) 

PIT4 

88.8 

k  =  145 

Argon  40 

TppS ,6 

120.0 

k  *  140 

Teflon  90 

MPD7 

28.0 

k-i  -  1206 
k2  -  68,600 

Argon  40 

4 

Lovberg.  R.  H.  and  Dailey.  C.  L. .  'Large  Inductive  Thruster  Performance 
Measurement',  AIAA-81-0708R,  1981. 

^Huberman,  M.  N.  and  Zafran,  S. ,  'Pulsed  Plasma  Propulsion  System/  Spa  ce  cr  af  t 
Design  Guide1,  AFRPL-TR-80-38,  Final  Report,  September  1980,  ADA091006. 

^Personal  Communication  -  Dominic  Palumbo  -  Fairchild  Industries,  Inc. 

^Personal  Communication  -  Kenn  Clark  -  Princeton  University. 
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The  efficiency  of  the  thruster  is  defined  is  the  kinetic  energy  divided 
by  the  pulse  energy.  For  the  PIT  and  TPP,  the  efficiency  is: 

n  -  1/(1  +  2  ei/T^  +  2k/(m£/3ve))  (32) 

The  efficiency  for  the  MPD  is: 

H  -  1/(1  +  2 ei/y^  +  2k1/(iJ/2ve)  +  2k2/(iJ/2v^/2) )  (33) 

The  efficiency  for  the  EEG3  is: 

n  -  13.67  m0-104y-0.299  (34) 

3. 2. 5. 2  Current/Pul  se  Time 

An  average  current  can  be  defined  as: 

Javg  =  [/oP  j2dtlt v]1'2  (35) 

which  can  be  further  reduced  by  using  the  electraaagne  tic  impulse  equation: 

■pve  *  (L'/j2dt)/2,  (36) 

where  L'  is  the  axial  inductance  gradient  for  the  device.  Canbining  Equa¬ 
tions  32  and  36  we  obtain: 

Jsvg  *  <2  »pve/(L'tp))1/2  (37) 

Equation  37  must  be  reduced  further  to  obtain  the  nodel  of  interest  (a  func¬ 
tion  of  m  and  v  only).  Therefore  L'  and  t  have  to  be  expressed  as 

tr  C  P 

functions  of  mp  and  ve.  This  is  discussed  in  Appendix  B.  The  resulting 
average  current  and  pulse  durations  are 

tp  «  13  mp1/3/%  (38) 

and 

Jtyg  *  148.25  »p1/3%  (39) 


for  the  PIT. 
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For  the  TPP,  the  expressions  are : 


tp  =  1.5xl0'11ap1-'3ve 


(40) 


and 


Javg  =  1.07xl07®p1/3  (41) 

For  the  MPD,  the  pulse  duration  (tp)  is  an  independent  par  sine  ter  since  the 
MPD  is  operated  in  a  quasi-steady  mode.  The  MPD  current  is: 

Jayg  =  2.42xl03  (®p%)1/2  (42) 

For  the  ERG3,  the  expressions  are: 

tp  =  4.45xl0-6mp1/3ve  (43) 


and 


J 


avg 


2 .2xl06mp1/3 


(44) 


3  .2  .5  .3  Vol tage 

The  thrusters  considered  in  this  study  are  all  el ectromagne ti c  and  are 
'current  driven'  devices  (i.e..  the  acceleration  forces  are  related  to  cur¬ 
rent,  not  voltage).  These  devices  do  not  require  a  high  voltage  to  accel¬ 
erate  the  pulse  mass  (although  sene  require  a  relatively  high  voltage  for 
the  formation  of  the  current  carrying  plasma  and  to  produce  the  required 
short  pulses).  The  voltage  is,  therefore,  a  derived  parameter  in  electro¬ 
magnetic  thrusters  and  is  indicative  of  the  impedance  of  the  device. 

The  average  voltage,  Vavg.  i*  computed  from  the  pulse  energy,  the  pulse 
time,  and  the  average  current  and  can  be  written  as: 

V,  *  v'V-i1  <45> 

Substituting  Ep,  tp  and  Javg  ^  c°rr*sP°nding  equation  for  each 

thruster  we  obtain  the  following. 
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For  the  PIT, 


V4vg  =  2.59x  10_4Bp1/3ve2  +  4.61xl04mp1/3  +  7.52x10  2ve  (46) 


For  the  TPP, 


V*vg  *  0.311  mp1/3ve  +  7.5xl07»p1/3/v#+87 .2 


(47) 


For  the  MPD, 


V«vg  *  2.07xl0-4«p1/2ve3/2  +  1.16xl04ap1/2ye"1/2+0.498  ve1/2-*-28.3  (48) 


For  the  ERG, 


V 


av  g 


0.032 


0.8 96 v  2.229 
p  « 


(49) 


3  .2  .5  .4  Summary 

Thruster  characteristi cs  and  performance  were  related  to  the  propulsion 
system  outputs.  The  models  deyeloped  for  the  thrusters  (except  for  the  ERG) 
were  based  on  an  understanding  of  the  physical  principles  of  operation  and  a 
number  of  assumptions  concerning  the  performance  of  optimized  thrusters 
combined  with  available  performance  data.  The  models  are  extremely  versa¬ 
tile,  permitting  rapid  evaluation  of  thruster  input  requirements  over  a  wide 
range  of  output  conditions. 

3.3  PROPELLANT  SYSTEM  CHARACTERISTICS 

The  propellant  system  mass  is  expressed  as 

“prp  *ys  =  K  Bpftd  (12) 

where  E  is  a  constant  greater  than  one  and  accounts  for  the  storage  contain¬ 
ers  and  propellant  feed  system.  The  storage  container  mass  varies  over  a 
wide  range,  depending  on  the  propellant  type  (gas/solid)  and  the  storage 
media  (gas  propellant  may  be  stored  c  ry  oge  ni  ca  1  1  y  or  in  high  pressure 
tanks).  Thus,  the  exact  storage  container  mass  cannot  be  determined  until 
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the  mission  is  specified.  We  assumed  that  the  container  mass  is  20%  of  the 
propellant.  The  propellant  system  mass  can  be  written  as: 


Mprp  sys  ~  1>^B1pftd 


(50) 


3.4  ES/PC  CHARACTERISTICS 

The  generic  ES/PC  functions  which  were  analyzed  are  presented  in 
Figure  10.  The  approach  used  was  as  follows: 

(1)  Define  the  input  constraints  and  output  requirements  of  the  ES/PC 

(2)  Conduct  a  literature  search  and  identify  ES/PC  concepts  and 
components  capable  of  performing  the  defined  functions. 

(3)  Analyze  each  component  independently  and  identify  the  critical 
parameters  and  any  operational  constraints. 

(4)  Generate  'skeletal'  performance  models  based  on  fundamental  under¬ 
standing  of  how  each  component  operates. 

(5)  'Calibrate'  the  skeletal  models  by  nsing  existing  data  to  form 
parametric  models  which  describe  the  operational  capabilities  and 
limits  of  each  component. 

(6)  Combine  the  component  models  to  develop  ES/PC  concept  models. 

(7)  Integrate  the  concept  models  with  the  propulsion  system  selecting 
and  sizing  technique. 


Figure  10.  Generic  ES/PC  Functions. 

The  results  of  this  approach  are  presented  in  the  following  paragraphs. 
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3.4.1 


ES/ PC  Input/Output  Constraints 


The  input  constraints  and  output  requirements  of  the  ES/PC  are  essen¬ 
tial  in  identifying  the  concepts  and  components  capable  of  performing  the 
required  functions.  The  constraints  and  requirements  are  presented  in 
Table  6.  These  were  computed  from  the  primary  power  output  and  thruster 
input  requirements.  More  than  one  mass,  size,  and  efficiency  model  may  be 
required  for  the  ES/PC  concepts  and  components,  because  the  parameter  values 
span  such  a  vide  range  (up  to  5  order  of  magnitude). 


TABLE  6 

ES/PC  INPUT  CONSTRAINTS  AND  OUTPUT  REQUIREMENTS 


PARAMETER 

VALUE 

Output 

EP  <kJ) 

0.1  -  400 

V  (kV) 

0.06  -  27 .5 

I  (kA) 

15.0  -  1170 

tp  Us) 

1.0  -  5800 

f  (Hz) 

0.04  -  19 

Input 

ppp  <kW) 

0.2  -  7800 

V  (V) 

28  -  1000 

з . 4 .2  ES/PC  Concepts 

The  energy  store  in  an  ES/PC  is  the  main  component  of  the  system.  The 
input  conditioners  are  selected  to  charge  it  efficiently  and  the  output 
conditioners  are  selected  to  deliver  the  stored  energy  in  the  required  form 

и. e.,  to  shape  the  pulse).  Three  types  of  energy  stores  were  considered: 
capacitive,  inductive,  and  inertial. 

Three  capacitive-based  ES/PC  concepts  were  considered.  The  first  con¬ 
cept  consisted  of  a  simple  capacitor  with  a  closing  switch  as  the  output 
conditioner  (Figure  11).  This  concept  is  ideal  for  thrusters  requiring 
short  pulses  ion  the  order  of  few  microseconds).  The  second  concept  was  an 
LC  circuit  (Figure  12).  This  concept  can  be  used  with  thrusters  requiring  a 
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assumed  the  total  compulsator  mass  to  be  2.5  times  the  rotor  mass,  which  can 
be  e  xpr e  s sed  as : 


the 

the 


The  campulsator  size 
same  length  and  twice 
rotor  tip  speed  by: 


comp 


=  10  Ep / ( av2 ) 


can  be  derived  from  the 
the  radins.  The  rotor 


rotor  size, 
radins,  r. 


(69) 

It  will  have 
is  related  to 


*  =  2J!fnr  (70) 

Combining  Equations  65,  67,  and  70  we  obtain: 

r  =  8Pvtp/n  (71) 

The  length  (b)  of  the  rotor  can  be  calculated  from  the  rotor  mass  and  radius 
and  is  given  by. 


b  =  Mrot/(pnr2) 


(72) 


These  equations  can  be  used  to  calculate  the  mass  and  size  of  a  compulsator 
optimized  for  a  set  of  independent  parameters  (tp>  Ep/  and  f)  specified  by 

the  mission.  During  selection  and  sizing  of  a  campulsator,  it  should  be 

kept  in  mind  that: 

(1)  a  short  pulse  requirement  may  lead  to  an  impr  acti  cal  ly  small  rotor 
radins.  A  practical  minima  pulse  duration  is  100  ps  to  500  ps. 

(2)  a  long  pulse,  on  the  other  hand,  may  lead  to  a  short  rotor  (thin 
disk)  which  is  dynamically  unstable.  In  this  case,  the  tip  speed 
must  be  decreased  until  the  1  ength- to-di  ame  ter  ratio  is  satis¬ 
factory  (say  1)  and  the  mass  and  radius  of  the  rotor  must  be 
recal cul ated. 

(3)  Finally,  surface  current  density  must  not  exceed  the  maximum 

allowable  limit  of  10  MA/m  in  order  to  control  the  armature  wind¬ 

ing  bond  shear  stress.  If  the  current  density  exceeds  this  value, 
the  length  of  the  rotor  must  be  increased  and  the  mass  and  rotor 
radius  must  be  recalculated. 

We  exercised  this  methodology  over  the  range  of  independent  parameters 
of  interest  and  some  of  the  results  are  presented  in  Figures  22  and  23.  The 
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The  pal  se  width  is  related  to  the  electrical  frequency,  f  e#  0f  the 
al  ter na t or  by : 

tp  =  l/(l6fe)  (65) 

because  the  ccmpulsator  compresses  the  half  cycle  pulse  width  of  the  alter¬ 
nator  by  a  factor  of  8  (assuming  electromagnetic  thrusters  act  r  esi  stiv  el  y) . 
The  electrical  frequency  is  related  to  rotational  speed  by: 

fe  =  f«  P/2  (66) 

where  fm  is  the  mechanical  rotor  speed  and  P  is  the  number  of  field  poles. 
The  lector  of  2  can  be  removed  from  Equation  66  if  either  polarity  can  be 
used.  The  electrical  frequency  becomes: 

fe  =  fffl  P  (67) 

The  minimum  pulse  width  can  be  decreased  by  increasing  either  the  rotor 
speed  or  the  number  of  poles  or  both.  The  maximan  rotor  speed  is  stress 
limited  to  150  m  /  s.  The  maximum  number  of  poles  is  limited  by  the  air  gap 
distance  between  the  rotor  and  stator.  As  the  pole-to-pole  separation 
approaches  the  air  gap  distance,  the  excitation  field  leakage  increases.  At 
some  point,  the  leakage  becomes  so  great  that  the  addition  of  poles  is 
futile. 

The  mass  components  of  a  compulsator  are:  the  rotor  (the  energy 

store),  the  back  iron  (return  path  of  the  magnetic  flux)  and  the  field  coils 
(generate  the  magnetic  field).  The  mass  of  the  rotor  is  given  by: 

Mrot  =  4  Es/y2  (68) 

where  v  is  the  rotor  tip  speed.  The  energy  stored  is  calculated  from  the 
ratio  of  pulsed  to  stored  energy  (a)  and  is  limited  to  less  than  5  percent 
by  peak  power  constraints.  The  back  iron  has  the  same  cross-sectional  area 
as  the  rotor  and  its  mass  should  be  approximately  equal  to  the  rotor  (both 
are  ferromagnetic  and  have  similar  densities).  The  field  coils  are  sized  to 
deliver  the  required  output  voltage  and  their  mass  is  relatively  small.  We 
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Figure  21.  Compulsator  Cross  Section.^-4 

We  conducted  s  literature  survey  and  developed  a  methodology  to  calcu¬ 
late  the  size  and  mass  of  a  compulsator.  The  methodology  and  the  limita¬ 
tions  of  the  campulsator  are  discussed  below. 

The  fundamental  limitations  to  the  compulsator  were  investigated  by 
Weldon  et  al.14  The  limits  can  be  divided  into  three  groups:  those  dealing 
with  the  effects  of  load  characteristics,  those  limiting  output  power,  and 
those  limiting  minimum  pulse  width. 

The  load  affects  the  ability  of  the  compulsator  to  compress  flux  and 
produce  short  pulses.  The  compulsator  was  found  to  reduce  the  basic  alterna¬ 
tor  half  cycle  pulse  width  by  a  factor  of  8  for  a  resistive  load  and  by  a 
factor  of  4  for  a  capacitive  load. 

The  peak  output  power  is  limited  by  several  factors.  The  thickness  and 
dielectric  strength  of  the  air  gap  between  the  rotor  and  stator  windings 
have  a  significant  effect  on  peak  power.  The  air  gap  must  be  minimized  in 
order  to  minimize  the  inductance  and  thereby  maximize  output  current.  This 
reduces  the  voltage  breakdown  strength  of  the  gap  (and  the  maximum  permis¬ 
sible  output  voltage  of  the  machine).  A  second  limitation  on  the  output 
power  is  the  shear  strength  of  the  insulation  system  used  to  bond  the  stator 
windings  to  the  stator  and  rotor  windings  to  the  rotor.  The  interaction 
between  the  compulsator  discharge  current  and  the  excitation  field  causes  a 
tangential  force  on  the  conductors.  This  force  slows  the  rotor,  converting 
inertial  energy  to  electrical  energy,  and  results  in  a  tangential  shear 
stress  on  the  insulation  bond.  The  maximise  strength  of  the  insulation  was 
found  to  be  28  MPa  (4  ksi)  which  corresponds  to  a  maximum  surface  current 
density  of  10  MA/ m. 
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where  P0  =  refrigeration  power. 

This  does  not  include  the  mass  of  radiators  required  to  reject  the  heat. 


3  .4  .3  .4  The  Coup  ul  sat  or 

The  compensated  pulsed  alternator  (compul  sa  tor)  is  an  inertial  energy 
store.  It  combines  both  the  energy  store  and  the  output  conditioning  func¬ 
tions  in  one  device.  It  provides  high-voltage,  high-current  pulses  by 
utilizing  the  principles  of  magnetic  induction  and  flux  compression.  The 
compul  sa  tor  is  a  single-phase  alternator  with  a  stationary  field  and  a 
rotating  armature.  The  armature  winding  and  an  identical  stationary  (or 
compensating)  winding  are  connected  in  series.  At  at  least  one  point  per 
cycle  the  inductance  of  the  armature  circuit  is  minimized.  The  variable 
armature  inductance  leads  to  a  flux  compression  action.  When  coupled  with 
alternator  action  this  produces  high  current  pulses  into  the  load.  The 
basic  operation  is  illustrated  in  Figure  20.  A  simplified  cross-sectional 
view  of  a  compul  sator  is  presented  in  Figure  21. 


Maonetlc 

Field 


14Wel  don,  W.  F. ,  Bird,  W.  L. ,  Driga,  M.  D. ,  Tol  k,  K.  M. ,  Rylander,  H.  G. , 
and  Woodson,  H.  H. ,  'Fundamental  Limitations  and  Design  Considerations  for 
Compensated  Pulsed  Alternators',  (XM  University  of  Texas  at  Austin. 
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luil  i«xy  equipment.  He  analysis  end  characterization  of  superconducting 
systems  is  described  in  detmil  in  Appendix  D. 

Superconducting  systems  have  four  major  components,  the  inductor,  the 
current  leads,  the  dewar,  and  the  refrigerator.  He  mass  of  the  inductor  is 
determined  by  the  choice  of  geometry,  the  materials,  the  magnetic  field,  and 
the  structural  support  system.  Then  suitable  choices  are  made  for  these 
parameters  the  specific  energy  of  the  inductor  is  found  to  be: 

wind  =  3-33  x  103J/kg  (61) 

He  current  leads  carry  current  from  the  normal  environment  into  the 
cryogenic  environment  of  the  inductor.  He  current  leads  are  a  major  mass 
component.  He  size  of  current  leads  is  determined  by  a  tradeoff  between 
lead  mass  and  resistive  losses.  He  lead  mass  is  then  related  only  to  the 
current  carried  and  is  given  by: 

“lead  *  200  1  (62) 

He  dewar  provides  the  cryogenic  environment  for  the  inductor  and  is 
designed  to  limit  heat  flow  from  ambient  to  cryogenic.  The  mass  of  the 
dewar  is  determined  by  the  size  of  the  inductor  and  by  the  loads  it  must 
carry.  He  mass  of  the  dewar  and  insulation  is  given  by: 

“dew  “  ^*2  Adew  (63) 

where  Adev  *  the  surface  area  of  the  dewar. 

He  remaining  component  in  a  superconducting  system  is  the  refrigera¬ 
tor.  He  long  duration  of  propulsion  missions  requires  that  a  closed  cycle 
refrigeration  system  be  used.  The  size  of  the  refrigeration  system  is 
primarily  dependent  on  the  heat  losses  into  the  cryogenic  system.  He 
losses  result  from  the  'downleads,'  the  dewar,  and  the  pulse  losses  in  the 
superconductor.  His  energy  must  be  removed  by  the  refrigeration  system  at 
cryogenic  temperature  and  rejected  at  ambient  temperature.  He  mass  of  the 
refrigerator  is  primarily  related  to  the  refrigeration  power  required  at  the 
cryogen  temperature.  He  mass  is: 
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program  are  the  mass,  resistance,  number  of  turns,  dimensions,  and 
efficiency  of  the  inductor. 

The  minimum  mass  constraint  results  in  maximum  resistance  and  minimon 
efficiency.  For  same  applications  the  efficiency  may  be  unacceptably  low. 
Therefore  a  tradeoff  study  must  be  conducted  to  find  the  optimm  size.  A 
higher  efficiency  can  be  obtained  but  there  is  a  mass  penalty.  The  tradeoff 
in  an  integrated  propulsion  system  is  between  the  additional  inductor  mass, 
radiator  ma  ss,  and  primary  power  mass.  Our  electromagnetic  propulsion 
system  optimization  and  sizing  computer  program  conduct  this  tradeoff 
analysis  every  time  an  inductor  size  is  required. 

An  example  of  specific  energy  of  minimum  mass  (minimum  efficiency) 
inductors  is  presented  in  Figure  19.  At  1  ow  current,  the  specific  energy  is 
proportional  to  the  current  and  has  a  slight  dependence  on  the  inductance. 
At  high  current,  the  energy  density  is  independent  of  current  and 
inductance . 


103  104  105  106 
Current  (A) 

Figure  19.  Specific  Energy  of  Minimum  Mass  Inductor',. 


3. 4. 3. 3  Superconducting  Inductors 

The  high  resistive  losses  characteristic  of  normal  inductors  can  be 
overccme  by  using  superconducting  materials.  Superconducting  energy  storage 
systems  are  more  complex  than  normal  systems  and  require  substantial 
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The  characteristic  dimensions  of  a  capacitor  can  be  calculated  from  the 
mass.  He  mass  of  the  dielectric  is  half  the  capacitor  mass.  The  density 
and  volume  of  the  dielectric  are  known.  Assuming  the  capacitor  is  a  cube, 
the  characteristic  dimension  wil  1  be  the  cube  root  of  the  volume. 

He  dielectric  thickness  can  be  calculated  frcm  Equations  S3.  56,  57. 
He  minimus  possible  dielectric  thickness  was  assumed  to  be  2.5  microns. ^ 
If  the  calculated  dielectric  thickness  is  less  than  the  minimum  specified, 
then  it  is  set  to  the  minimum  and  the  mass  of  the  capacitor  recalculated. 

He  efficiency  of  capacitors  was  assumed  to  be  90%,  that  is  10%  of  the 
energy  stared  is  dissipated  in  the  teminations  and  bus  bars. 

3. 4. 3. 2  Normal  Inductors 

Inductors  are  used  in  pulse  power  systems  either  as  energy  stores  or  as 
pul  se-f  oxming  ccmponents.  When  used  as  energy  stores,  inductors  are  attrac¬ 
tive  because  they  provide  a  constant  current  source  characteristic.  His 
characteristic  is  very  important  for  same  electromagnetic  thrusters.  He 
disadvantage  of  inductors  is  that  they  are  'lossy'  and  inefficient.  When 
used  as  pulse-forming  components,  inductors  control  the  shape  and  duration 
of  pulses  delivered  from  capacitors.  Hey  can  be  used  either  in  a  normal 
'LC'  combination  or  in  a  multi-staged  pulse-forming  network.  As  pul  se- 
forming  components,  inductors  carry  current  for  a  short  time  and  therefore 
the  resistive  losses  inherent  in  inductors,  are  not  so  important. 

Sizing  inductors  for  either  function  (energy  storage  or  pulse  forming) 
is  complex  and  must  account  for  thermal,  mechanical,  and  geometrical  con¬ 
straints.  He  effects  of  these  constraints  on  the  size,  mass,  and  effici¬ 
ency  of  normal  inductors  are  discussed  in  detail  in  Appendix  C. 

We  are  primarily  interested  in  minimum-mass  inductors.  He  minimum 
mass  configuration  is  found  by  simultaneous  solution  of  the  equations  de¬ 
scribing  the  constraints  (thermal,  mechanical,  and  ge  erne  trical ).  A  closed 
form  solution  is  not  possible  and  a  numerical  approach  must  be  used.  A 
computer  program  was  developed  to  find  the  solution.  He  inductance,  cur¬ 
rent,  and  conduction  time  were  treated  as  the  independent  variables  and 
input  to  the  computer  program.  He  material  properties  such  as  resistivity, 
density,  and  maximum  stress  must  also  be  supplied.  He  outputs  of  the 
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Equation  59  does  not  account  for  the  aits  of  the  capacitor  case,  the 
terminations,  or  the  impregnating  fluid  (if  required).  The  mass  of  these 
elements  was  estimated  to  range  from  20%  to  80%  of  the  total  capacitor 
m*ss,12.13  we  selected  a  factor  of  50%  to  prowide  the  specific  energy  of  a 
capacitor  as: 

*  (capacitor)  -  3985  N-**-242 

and  for  R<20%  and  (60) 

w  (capacitor)  =  3985  (N.  e°  -046/ 2 .5 ) “0 .241 

for  R>20%. 

The  specific  energy  is  presented  in  Figure  18. 


Life  (Number  of  Charge/Di scharge  Cycles) 

Figure  18.  Specific  Energy  of  Capacitors  Using  Polypropylene  Film. 


^Personal  Communication  -  A.  S.  Gilmour,  SUNTAB,  April  1983  . 

^Personal  Communication  -  Howard  Mauldin,  Sandia  National  Laboratory, 
August  1983. 


N  -  (E/Ebr8-3(2.5e-0-046R)  (57) 

Substituting  Equation  57  into  Equation  55,  ve  get  the  specific  energy 
of  dielectric  for  voltage  reversal  greater  than  20  percent: 

w  -  eo(k/p)(N.e°*046E/2.5)'°-241E§/2 

and  for  voltage  reversal  less  than  20  percent:  (58) 

w  -  so(k/p)>T0-241E^/2 

To  maximize  the  specific  energy  (minimize  tie  mass).  Equation  58  shows 
that  kE^/ P  should  be  maximum.  An  evaluation  of  the  commonly  used  dielectric 
materials  in  capacitors  is  presented  in  Table  8.  The  most  attractive 
dielectric  is  polypropylene  film. 


TABLE  8 

DIELECTRIC  PROPERTIES 


MATERIAL 

k 

Kg/m3 

h 

V/m 

kE^/p 

V-m/kg 

K-Film 

10 

1750 

2.75xl08 

4.32X1014 

PCLYMIDE 

3.6 

1420 

2  .75xl08 

1.9xl014 

POLYSULFONE 

3  .1 

1240 

2.95xl08 

2.17X1014 

POLYCARBONATE 

2  .8 

1210 

2  .75x10s 

1.75xl014 

POLYPROPYLENE 

2.2 

900 

8.66xl0S 

1.8xl015 

TEFLON 

2.1 

2300 

5  , 9xl07 

3.2xl012 

Substituting  values  for  polypropylene  in  Equation  58  yields  the  optimum 
specific  energy  of  the  dielectric  which  is  expressed  as: 

w  -  7970  (N. e°-046R/2 .5)-0,241 


and  for  £>20%  and 


(59) 


w 


7970  N 


,-0.2  41 


for  R<20%. 
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The  parameters  t  and  p  are  material  properties.  For  a  given  dielectric, 
the  energy  density  is  directly  proportional  to  the  square  of  the  applied 
electric  field  on  the  dielectric.  To  maximize  the  energy  density,  the 
dielectric  must  be  stressed  to  the  maximum  allowable  electric  field  deter¬ 
mined  by  the  breakdown  strength  (Eb). 

The  maximum  allowable  field  is  a  quantity  which  depends  on  the  required 
number  of  cycles,  N,  (number  of  charge/discharge  cycles),  and  the  voltage 
reversal  (R)  experienced  during  the  capacitor  operation.8  The  allowable 
ratio  of  charging  to  breakdown  voltage  reflects  the  limit  on  capacitor  life 
due  to  corona  formation.?  When  the  applied  electric  field  is  greater  than 
the  allowable  field,  small  currents  may  arise  in  the  dielectric  and  may 
degrade  the  dielectric  or  lead  to  a  complete  breakdown.1^  Voltage  reversal 
affects  the  capacitor  life  because  it  leads  to  degradation  of  the  dielectric 
and  corona  formation. 

Data  for  high- energy-density  capacitors11  were  used  to  determine  the 
relationship  between  capacitor  life  and  the  applied  electric  field.  This 
rel ationship  is 

N  *  <E/Eb)-8-3  (56) 

Equation  56  shows  that  the  life  of  the  capacitor  is  equal  to  one  cycle  when 
the  applied  field  equals  the  breakdown  strength  and  that  life  increases 
exponentially  as  the  applied  field  is  decreased. 

The  effect  of  voltage  reversal  on  the  capacitor  life  was  analyzed  in 
Reference  9  and  was  found  to  have  a  negligible  effect  if  the  voltage 
reversal  is  less  than  20%.  Voltage  reversal  greater  than  20%  reduces  the 
life  by: 


a 

Sargent,  W.  J.  'Energy  Storage  Capacitors',  Los  Alamos  National  Scientific 
Laboratory. 

^Dailey,  C.  L.,  White,  C.  W.,  'Capacitors  for  Aircraft  High  Power',  AFAR,- 
TR-74-79 . 

i,J£unhardt,  E.  E.,  'Electric  and  Magnetic  Properties  of  Materials',  Texas 
Tech  University. 

^Maxwell,  Data  Sheets  -  Series  M  Capacitors 
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where  C  is  the  capacitance  sad  V  is  the  voltage  to  which  the  capacitor  is 
charged.  For  a  parallel-plate  capacitor,  the  capacitance  can  be  written  as 


C  -  ke0(A/d) 


(52) 


Figure  17.  Parallel  Plate  Capacitor. 

where  *0  is  the  permittivity  of  free  space,  k  is  the  relative  dielectric 
constant,  A  is  the  surface  area  of  one  of  the  plates,  and  d  is  the  distance 
separating  the  plates.  (This  assumes  that  the  linear  dimensions  of  the 
plates  are  much  greater  than  the  plate  separation  distance.) 

The  voltage  can  be  written  as 

V  -  Ed  (53) 

where  E  is  the  electric  field  between  the  plates. 

Substituting  C  and  V  into  equation  51  we  obtain: 

E,  -  s0kAdE2/2  (54) 

'Ad*  is  simply  the  volume  of  the  dielectric.  Using  the  density  of  the 
dielectric  (p),  the  specific  energy  ( w)  of  the  dielectric  can  then  be 
written  as 

»  *  e0(k/p)E2/2 
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(55) 


Pe  rant  a  ti  ons  of  the  input  conditioners,  energy  stores,  and  output 
conditioners  will  give  the  possible  ES/PC  concepts.  This  is  illustrated  in 
Table  7.  The  type  of  input  and  output  conditioner  is  also  constrained  by 
the  primary  power  source  selected  and  the  thruster  (e.g.,  Series/Parallel 
switching  of  the  primary  power  sources  is  only  feasible  when  modular  power 
sources  are  used).  No  input  conditioner  may  be  required  when  thermocycles 
such  as  Rankine  or  Bray  ton  are  the  primary  power  sources  and  they  are  used 
with  an  inertial  energy  store.  The  thruster  requirements  will  limit  the 
type  of  energy  store  and  output  conditioner  which  can  be  selected. 


TABLE  7 

SUMMARY  OF  ES/PC  CONCEPTS 


INPUT 

CONDITIONER 

ENERGY 

STORE 

OUTPUT 

CONDITIONER 

DC/DC  Converter 

CAPACITOR 

Thyratron 

S/P  Switching 

Spark  Gap 

SCR  (as  closing 
switch) 

LC  Circuit 

PFN  Circuit 

Normal  Inductor 

SCR  (As  opening 

Superconducting 
Induct  or 

switch) 

DC/DC  Converter  and 

Cota  pul  sa  tor 

SCR  (As  closing 

Electric  Motor 

switch) 

o  None 

• 

The  next  step  is  to  characterize  the  ES/PC  components  and  develop  the 
mass,  size,  and  efficiency  models.  The  models  will  be  presented  in  para¬ 
graphs  3.4.3  through  3.4.6  and  they  will  be  grouped  by  generic  function. 

3.4.3  Energy  Store  Models 

3. 4. 3.1  Capacitors 

Two  conducting  bodies  (electrodes)  are  separated  by  a  dielectric  to 
form  a  capacitor.  The  energy  is  stored  electrostatically  in  the  dielectric. 
A  parallel  plate  capacitor  is  shown  in  Figure  17.  The  energy  stored  in  an 
ideal  capacitor  is  given  by: 
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Figure  14.  ES/PC  Concept  Using  a  Normal  Inductor  as  the  Energy  Store. 
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Figure  15.  ES/PC  Concept  Using  a  Superconducting  Inductors  as  the  Energy 
Store. 

store  (Figure  16).  A  hcnopolar  generator  was  also  considered  but  was  not 
feasible  because  of  the  relatively  short  pulses  required  by  the  thrusters. 


Figure  16.  ES/PC  Concept  Using  a  Compulsator. 

The  input  conditioners  considered  for  these  concepts  were:  DC/DC  < 

converters,  series/parallel  switching  of  primary  power  modul es,  and  electric 
motors. 

The  output  conditioners  were  switches  and  inductors.  Inductors  are  a 
part  of  the  capacitive-based  concepts  (LC  and  PFN).  The  switches  considered  j 

were  SCR's,  Thyratrons,  and  spark  gaps.  ; 

I 

l 

0  I 


A 
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pulse  duration  of  a  few  hundred  microseconds.  The  third  concept  was  a 
multi-stage  pulse-forming  network  (PFN)  (Figure  13).  This  concept  is  ideal 
for  thrusters  requiring  long  pulses  (on  the  order  of  milliseconds). 


Figure  13.  ES/PC  Concept  Dsing  a  Multi-Stage  Pulse  Forming  Network. 

Two  inductive-based  ES/PC  concepts  were  considered:  one  using  a  normal 
inductor  as  the  energy  stare  (Figure  14)  and  the  second  using  a  superconduc¬ 
ting  inductor  as  the  energy  store  (Figure  15).  Both  systems  are  attractive 
for  pulse  durations  on  the  order  of  milliseconds.  The  major  difference 
between  the  two  is  that  the  superconducting  inductor  offers  a  higher  effi¬ 
ciency  than  a  normal  inductor.  Superconducting  inductors  require  a  refrig¬ 
eration  system  which  offsets  the  saving  from  higher  efficiency.  A  tradeoff 
evaluation  must  be  conducted  before  a  selection  can  be  made. 

The  third  type  of  energy  storage  is  inertial.  Inertial  stores  have  the 
advantage  of  storing  large  amounts  of  energy  in  a  small  volume.  The  system 
considered  used  a  c  ampul  sated  pulsed  alternator  ( c  ampul  sat  or )  as  the  energy 
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results  are  displayed  as  specific  energy  (pulse  energy  divided  by  machine 
mass)  to  permit  comparison  with  alternatives.  Figure  22  shows  the  specific 
energy  as  a  function  of  pulse  energy  for  a  frequency  of  10  Hz  and  pulse 
width  of  1  ms.  Specific  energy  increases  with  pulse  energy  and  is  an  order 
of  magnitude  lower  than  that  of  dielectric  capacitors. 

Figure  23  shows  the  specific  energy  as  a  function  of  the  pulse  duration 
for  a  pulse  energy  of  10~  J  and  a  frequency  of  10  Hz.  The  specific  energy 
decreases  with  increasing  pulse  duration  for  a  fixed  pulsed  energy.  For 
long  pulse  durations  the  rotor  length  becomes  too  short  unless  the  tip 
speed  is  decreased.  When  the  tip  speed  decreases  the  rotor  mass  must 
increase  for  a  fixed  stored  energy  (Equation  68).  The  result  is  increased 
rotor  mass  and  decreased  specific  energy. 

3.4.4  Input  Conditioner  Models 
3  .4  .4.1  DC -DC  Converters 

DC-DC  converters  are  used:  (1)  to  perform  voltage  scaling  to  allow 

high  voltage  power  transmission.  (2)  to  operate  with  near-constant  input 
impedance  to  provide  maximum  power-point  operation  of  the  source,  (3)  to 
provide  constant  power,  variable  voltage  and  current,  energy  store  charging, 
(4)  to  protect  the  system  from  electrical  transients  generated  by  the 
source,  and  (3)  to  perform  all  above  functions  for  the  entire  range  of  power 
level  of  interest. 

In  a  previous  Air  Force  contract,  we  modeled  a  DC-DC  converter  that 
satisfies  all  these  requirements.^  An  actively  controlled  series  resonant 
inverter-converter  is  required.  A  block  diagram  of  a  DC-DC  converter  is 
presented  in  Figure  24.  The  electric  energy  frcm  the  DC  source  is  passed 
through  a  high  frequency  input  filter  which  protects  the  source  and  the  DC- 
DC  converter  from  electrical  transients.  The  next  stage  of  the  converter  is 
a  series  resonant  inverter  with  active  pulse  modulation.  The  inverter 
generates  high  frequency  power  which  is  modulated  in  frequency  and  ampli¬ 
tude  during  generation  Active  modulation  control  provides  system  stabili¬ 
zation,  maximma  power  point  tracking,  and  load  voltage  tracking.  The  high 
frequency  transformer  performs  the  voltage  scaling  function  and  a  diode 
bridge  rectifies  the  power.  The  high  frequency  output  filter  removes  the 
harmonic  content  of  the  internally  generated  AC  power. 
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Figure  24.  DC -DC  Converter  Block  Diagram. 


Hie  mass  of  this  DC-DC  converter  was  found  to  be: 

m  =  10  +  0.057P0*5  +  0.014P0,75  +  10-4P  (73) 

where  m  is  the  mass  in  kg,  and  P  the  inpat  power  in  W.  This  model  is  valid 
for  power  levels  ranging  from  5  kW  to  1  MV.  The  mass  and  specific  mass  of 
DC-DC  converters  are  presented  in  Figures  25  and  26  respectively.  The 
efficiency  of  DC-DC  converters  was  reported  to  range  from  88%  to  97%  and  to 
depend  on  the  specific  design,  not  on  the  power  level  handled.  We  assomed 
an  efficiency  of  92%  for  all  power  levels. 


Figure  25.  DC-DC  Converter  Mass. 
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Figure  26.  DC -DC  Converter  Specific  Mass. 


3.4  .4  .2  Electric  Motors 

Electric  motors  are  required  to  convert  electrical  energy  from  the 
primary  power  source  to  mechanical  energy  for  inertial  energy  stores. 

We  investigated  electric  motors  to  determine  which  type  is  best  suited 
for  space  applications,  and  to  develop  mass  models  over  the  power  range  of 
interest.  We  conducted  a  literature  search  and  contacted  motor  manufac¬ 
turers.  Only  motors  in  the  low  power  range  (150  W  -  75  kW)  have  been 
developed  for  space /aero space  applications.  No  data  was  found  for  power 
levels  above  75  kf.  At  these  low  power  levels,  permanent  magnet  brushless 
DC  motors  and  variable  frequency  AC  induction  motors  have  the  highest  spe¬ 
cific  power.  The  available  perfoxaance  data  for  space-qualified  motors 
were  curve-fitted  and  the  resulting  equation  is: 

w  »  P(0.046  +  0.21  »/P)  (74) 

where  w  is  the  specific  power  in  W/kg,  u  is  the  rotational  speed  in  rad/s 
and  P  is  the  average  power  in  Watts.  The  efficiency  of  these  motors  was 
estimated  to  be  90%.  The  efficiency  of  electric  motors  generally  increases 
with  power  level  and  rotational  speed. 


Electric  notor  performance  data  for  terrestrial  application  were  used 
to  project  the  performance  of  space  qual  if  ied  motors  to  power  levels  higher 
than  75  kf.  The  available  performance  data  for  motors  for  terrestrial 
applications  were  curve-fitted  to  obtain  the  specific  power  and  efficiency. 
The  specific  power  is: 


w  =  4.706P0*2926 


(75) 


and  the  efficiency: 


r\  =  67.96P0*0252  (76) 

The  projected  specific  power  for  space-qualified  motors  is: 

w  =  19.18P0-2368  (77) 

This  model  has  the  same  'slope'  as  the  data  for  motors  for  terrestrial 
applications  and  is  presented  as  a  dashed  line  in  Figure  27.  This  model 

results  in  a  conservative  estimate  for  the  mass  of  electric  motors  for  space 

appl ica  tions. 

3  .4.4.3  Series-Parallel  Switching  of  Primary  Power  Nodules 

Series-Parallel  Switching  (S-P)  of  primary  power  modules  is  an  input 
conditioning  method  which  allows  constant  power  charging  of  either  capaci¬ 
tive  or  inductive  energy  stores,  and  eliminates  the  need  for  voltage 
scaling.  S-P  can  be  designed  to  deliver  any  required  current/ vol  tage  level 
to  the  energy  store.23  S-P  can  only  be  used  with  modular  power  sources  such 
as  photovoltaic,  thermoelectric,  and  thermionic. 


13Ijichi,  L. ,  and  Billerbeck,  V.  J. ,  'Capacitor  Bank  Charging  by  Series- 
Parallel  Switching  of  Solar  Arrays',  16th  International  Electric  Propulsion 
Conference,  1982,  AIAA-82-187  9. 
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Figaro  27.  Specific  Power  of  Electric  Motor: 


To  illustrate  this  concept,  consider  s  solar  array  which  is  nsed  to 
charge  a  capacitor  bank.  The  solar  array  can  be  divided  in  N  modal  es  where 
N  is  an  even  integer.  All  the  modules  are  connected  in  parallel  at  the 
beginning  of  the  charging  cycle  (configuration  A  in  Figure  28).  The  I-V 
characteristic  of  configuration  A  is  curve  'a'  in  Figure  29.  The  charging 
proceeds  from  point  A  to  point  B  along  curve  'a1.  When  the  voltage  reaches 
point  B.  the  modules  are  switched  to  configuration  B  in  Figure  28.  The  I-V 
characteristic  of  the  array  is  then  curve  'b'  in  Figure  29.  The  charging 
proceeds  along  curve  'b'  to  point  C.  Similar  switchings  are  repeated  until 
the  voltage  reaches  the  required  value  with  the  solar  array  in  configuration 
D.  Depending  on  the  number  of  modules  and  switches,  S-P  will  result  in 
capacitor  charging  which  closely  approaches  the  constant  power  curve,  'e'. 

We  assumed  that  the  modular  power  source  had  eight  modules  and  that  the 
switches  used  were  SCB's.  To  model  the  mass  we  need  an  SCR  mass  model  and 
a  method  to  determine  the  amber  of  SCB's  required.  For  eight  modules,  the 
switches  required  and  their  characteristics  are:  8  S(S's  capable  of  1/8  the 
capacitor  voltage  and  eight  times  the  capacitor  current,  4  SCB's  capable  of 
1/4  the  voltage  and  4  times  the  current,  2  SCB's  capable  of  1/2  the  voltage 
and  twice  the  current,  and  7  SCB's  capable  of  the  full  voltage  and  current. 
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Figure  28.  Series-Parallel  Switching  Concept. *5 
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29.  I-V  Characteristics  of  S-P  Switched  Solar  Arrays. 
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His  is  •  total  of  21  switches.  He  mass  and  efficiency  of  each  switch  is 
calculated  separately.  He  mass  of  the  S-P  system  is  the  sum  of  the 
individual  masses  and  the  efficiency  is  the  product  of  the  individual 
eff iciencies. 

3 .4 .5  Output  Conditioner  Models 
3  .4  .3.1  SCSI'S 

SCR' s  are  semiconductor  switches  which  have  demonstrated  long  life  and 
reliability.  He  mass,  efficiency  and  operating  constraints  of  SCR's  are 
discussed  below. 

Fast  switching  SCE' s  capable  of  switching  up  to  17000  amperes  and  2000 
volts  in  a  pulsed  mode  are  available.  For  higher  power  requirements,  SCR's 
can  be  placed  in  a  series  or  parallel  arrangement  to  switch  the  required 
current  and  voltage.  SC&'s  have  a  life  higher  than  the  10^  cycles  and  an 
efficiency  of  better  than  98%. 

SCR  current  and  voltage  capabilities  are  limited  by  the  breakdown 
voltage  of  the  junction  and  by  power  dissipation  in  the  on-state.  Assem¬ 
blies  with  cooling  systems  capable  of  keeping  the  SCE’s  at  a  safe  tempera¬ 
ture  (175°  C)  are  available. 

Manufacturer's  specification  data^  were  used  to  model  the  mass  and 
efficiency  of  the  SCR's.  Two  mass  models  were  generated:  one  for  SCR's 
used  as  closing  switches  (pulse  rated  data  were  used),  and  one  for  SCE's 
used  as  opening  switches  (RMS  data  were  used).  He  mass  models  are  pre¬ 
sented  in  Figure  30.  He  mass  of  SCE1  s  can  be  best  modeled  as  normalized  to 
the  maximum  safe  conduction  current  (m/I).  By  normal  izing  the  mass,  we  were 
able  to  account  for  the  modularity  of  SCS's.  For  s  specific  application, 
the  mass  of  the  SCR  (and  assembly)  is  the  product  of  the  normalized  mass  and 
the  required  current.  He  normalized  mass  models  for  the  SCR's  as  closing 
switches  are: 

m/I  »  1.36  x  10~3  for:  V< 1200V 

and  (78) 

m/I  -  0.67  x  10_6V  x  5.6  x  10“4  for:  V>1200V 
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Figure  30.  SCR'  s  Mass  Models. 


The  models  for  the  SCK's  as  opening  switches  are: 

m/ I  *  2.51  x  10-3  for:  V<1200V 

and  (79) 

m/I  =  2.47  x  10~6V  -4.5  x  10“4  for:  V>1200V 

The  efficiency  of  the  SCR' s  as  modelled  was  about  99%. 

3 .4.5.2  Thyratrons 

A  thyratron  is  a  gas-filled  closing  switch  nsed  in  many  applications 
which  require  high  peax-power  levels  (MW  or  greater)  for  very  short  dura¬ 
tions  (10's  of  jis  or  less).  A  thyratron  is  shown  schematically  in 
Figure  31.  It  consists  of  two  electrodes:  an  anode  and  a  cathode.  A  third 
electrode,  the  control  grid,  is  placed  between  the  main  electrodes.  These 
three  elements  are  contained  in  a  sealed  tube  (typically  ceramic)  which  is 
f  il  1  ed  with  a  gas. 


Anode 


<* 


Sea  1  ed 
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Control 

Grid 


(1017) 


-  Cathode 


Figure  31.  Ihyratron  Scheme  tic. 

Tube  conduction  is  controlled  by  the  grid  voltage.  When  a  positive 
signal  is  applied  to  the  grid,  the  gas  in  the  cathode-grid  space  becomes 
ionized.  Ions  near  the  grid  are  then  accelerated  by  the  anode  field,  caus¬ 
ing  the  entire  volume  to  become  ionized.  The  tube  then  acts  as  a  closed 
switch.  A  thyratron  will  continue  to  conduct  as  long  as  the  gas  remains 
ionized.  During  conduction,  a  positive  space  charge  forms  around  the  grid 
and  effectively  removes  it  from  the  circuit.  Additional  potentials  applied 
to  the  grid  will  not  control  the  anode- to-cathode  current.  In  order  to 
restore  the  tube  to  its  'open'  condition,  the  anode  voltage  must  be  removed 
or  made  slightly  negative  for  a  time  sufficient  to  allow  the  gas  to 
deioniz  e. 

Thyratrons  are  used  with  mercury,  rare  gases,  or  hydrogen  as  the  filler 
gas.  Hydrogen  is  the  most  commonly  used  filler  gas.  Hydrogen  provides  a 
lower  recombination  (recovery)  time  and  the  voltage  drop  across  its  plasma 
is  low  enough  to  keep  the  ion  velocities  low.  Ions  with  high  velocities  can 
damage  the  electrodes. 

The  parameters  affecting  thyratron  performance  are:  peak  current,  peak 
voltage,  pulse  width,  recovery  time,  pulse  frequency  (rep-rate),  forward 
voltage  drop,  and  lifetime  (total  saber  of  pulses). 

The  peak  current  is  limited  by  cathode  arcing  and  grid  quenching.  When 
the  peak  current  exceeds  some  allowable  limit,  an  arc  forms  between  the 
anode  or  grid  and  the  cathode  and  dimages  the  cathode.  High  peak  currents 
can  also  cause  localized  heating  (only  a  small  part  of  the  cathode  is  used 
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for  conduction  in  each  pulse)  which  degrades  the  performance  of  the  thyra- 
tron.  Thyratrons  with  a  peak  current  capacity  of  120  kA  have  been  bail  t.^ 

The  peak  voltage  which  a  thyratron  is  capable  of  'standing  off1  is 
detemined  by  the  distance  between  the  anode  and  the  control  grid.  Under 
repetitive  pulse  operation,  more  complex  mechanisms  begin  to  limit  the 
stand-off  voltage,  such  as  charge  accumula tion  on  the  tube  walls.  One  way 
to  increase  the  stand-off  capabilities  is  to  increase  the  number  of  control 
grids.  Manufacturing  procedures  and  tube  packaging  constraints  place  a 
practical  limit  of  40-50  kV/gap  on  the  stand-off  voltage. 

The  pulse  width  is  limited  by  the  peak  current  and  the  desired  life  of 

the  tube.  High  currents  for  long  times  can  result  in  arc  formation  and  a 

high  rate  of  cathode  vaporization.  Typical  pulse  lengths  at  peak  currents 
are  on  the  order  of  tens  of  microseconds. 

The  pulse  rep-rate  and  the  recovery  time  are  inter— related.  The  maxi¬ 
mum  rep-rate  is  limited  by  the  recovery  time  (time  required  for  a  gas  to 
recombine  and  regain  its  stand-off  capabilities).  Thyratrons  with  rep-rates 
on  the  order  of  kHz  have  been  built.  High  rep-rates  increase  the  average 
heating  of  the  tube  and  affect  the  tube  performance  and  life. 

The  voltage  drop  in  hydrogen  thyratrons  must  be  limited  to  keep  the 

ion  velocity  below  the  critical  velocity  which  damages  the  electrodes. 

Typical  voltage  drops  are  about  100  volts. 

Thyratron  life  is  strongly  related  to  the  allowable  peak  current,  pulse 
width,  and  rep-rate  in  high  frequency  applications.  The  end  of  the  tube 
life  occurs  when  the  cathode  can  no  longer  emit  electrons  into  the  neutral 
plasma  to  begin  conduction.  Typical  hydrogen  thyratron  life  is  greater  than 
1011  pulses. 

A  mass  model  of  the  hydrogen  thyratron  was  developed  by  curve-fitting 
available  data.  The  data  was  derived  from  manufacturers'  specification 
sheets.  The  mass  of  thyratrons  can  best  be  presented  as  specific  power  (w) 
versus  peak  power  (P).  The  relationship  is: 

w  =  5.31  x  108  p-0-21  (80) 


^EG  and  G  Ceramic  Metal  Grounded  Grid  Thyratrons  Data  Sheet,  H5003-1. 
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and  is  illustrated  in  Figure  32. 


Figure  32.  Thyrstron  Specific  Power. 


The  efficiency  of  the  thyratron  can  be  calculated  using  the  voltage 
drop  (V^op).  The  energy  dissipated  is: 

Edis  «  Vdropf/S  (81) 

and  the  efficiency  is  given  by: 


H  -  (Ep  “Edis)/Ep  (82) 

In  most  cases  the  efficiency  is  better  than  93%. 

3  .4  .3  .3  Spark  Gaps 

A  spark  gap  is  a  high-voltage,  high-peak-current  switch  used  in  many 
pulse  power  applications.  The  major  feature  of  this  switch  is  that  it  is 
capable  of  transferring  many  coulombs  of  charge  per  pulse.  The  three  most 
common  types  of  spark  gaps  are  the  direct  over-volted  spark  gap,  the  trig¬ 
gered  spark  gap,  and  the  field  distortion  gap.  All  spark  gaps  consist  of 
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two  electrodes  (an  inode  and  cathode)  separated  by  some  distance  (called  the 
inter-electrode  spacing).  The  elements  are  enclosed  in  a  sealed  envelope 
which  is  either  evacuated  or  filled  with  a  pressurized  gas. 

The  direct  ovei—volted  spark  gap  breaks  down  when  the  stress  in  the 
inter-electrode  gap  exceeds  some  critical  value.  This  critical  value  is 
called  the  self-breakdown  voltage  (SBV).  The  self-breakdown  voltage  is  a 
very  strong  function  of  the  inter—  el  ectrode  spacing,  the  electrode  shape, 
and  the  filler  gas  conditions  (if  a  gas  is  present). 

The  triggered  spark  gap  has  a  third  electrode  placed  inside,  but  sepa¬ 
rate  from,  one  of  the  main  electrodes.  An  electrical  pulse  is  applied  to 
the  third  electrode  generating  an  arc  between  it  and  the  nearby  electrode. 
The  resultant  'mini  discharge1  causes  the  main  electrodes  to  break  down. 
The  breakdown  time  is  a  function  of  the  voltage  across  the  gap.  The  closer 
the  gap  voltage  is  to  the  SBV,  the  faster  the  gap  will  break  down. 

The  field  distortion  gap  has  a  third  element  (usually  a  ring  shaped 
device)  placed  between  the  two  main  electrodes.  This  'trigger'  is  placed  at 
a  point  where  it  is  not  normally  distorting  the  electric  field  lines  in  the 
gap.  When  an  electrical  pul  se  is  applied  to  the  trigger  electrode,  the 
field  in  the  gap  is  distorted  and  a  breakdown  occurs. 

The  parameters  which  affect  spark  gap  performance  are:  peak  current, 

peak  voltage,  pulse  width,  pulse  rep-rate,  total  coulomb  transfer,  recovery 

O 

time,  and  lifetime. 

The  peak  current  is  limited  by  the  required  pulse  width,  the  resultant 
charge  transfer  per  pulse,  and  the  total  umber  of  pulses.  Peak  currents  of 
hundreds  of  kA  are  not  uncommon  in  typical  spark  gap  applications. 

The  voltage  which  a  spark  gap  can  stand  off  is  a  complex  function  of 
the  filler  gas  type  (if  a  gas  is  present),  the  electrode  material,  and  the 
electrode  spacing.  For  extremely  high  voltage  applications,  the  maximxm 
standoff  voltage  is,  in  practice,  limited  by  surface  flashover  outside  of 
the  gap  rather  than  by  the  gap  design  itself.  Peak  voltages  of  millions  of 
volts  have  been  achieved  and  values  as  high  as  100  kV  are  typical. 

The  maximmn  pulse  width  that  a  spark  gap  is  capable  of  is  a  function  of 
the  required  current  and  life.  Pulse  widths  of  tens  of  microseconds  are 
c omm o n.  The  maximum  single-shot  pulse  width  is  limited  by  electrode 
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cooling.  The  minimum  pulse  width  is  set  by  external  circuit  parameters. 
Pulse  widths  less  than  1  ns  have  been  achieved. 

The  allowable  pulse  rep-rate  is  limited  by  the  recovery  time  and  the 
rated  voltage  standoff.  The  recovery  time  is  a  fundamental  limit  to  the 
maximum  pulse  rep-rate.  The  recovery  time  is  limited  by  the  plasma  recombi¬ 
nation  time  and  the  gas  and  electrode  cooling  rates.  For  100%  voltage 
standoff,  a  practical  limit  on  pulse  rep-rate  is  a  few  hundred  pulses  per 
se  cond. 

The  total  charge  transferred  during  the  life  of  a  switch  is  equal  to 
the  charge  transferred  per  pulse  multiplied  by  the  total  number  of  pulses. 
The  charge  transferred  per  pulse  is  directly  related  to  the  current  and 
pulse  width.  Spark  gaps  transfer  a  relatively  large  charge  per  pulse.  This 
changes  the  surface  m  icr  ostr  net  or  e  of  the  main  electrodes  and  affects 
opera  tion.  Very  large  charges  per  pulse  have  been  transferred  in  single¬ 
shot  devices. 

The  life  of  spark  gaps  is  an  area  of  considerable  uncertainty.  Nearly 
all  manufacturers  estimate  spark  gap  life  at  10,000  pulses  at  maximum  rated 
conditions.  Very  little  work  on  extrapolation  of  life  data  has  been  per¬ 
formed.  The  charge  transferred  per  pulse  must  be  decreased  to  increase  life 
(a  spark  gap  can  only  transfer  a  finite  total  charge  during  its  life). 
Increased  life  can  also  be  obtained  by  preionizing  the  gas  with  an  external 
source  to  reduce  the  charge  transfer  load  imposed  on  the  electrodes.  Al¬ 
though  this  technique  has  not  been  thoroughly  tested,  it  is  expected  to 
increase  the  life  to  greater  than  10^  pulses. 

Many  factors  affect  the  performance  and  mass  of  spark  gaps.  After 
analyzing  the  available  data,  we  determined  that  the  mass  of  the  spark  gaps 
is  primarily  a  function  of  the  total  charge  transferred.  With  this 
approach,  we  account  for  the  peak  current,  pulse  width,  and  life  of  the 
spark  gap.  The  resulting  mass  model  is  presented  in  Figure  33  and  is 
written  as: 


m  =  0.0525  cO-4477 


(  83  ) 


where  m  is  the  mass  in  kg  and  C  is  the  total  charge  transferred  in  coulombs. 
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4 .2 .3  Velocity  Increment 

The  optiana  propulsion  system  mass  was  proportional  to  the  velocity 
increment  for  s ta t J onkeeping  missions  (Figure  41).  We  assumed  that  the 
velocity  increment  for  orbit  transfer  missions  is  constant. 


LEGEM 
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5  ?oy*r — Solar  Call! 

6  ES/PC— fVH/X-DC 

~0  10  20  30  40  50  &0 

Exhaust  Velocity  (km/s) 

Figire  41.  Velocity  Increment  Effects  on  Optimum  Propulsion  System  Ma  '  * . 

4 .2 .4  Primary  Power  Options 

Photovoltaic  power  sources  resulted  in  the  lowest  mass  of  optimum 
propulsion  system  for  stationkeeping  systems.  This  is  illustrated  in 
Figure  42. 

Nuclear-based  primary  power  sources  become  competitive  with  photo¬ 
voltaic  for  orbit  transfer  missions.  A  typical  example  is  illustrated  in 
Figure  43. 

4.2.5  ES/PC  Options 


The  capacitive-based  ES/PC  options  resulted  in  the  lowest-mass  optimise 
propulsion  systems  for  stationkeeping  missions.  The  normal  inductive-based 
ES/PC  options  resulted  in  comparable  optimise  propulsion  system  mass.  The 
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In  Table  12  we  are  presenting  two  propulsion  systems  designed  by 
EPSDES.  The  first  has  an  Isp  of  1600  s,  similar  to  the  Martin  Marietta 
system,  and  the  second  has  an  Isp  of  3100  s.  The  propulsion  system  designed 
by  EPSDES  to  operate  at  an  Isp  of  1600  s  has  the  same  mass  as  the  Martin 
Marietta  system.  We  cannot  compare  other  characteristics  of  the  system 
because  of  lack  of  information  on  Martin  Marietta's  system.  But,  we  believe 
that  the  EPSDES-selected  system  operates  with  lower  pulse  mass  and  higher 
frequency  than  Martin  Marietta's.  The  second  system  presented  in  Table  12 
is  the  optimm  system  selected  by  EPSDES.  The  optimum  system  has  30%  less 
mass  than  the  Martin  Marietta  system  and  operates  at  an  Isp  of  3100  s.  Ve 
believe  that  the  optimum  system  has  a  higher  frequency  and  lower  pulse  mass 
than  Martin  Marietta's  system 

4.2  SENSITIVITY  ANALYSIS 

We  analyzed  the  sensitivity  of  optimum  propulsion  system  mass  to  varia¬ 
tions  in  mission  parameters,  primary  power  options  and  ES/PC  options  for 
stationkeeping  and  orbit  transfer  missions.  This  was  done  by  exercising 
EPSDES  over  the  entire  range  of  interest. 

4.2.1  Payload  Mass 

The  optimum  propulsion  system  mass  was  proportional  to  payload  mass  for 
both  stationkeeping  and  orbit  transfer  missions.  This  is  illustrated  in 
Figures  37  and  38  respectively.  The  mass  of  the  propulsion  systems  for 
stationkeeping  missions  ranged  from  5  to  15%  of  the  payload  mass  depending 
on  the  thruster  used.  The  EBG  had  the  highest  propulsion  system  mass  and 
the  other  thrusters  had  similar  masses.  The  mass  of  the  propulsion  systems 
for  orbit  transfer  missions  ranged  froa  2.5  to  3.5  times  the  payload  mass. 
The  EBG  had  the  highest  propulsion  system  mass. 

4 .2 .2  Mission  Duration 

The  optimum  propulsion  system  mass  was  proportional  to  mission  duration 
for  stationkeeping  missions  (Figure  39),  and  increased  exponentially  with 
decreasing  mission  duration  for  orbit  transfer  missions  (Figure  40).  The 
optimum  exhaust  velocity  decreased  with  decreasing  mission  durations  for 
orbit  transfer  missions. 
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TABLE  11 

ORBIT  TRANSFER  PERFORMANCE  COMPARISON 


Parame  ter 

IAP  MPD 
System 
(OPTIMUM) 

IAP  MPD 
System 

MMlS  MPD 
System 

Chem ical 

Advanced 

Systems 

Centaur 

Payload  Mass,  kg 

5000. 

5000. 

5000. 

5000. 

5000. 

AV,  m/  s 

5881 . 

5881 . 

5881  . 

5881  . 

5881 . 

Mission  Time, days 

408. 

408. 

408. 

0.25 

0.25 

Thrust  Eff. ,  % 

50.4 

33.5 

40. 

Total  Eff. .  % 

42.7 

27.6 

36  . 

Average  Power,  kw 

63.* 

47 .3* 

50.* 

I  sp,  s 

3700. 

1600. 

1600. 

482  . 

443  . 

Power  Mass, kg 

859 .6 

655 .5 

OTV  Dry  Mass,  kg 

1709. 

1745 

2958. 

2295  . 

2590. 

Propellant  Mass, kg 

1141 . 

2985 

35  82. 

12145 . 

14120 . 

Propulsion  Mass, kg 

2850. 

4730. 

6540  . 

14440 . 

16710. 

•BOL 


TABLE  12 

ORBIT  TRANSFER  PERFORMANCE  COMPARISON 
(The  specific  energy  of  capacitors  was  assumed  to  be  1.7  J/kg. ) 


Parame  ter 

IAP  MPD 

IAP  MPD 

MM18  MPD 

Chemical  Systems 

System 

(OPTIMUM) 

System 

Sy  s  t  em 

Advanced 

Centaur 

Payload  Mass,  kg 

5000. 

5000. 

5000  . 

5000. 

5000. 

AV,  m/s 

5881. 

5881 . 

5  881  . 

5881  . 

5881  . 

Mission  Time, days 

408. 

408  . 

408  . 

0.25 

0.25 

Th  rus  t  Eff.,  % 

26.6 

18.3 

40. 

Total  Eff. ,  % 

23  .1 

15.9 

36. 

Average  Power, kw 

119.* 

98,6* 

50.* 

Isp,  s 

3100. 

1600. 

1600  . 

482. 

443  . 

Power  Mass, kg 

1550. 

1290  . 

OTV  Dry  Mass, kg 

2973  . 

3055  . 

2958  . 

2295  . 

25  90  . 

Propellant  Mass, kg 

1666. 

3575. 

3582  . 

12145  . 

14120 . 

Propulsion  Mass, kg 

4640  . 

6630  . 

6540. 

14440 . 

16710  . 

•BOL 
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TABLE  10 


STATIONKEEPING  PERFORMANCE  COMPARISON 
(The  specific  energy  of  capacitors  was  assumed  to  be  1.7  J/kg) 


Chemical 

System 


Payload  Mass,  kg 

20978. 

20978. 

19  878 

19878 

19  878 

AV,  m/s/yr 

30.5 

19.5 

100. 

100. 

100. 

Mission  Time,  yr 

5. 

5. 

5. 

5. 

5. 

Thruster  No. 

1 

1 

4 

19 

Average  Power,  kw 

6.39 

7.2 

27.2 

4.9 

Thruster  Eff . ,  % 

11.28 

11.35 

11.3 

54. 

Total  Eff.  .  % 

9.79 

9.73 

9.76 

43. 

Isp,  s 

4000. 

3500. 

3700. 

5000. 

285. 

Frequency,  hz 

80.6 

73.7 

77.1 

2.7 

Pulse  Energy,  J 

69. 

83  .7 

625. 

Pulse  Time,  ms 

0.1 

0.11 

1. 

ES/PC  Mass,  kg 

102. 

114. 

432. 

1476. 

Primary  Power,  kg 

121. 

132. 

506. 

Dry  Mass,  kg 

239. 

258. 

994. 

2579. 

1142. 

Propellant  Mass, kg 

81. 

59. 

280. 

233  . 

3984. 

Propulsion  Mass, kg 

320. 

317. 

1274 

2812. 

5126. 

The  EPSDES  selected  a  propulsion  system  that  has  a  mass  less  than  1/2 
the  Martin  Marietta  system.  The  selected  system  has  an  Isp  of  3700  s  which 
is  twice  the  Isp  of  the  Martin  Marietta  system.  This  led  to  higher  system 
efficiency  and  lower  mass. 

We  suspect  that  the  higher  Isp  was  selected  because  of  a  difference  in 
the  specific  energy  of  the  energy  store  (capacitor).  Therefore,  we  reran 
EPSOES  with  a  specific  energy  of  1.7  J/kg  which  we  believe  that  Martin 
Marietta  used  in  their  analysis.  The  results  of  the  new  run  of  EPSDES  and 
Martin  Marietta  results  are  presented  in  Table  12. 
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TABLE  9 


STATIONKEEPING  PERFORMANCE  COMPARISON 


Parame  ter 

IAP  MPD 

System 

N-S 

Module 

IAP  MPD 

System 

E-f 

Module 

IAP  MPD 
System 

MM18  MPD 
System 

Chemical 

Sy  stem 

Payload  Mass,  kg 

20478. 

20478. 

19878 

19878 

19878 

AV,  m/s/yr 

30.5 

19.5 

100. 

100. 

100. 

Mission  Time,  yr 

5. 

5. 

5. 

5. 

5. 

Thruster  No. 

1 

1 

4 

19 

Average  Power,  k» 

3.52 

4.09 

15.22 

4.9 

Thruster  Eff. ,  % 

26.2 

28.72 

27.5 

54  . 

Total  Eff.  .  % 

21.59 

23  .7 

22.6 

43  . 

Isp,  s 

5000. 

5000. 

5000. 

5000. 

285  . 

Frequency,  Hz 

6. 

5.21 

5.6 

2.7 

Pulse  Energy,  J 

484. 

650. 

625  . 

Pul  se  Time ,  ms 

0.21 

0.25 

1. 

ES/PC  Mass,  kg 

73.2 

81 .5 

309.5 

1476. 

Power  Mass,  kg 

86.4 

93  .8 

360.4 

Dry  Mass,  kg 

171. 

181.75 

705  .5 

2579. 

1142  . 

Propellant  Mass, kg 

63. 

40.25 

206.5 

233  . 

3984. 

Propulsion  Mass, kg 

234. 

222  . 

912 

2812. 

5126. 

4.1.3 .2  The  Orbit  Transfer  System 

Hie  orbit  transfer  mission  specifications  vere: 

Payload  Mass  =  5000  kg 

AV  *  5881  m/s 

Mission  time  *  408  days 

Using  the  EPSDES  computer  code,  ve  sized  the  optimum  propulsion  system  for 

this  mission.  The  EPSDES  and  Martin  Marietta  results  are  presented  in 
Table  11. 
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4.1.3  .1  Hie  St* tionkeeping  System 

Hie  station-keeping  mission  specifications  were: 

Payload  mass  *  19878  kg 

Mission  time  =  5  years 

AV  =  500  m/s 

Martin  Marietta  assmned  f oar  thrusting  modules.  Three  modules  had  five 
thrusters  and  one  module  had  four  thrusters.  Each  module  had  one  ES/PC 
system.  N-S  stationkeeping  required  a  AV  of  61  m/s/yr  and  a  duty  cycle  of 
0.666.  E-W  stationkeeping  required  a  AV  of  39  m/s/yr  and  a  duty  cycle  of 
0.333.  For  EPSDES  modeling  purposes,  we  asstned  that  two  modules  were  used 
for  N-S  and  two  for  E-W  stationkeeping,  that  each  module  had  one  thruster, 
and  that  each  thruster  was  required  to  deliver  half  the  required  AV.  The 
payload  mass  was  also  increased  by  a  mass  equivalent  to  3  thrusting  modules. 
Using  EPSDES  computer  code,  both  E-W  and  N-S  MPD  propulsion  systems  were 
sized.  The  EPSDES  and  Martin  Marietta  results  are  presented  in  Table  9. 

The  system  selected  by  the  EPSDES  has  a  propulsion  system  mass  of  less 
than  1/3  of  the  Martin  Marietta  system.  The  main  difference  between  EPSDES 
and  the  Martin  Marietta  results  was  the  mass  of  the  ES/PC  EPSDES  predicted 
an  ES/PC  which  was  one-fifth  the  mass  of  the  ES/PC  predicted  by  Martin 
Marietta.  This  was  caused  by  the  conservative  specific  energy  of  capacitors 
used  by  Martin  M*':'*tta,  which  was  1.7  J/kg.  As  a  result  of  this  compari¬ 
son,  we  reran  EPSDES  for  the  same  mission  but  we  assmed  that  the  capacitors 
had  the  same  specific  energy  used  by  Martin  Marietta.  The  EPSDES  results  of 
the  new  run  and  Martin  Marietta  results  are  presented  in  Table  10. 

The  new  system  selected  by  the  EPSDES  has  a  propulsion  system  mass  of 
less  than  a  half  of  the  Martin  Marietta  system.  EPSDES  selected  a  propul¬ 
sion  system  that  delivers  small  pulses  one-eighth  the  size  of  Martin 
Marietta's  at  a  lot  Higher  frequency  than  Martin  Marietta  (25  times).  This 
is  to  be  ezpected  since  a  high  price  must  be  paid  to  store  energy.  In 
summary,  the  system  selected  by  EPSDES  uses  the  ES/PC  more  efficiently  which 
leads  to  a  saving  of  half  the  propulsion  system  mass. 


4.1.2  Optimum  Propulsion  Systems 


We  ased  EPSOES  to  examine  the  general  characteristics  of  optimized 
propulsion  systems  for  both  stationkeeping  and  orbit  transfer.  The  mass  of 
the  optimum  propulsion  system  for  stationkeeping  missions  ranges  from  5%  to 
15%  of  the  payload  mass  depending  on  the  thruster  used.  Hie  EEC  had  the 
highest  propulsion  system  mass.  The  minimum  allowable  pulse  mass  of  the  ERG 
is  0.1  g,  3  orders  of  magnitude  higher  than  that  for  the  other  thrusters. 
This  results  in  high  stored  energy  and  a  corresponding  high  ES/PC  mass.  The 
largest  mass  component  of  the  propulsion  system  for  stationkeeping  was  the 
propellant  at  about  50%.  The  optimization  tradeoff  for  stationkeeping  is 
primarily  between  energy  storage  mass  in  the  ES/PC  and  primary  power  mass. 
As  expected  the  ES/PC  and  primary  power  have  about  equal  mass  in  optimized 
systems. 

The  mass  of  the  optimum  propulsion  system  for  orbit  transfer  missions 
ranges  from  2.5  to  3.5  times  the  payload  mass,  with  the  ERG  being  the 
heaviest.  The  ES/PC  mass  was  the  least  massive  component  at  about  15%  of 
propulsion  system  mass.  In  orbit  transfer  missions  the  optimization 
tradeoff  is  primarily  between  propellant  and  primary  power.  As  expected 
these  two  components  had  about  equal  mass  in  optimized  orbit  transfer 
systems. 

Within  the  ES/PC  system,  the  major  components  of  optimized  station¬ 
keeping  systems  had  approximately  equal  mass.  The  ERG  was  an  exception.  The 
major  mass  component  of  the  ES/PC  for  the  ERG  was  the  energy  store  at  about 
80%  of  the  ES/PC  mass.  The  major  mass  component  of  the  ES/PC  for  orbit 
transfer  missions  was  the  radiator  at  about  45%.  The  energy  store  was  30%, 
the  input  conditioner  15%,  and  the  output  conditioner  10%  of  ES/PC  mass. 

4.1.3  Comparison  of  EPSDES  with  Martin  Marietta  Results 

Martin  Marietta18  conducted  two  propulsion  system  point  designs  in 
their  '  MPD  Definition  Study  Program1.  The  first  point  design  was  for  a 
stationkeeping  mission  and  the  second  was  an  orbit  transfer  mission.  We 
exercised  EPSDES  for  the  same  missions.  Both  designs  were  MPD  thruster- 
based.  Our  objective  was  to  compare  the  EPSDES  approach  to  system  sizing 
with  the  more  traditional  approach  used  by  Martin  Marietta. 

18,MPD  Thruster  Definition  Study,'  Martin  Marietta,  1983. 
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MASS  (ton)  «  SYSTEM  MASS  (ton) 


igure  35. 


Propulsion  System  Optimization  (Orbit  Transfer,  TPP) . 


Figure  36.  Optimum  Prooulsion  System  Mass  (Orbit  Transfer,  TPP) 
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SECTION  4 


RESULTS 

Ve  exercised  the  EPSDES  program  extensively.  The  first  series  of 
computations  were  designed  to  examine  the  performance  of  the  program,  to 
establish  the  validity  of  the  solutions  vhich  it  finds,  and  to  determine  the 
general  characteristics  of  optimally  sized  propulsion  systems.  Ve  then  used 
EPSDES  to  conduct  a  sensitivity  analysis.  Ve  examined  the  sensitivity  of 
propulsion  system  size  to  variations  in  mission  parameters  over  the  entire 
range  of  mission  parameters  defined  in  Section  2. 

4.1  EPSDES  VALIDATION 

4.1.1  EPSDES  Optimization 

The  basic  assumption  underlying  the  development  of  EPSDES  was  that  for 
each  mission  and  thruster  combination,  there  is  one  set  of  propulsion  system 
output  parameters  (mpt  and  f)  which  minimizes  the  propulsion  system  mass. 
The  first  task  was  to  verify  this  assmption. 

Ve  selected  an  orbit  transfer  mission  and  a  TPP-based  propulsion  system 
to  check  the  validity  of  EPSDES.  EPSDES  was  run  with  a  fixed  exhaust 
velocity  and  the  pulse  mass  (and  frequency)  was  varied.  The  results  at 
several  exhaust  velocities  are  shown  in  Figure  35.  Each  point  on  the  graph 
represents  an  integrated  TPP-based  propulsion  system  which  will  be  capable 
of  doing  the  required  mission.  There  is  a  clear  minimum  in  propulsion 
system  mass  at  each  exhaust  velocity.  Ve  conclude  that  for  each  exhaust 
velocity,  there  is  an  optimum  pulse  mass  which  minimizes  the  propulsion 
system  mass. 

Figure  35  also  shows  that  selection  of  the  exhaust  velocity  has  a 
pronounced  effect  on  propulsion  system  mass.  'Die  locus  of  the  minima  of  the 
curves  in  Figure  35  are  plotted  in  Figure  36.  There  is  a  combination  of 
exhaust  velocity  and  pulse  mass  (and  frequency)  which  results  in  minimum 
propulsion  system  mass.  This  combination  results  in  the  'optimized'  propul¬ 
sion  system.  We  ran  EPSDES  for  a  wide  range  of  missions  and  an  optimum 
propulsion  system  was  always  found. 
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Figure  34.  Flowchart  of  the  Electric  Propulsion  System  DESign  (EPSDES) 
Canputer  Program. 


3  .5  SYSTEM  INTEGRATION 


To  assess  the  propulsion  system  as  an  integrated  system,  the  subsystems 
and  components  characteristics  and  performance  models  must  be  incorporated 
into  a  system  level  description  of  an  electromagnetic  propulsion  system. 
The  system  integration  can  be  best  accomplished  by  using  numerical 
techniques  because  of  the  complexity  of  the  models  and  the  large  nnmber  of 
options . 

We  developed  an  Electromagnetic  Propulsion  System  Des ign  (EPSDES)  com¬ 
puter  program  to  assist  in  selecting  and  sizing  the  optimum  propulsion 
system  for  a  specific  mission.  EPSDES  is  written  in  FORHLAN.  EPSDES  calcu¬ 
lates  the  combination  of  mp,  v  e«  aud  f  which  minimizes  the  propulsion  system 
mass  and  satisfies  the  mission  requirements.  A  simplified  flow  chart  is 
presented  in  Figure  34. 

EPSDES  is  an  interactive  program  which  prompts  the  user  for  the  re¬ 
quired  input  data.  The  input  data  are:  (1)  the  mission  parameters  Mpay*  AV, 

and  u.  (2)  the  thruster  type,  (3)  the  primary  power  type,  (4)  the  ES/PC 
type,  (5)  the  exhaust  velocity  range  and  increment,  and  (6)  an  initial  guess 
for  the  pulse  mass.  The  output  data  are  the  optimum  subsystem  masses  and 
component  characteristics  at  each  exhaust  velocity  within  the  range 
spe  cif ied. 

The  validation  of  EPSDES  as  a  propulsion  system  optimization  program  is 
presented  in  Section  4.  Also,  we  are  presenting  in  Section  4  a  comparison 
between  EPSDES  results  and  known  point  designs,  and  the  results  of  sensi¬ 
tivity  analyses  of  the  optimum  propulsion  system  size  to  variations  in 
mission  parameters. 
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Total  Coulomo  Transferred  ( C ' 


Figure  33.  Spark  Gap  Mass. 

This  model  assumes  that  the  mass  is  independent  of  the  spark  gap  type  and 
that  the  total  charge  transfer  is  equal  to  the  product  of  the  number  of 
pulses  and  the  charge  transfer  per  pulse.  The  efficiency  of  spark  gaps  was 
estimated  to  be  90%. 

3.4.6  Heat  Rejection  System  Model 

The  ES/PC  system  produces  heat  which  must  be  rejected  to  maintain  the 
operation  and/or  to  keep  temperatures  within  acceptable  limits.  Ve  asstmed 
that  the  heat  is  rejected  through  a  low  temperature  radiator. 

In  a  previous  contract  to  the  Air  Forced  we  modelled  the  specific  mass 
of  a  low  temperature  radiator  for  space  applications.  The  specific  mass  is 
given  by: 


w  =  7.24xlO’V4  (84) 

where  w  is  in  kg/V  and  T  is  the  temperature  of  the  radiator  in  Equation 

84  is  valid  for  temperatures  ranging  from  250°K  to  500°I. 
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Figure  43.  Effects  of  Primary  Power  Options  on  Optimum  Propulsion  System 
Mass  for  Orbit  Transfer  Missions. 
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superconductor  inductive  ES/PC  resulted  in  a  much  heavier  optimum  propulsion 
system.  The  inertial  energy  stores  were  not  feasible  for  the  PIT  and  TPP 
thrusters  because  of  the  short  pulse  duration  required.  A  typical  example 
of  the  effects  of  ES/PC  options  on  propulsion  mass  is  presented  in 
Figure  44. 
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Figure  44.  Propulsion  System  Mass  for  Different  ES/PC  Options  for 
Stationkeeping  Missions. 


For  orbit  transfer  missions,  the  capacitive-based  ES/PC  options  re¬ 
sulted  in  the  least  massive  optimum  propulsion  system.  The  inductive— based 
ES/PC  options  (both  normal  and  superconductor)  resulted  in  comparable  re¬ 
sults.  The  inertial  energy  stores  resulted  in  much  more  massive  propulsion 
systems.  The  propulsion  system  masses  for  various  ES/PC  options  are  pre¬ 
sented  in  Figure  45. 


4.3  POINT  DESIGNS 

To  demonstrate  the  methodology  developed  in  this  program,  we  conducted 
two  point  designs,  one  was  for  a  stationkeeping  mission  and  one  for  an  orbit 
transfer  mission.  The  missions  specifications  and  the  results  of  the  point 
designs  are  discussed  below. 
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Figure  45.  Propulsion  System  Us  s  s  for  Orbit  Transfer  Missions  Using 
Different  ES/PC  Options. 


4.3.1  Stationkeeping  Mission 


Hie  stationkeeping  mission  specifications  mere: 


Urns  ter:  MPD 

Primary  power  source:  Solar  cell  arrays 

Payload  mass:  20,000  kg. 

AV:  875  m/s 

Mission  duration:  7  years 

To  analyze  tkis  mission  we  assumed  tkat: 

(1)  Four  thrusters  are  required.  Two  for  E-T  stationkeeping  and  two 
for  N-S. 

(2)  The  two  thrusters  for  each  function  are  identical. 

(3)  Each  thruster  has  an  independent  ES/PC. 

(4)  The  mission  AV  was  divided  equally  noug  the  four  thrusters. 

(5)  E-W  stationkeeping  was  done  continuously,  i.  e. ,  duty  cycle  -  1. 

(6)  N-S  stationkeeping  was  done  with  a  duty  cycle  “  0.3. 

The  optimal  propulsion  system  for  each  thruster  was  then  determined 
from  the  Propulsion  System  Optimization  Computer  Code.  The  optimum 


propulsion  system  is  presented  in  Table  13.  the  size  and  characteristics  of 
the  thrusters  in  Table  14,  and  the  size  and  characteristics  of  the  ES/PC 
components  in  Table  13. 

The  optimum  propulsion  system  output  parameters  are:  a  pulse  mass  of 
0.3  mg,  an  exhaust  velocity  of  30,000  m/s,  and  a  frequency  of  7  Hz.  These 
output  parameters  resulted  in  a  propulsion  system  which  has  a  mass  of  1200 
kg,  which  is  6%  of  the  payload  mass,  and  a  system  efficiency  of  22%.  The 
average  power  level  required  to  power  the  system  is  19  kV. 

The  low  system  efficiency  is  caused  by  the  low  efficiency  of  the 
thrusters  which  is  about  27%.  The  thrusters  have  a  low  efficiency  because 
of  the  low  pulse  duration  which  is  on  the  order  of  few  hundred  microseconds. 

The  optimw  ES/PC  mass  is  330  kg  which  is  30%  of  the  propulsion  system 
mass.  As  expected,  the  major  mass  components  of  the  ES/PC  have  equal 
weights. 


TABLE  13 

STATION  KEEPING  OPTIMJM  PROPULSION  SYSTEM  FOR  A  POINT  DESIGN 


Parameter 

MPD 

System 

N-S 

Module 

MPD 

System 

E-W 

Module 

Propulsion 

System 

Payload  Mass,  kg 

21000.* 

21000.* 

20000 . 

A  V,  a/s/yr 

31.23 

31.25 

125  . 

Mission  Time,  yr 

7. 

7. 

7. 

Thrust /Engine,  N 

0.07 

0.021 

0.045 

Average  Power,  br 

6.27 

2.80 

18.14 

Thrust  Eff.  ,  % 

30.90 

23  .06 

26  .98 

Total  Eff.  .  % 

25.32 

19  .00 

22.25 

Isp,  s 

4500. 

5000. 

4750. 

Frequency,  Hz 

5.40 

6.97 

6.18 

Pulse  Energy,  J 

938. 

331. 

644.5 

Pulse  Time,  ms 

0.31 

0.18 

0.24 

ES/PC  Mass,  kg 

110. 

62.7 

345.4 

Dry  Mass,  kg 

250.7 

157.5 

816.4 

Propellant  Mass, kg 

103  .3 

92.5 

391 .6 

Propulsion  Mass, kg 

354  . 

250. 

1208. 

The  payload  mass  was  increased  by  an  amount  equivalent  to  the  mass  of  3 
thrusting  modules. 
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TABLE  14 


THRUSTER  SIZE  AND  CHARACTERISTICS  FOR  STATIONKEEPING  POINT  DESIGN 


Parame  ter 

MPD 

for 

N-S 

MPD 

for 

E-f 

Outer  radius,  cm 

2.8 

1.9 

Length,  cm 

11.0 

5.9 

Efficiency,  % 

30.9 

23. 

Mass  flow  rate,  g/s 

0.95 

0.33 

Isp,  s 

4500. 

5000. 

Frequency,  Hz 

5.4 

6.97 

Pulse  duration,  ms 

0.31 

0.18 

Pulse  energy,  J 

957. 

330.7 

Average  Voltage,  V 

196.5 

183. 

Average  current,  A 

15810. 

9889. 

4 .3 .2  Orbit  Transfer  Mission 

Two  thrusters  were  assumed  for  this  mission.  The  required  AV  was 
divided  equally  among  the  two. 

The  optimnm  propulsion  system  for  the  orbit  transfer  mission  is  pre¬ 
sented  in  Table  16,  the  size  and  characteristics  of  the  thruster  in 
Table  17.  and  the  ES/PC  system  in  Table  18. 

The  selected  ES/PC  system  was  a  DC -DC  converter,  a  capacitor  and  a 
hydrogen  thyratron  closing  switch. 

The  optimon  propulsion  system  output  parameters  are:  a  pulse  mass  of 
35  mg,  an  exhaust  velocity  of  26000  m/ s  and  a  frequency  of  10.5  Hz.  The 
propulsion  system  mass  is  20,400  kg  which  is  twice  the  payload  mass.  The 
system  efficiency  is  55%.  The  average  power  level  required  to  power  this 
mission  is  450  kW. 

The  mass  breakdown  of  the  major  components  of  the  propulsion  system  is 
as  follows:  the  primary  power  source  mass  is  10,000  kg,  the  propellant  and 
containers  mass  are  7660  kg,  and  the  ES/PC  mass  is  2600  kg.  The  primary 
power  mass  is  the  heaviest  because  of  the  relatively  short  mission  time. 

The  heaviest  component  of  the  ES/PC  was  the  radiator  with  a  mass  of 
865  kg  (35%  of  the  ES/PC  mass).  The  energy  store  mass  was  about  30%  of  the 
ES/PC  mass,  and  the  input  and  output  conditioners  have  about  equal  mass  (15% 
of  the  ES/PC  mass). 
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e 


0 


Component 

Far ameter 

ES/PC 

ES/PC 

f  or 

f  or 

N-S* 

E-W* 

-  •• 

Capmci  tor 

Efficiency,  % 

98. 

98. 

Mass,  kg 

5.62 

2.76 

Volmne,  m3 

3.12E-3 

1.53E-3 

Dimension,  m 

0.15 

0.12 

Capacitance,  F 

2.48E-3 

9.87E-4 

Voltage,  V 

393  .1 

366. 

Energy  stored,  J 

191.4 

66.1 

Number  of  cycles 

3.58E8 

1.54E9 

Specific  energy. 

J/kg 

34. 

23  .9 

Inductor** 

Efficiency,  % 

92.5 

92.3 

Mass,  kg 

1 .75 

0.78 

Inductance,  H 

0.38E-6 

0.34E-6 

Mean  radius,  cm 

6.5 

5  . 

Length,  cm 

26.1 

19.9 

No.  of  turns 

2.8 

3. 

Thickness,  cm 

0.67 

0.51 

Resistance,  ohm 

5.05E-5 

7  .67  e— 5 

sex 

Efficiency,  % 

99. 

99. 

Mass,  kg 

21.5 

13.5 

Volume,  m3 

29.8E-3 

19.9E-3 

Dimension,  m 

0.31 

0.27 

No.  of  SQ'  s 

6 

4 

No.  in  parallel 

6 

4 

No.  in  series 

1 

1 

DC -DC 

Efficiency,  % 

92  . 

92. 

Converter 

Mass,  kg 

25  . 

18.7 

Average  power,  V 

6259. 

2797. 

V ol ume ,  m^ 

9 .25 E- 3 

6.92E-3 

Dimension,  m 

0.21 

0.19 

Radiator 

Temperature,  Deg. 

E 

300. 

300. 

Mass,  kg 

19.4 

8.7 

Surface  area,  m^ 

6. 

2.7 

Ave.  power  rejected,  W 

1092. 

492  . 

For  one  thruster  module. 
For  one  stage  of  the  PFN 


TABLE  16 


ORBIT  TRANSFER  PROPULSION  SYSTEM 


Pa r am e  ter 

Propul  si  on 
System 
for  1  PIT 

Total 

Propul sion 
System 

Payload  Mass,  kg 

20000.* 

10000 . 

AV,  m/s 

2  900. 

5880. 

Mission  Time,  days 

100. 

100. 

Thrust/Engine,  N 

9.60 

9  .60 

Average  Power,  kW 

225. 

450. 

Thrust  Ef f. ,  % 

62  .3 

62.3 

Total  Eff. ,  % 

55.5 

55  .5 

Isp,  s 

2600. 

2600. 

Frequency,  Hz 

10.5 

10.5 

Pulse  Energy,  kJ 

19  . 

19. 

Pulse  Time,  ps 

16.5 

16.5 

ES/PC  Mass,  kg 

1300. 

2600. 

Dry  Mass,  kg 

7010. 

14020. 

Propellant  Mass, kg 

3190. 

63  80. 

Propulsion  Mass, kg 

10200. 

20400 . 

The  payload  mass  was  increased  by  an  amount  equivalent  to  the  mass  of 
one  thrusting  module. 


TABLE  17 

THRUSTER  SIZE  AND  CHARACTERISTICS  FOR  ORBIT  TRANSFER  MISSION 


Parameter 

PIT 

Diameter,  m 

2.18 

Efficiency,  % 

62.3 

Isp,  s 

2600. 

Frequency,  hz 

10.5 

Pulse  duration,  ps 

16.5 

Pulse  energy,  kJ 

19  . 

Average  Voltage,  kV 

24.6 

Average  current,  kA 

126.6 
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TABLE  18 


ES/PC  SIZE  AND  CHARACTERISTICS  FOR  ORBIT  TRANSFER  MISSION 


Con  pone  at 

Par ame  ter 

ES/PC 
System 
for  1  PIT 

Capacitor 

Efficiency,  % 

98. 

Mass,  kg 

402.4 

Volxme,  m3 

0.22 

Dimension,  a 

0.61 

Capacitance,  (iF 

62.7 

Voltage,  kV 

24.6 

Energy  stored,  kJ 

19. 

Number  of  cycles 

9.1E7 

Specific  energy,  J/kg 

47.3 

Tbyra  tron 

Efficiency,  % 

98 .9 

Mass,  kg 

174.8 

Peak  power,  W 

1.16E9 

Pulse  Duration,  |is 

16.5 

Current,  kA 

126.6 

DC -DC 

Efficiency,  % 

92. 

Converter 

Mass,  kg 

204. 

Average  power,  kf 

225. 

Volume,  m3 

7.5E-3 

Dimension,  a 

0.42 

Radiator 

Temperature,  Deg.  K 

300. 

Mass,  kg 

432.5 

Surface  area,  a1 

134  .1 

Avg.  power  rejected,  kW 

24.3 
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SECTION  5 


CONCLUSIONS  AND  RECOMMENDATIONS 


S.l  CONCLUSIONS 

A  amber  of  specific  conclusions  can  be  drawn  from  this  study. 

(1)  Optimization  of  integrated  propulsion  systems  requires  that  the  primary 
power  sources,  the  ES/PC,  and  the  thruster  performance  be  related  to 
the  primary  propulsion  system  outputs  (pulse  mass,  exhaust  velocity, 
frequency,  and  total  thrusting  time). 

(2)  Models  which  meet  these  requirements  were  successfully  developed  for 
all  the  systems  of  interest. 

(3)  Most  of  the  models  developed  were  based  on  an  understanding  of  the 
physical  principles  of  operation  and  a  amber  of  assumptions  concerning 
the  performance  of  optimized  systems.  The  rest  of  the  models  were 
developed  by  curve-fitting  available  perfomance  data,  either  from 
manufacturer's  literature  or  from  point  design  studies. 

(4)  The  excellent  qualitative  and  quantitative  agreements  between  model 
prediction  and  observed  performance  lend  considerable  credence  to  the 
model  s. 

(5)  The  models  are  extremely  versatile,  permitting  rapid  evaluation  of 
proposed  alternatives  fox  specific  missions. 

(6)  The  models  can  easily  be  updated  when  new  data  becomes  available. 

Using  these  models  to  analyze  optimized  propulsion  systems  for  sta¬ 
tionkeeping  and  orbit  transfer  missions,  we  can  conclude; 

(1)  There  exists  one  set  of  propulsion  system  outputs  that  minimizes  the 
mass  of  the  system. 

(2)  The  optimum  propulsion  system  outputs  are  within  the  capabilities  of 
the  state-of-the-art  thrusters. 

(3)  The  propulsion  system  mass  is  about  10  percent  of  the  payload  for 
stationkeeping  and  ranges  from  2.5  to  3.5  times  the  payload  mass  for 
orbit  transfer  missions. 

(4)  Photovoltaic  power  sources  are  the  least  massive  for  stationkeeping 
missions  and  nuclear-based  power  sources  become  competitive  for  orbit 
transfer  missions. 

(5)  Average  power  levels  range  from  2  to  35  kW  for  stationkeeping  and  from 
200  kW  to  2  MV  for  orbit  transfer  missions. 

(6)  Capacitive  based  ES/PC  systems  are  the  least  massive,  inductive  based 
systems  are  comparable,  and  inertial  systems  are  much  heavier. 
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(7)  ES/PC  components  have  equal  weights  for  stationkeeping  systems,  and  the 
radiator  is  the  largest  component  in  the  system  for  orbit  transfer 
missions  (45%). 

5.2  RECOMMENDATIONS 

The  thruster  performance  plays  an  important  role  in  sizing  the  ES/PC 

and  the  primary  power  source.  It  is  necessary  to  validate  the  performance 

models  of  the  thrusters  over  the  wide  range  m  and  v  of  interest.  There 

P  c 

are  two  major  areas  of  uncertainty  in  the  thruster  models  developed.  The 
first  problem  arises  f rca  the  nature  of  the  thrusters  for  which  data  was 
available.  It  is  not  clear  that  the  experimental  thruster  perfomance  data 
available  to  us  was  generated  with  optimon  thrusters  (that  is,  the  thruster 
was  optimized  for  every  m^  and  ve  for  which  data  was  reported).  The  per— 
formance  data  must  be  evaluated  more  carefully  to  answer  this  question.  The 
second  problem  arises  from  the  limited  amount  of  data  available.  The 
missions  considered  in  this  study  require  the  thruster  models  apply  at 
operating  conditions  (mp  and  ve)  which  could  be  orders  different  from  the 
conditions  for  which  we  have  data.  Therefore,  we  recommend  that  performance 
data  for  optimum  thrusters  be  generated  over  a  much  mare  extensive  range  of 
®P  (or  mdotp)  and  ve. 

The  effects  of  the  component  life  (nimiber  of  cycles)  on  most  of  the 
components  forming  the  ES/PC  are  not  well  understood.  The  performance  of 
the  components  was  degraded  to  account  for  the  component  life  effects.  We 
recommend  that  a  more  extensive  analysis  be  conducted  to  evaluate  the  life 
effects  on  the  components  performance  of  the  most  promising  ES/PC.  This  can 
be  done  by  first  selecting  the  most  premising  ES/PC  Our  analysis  shear  that 
the  capacitive— based  ES/PC  is  the  most  promising.  Then  each  component  in 
the  ES/PC  should  be  tested  individually  to  study  the  life  effects  on  the 
performance  of  the  component.  After  testing  the  components  individually, 
they  should  be  integrated  and  tested  as  a  system  to  insure  that  no  com¬ 
pounded  effect  is  available. 
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APPENDIX  A 


PRIMARY  POWER  FOR  ELECTROMAGNETIC  PROPULSION 

Primary  power  source  performance  and  mass  models  are  required  to 
select,  size,  and  design  an  integrated  electromagnetic  propulsion  system.  A 
primary  pover  source  delivers  electrical  output  at  a  relatively  low  average 
level  to  an  Energy  Star age/Pover  Conditioning  system  (ES/PC)  which  condi¬ 
tions  the  pover  and  delivers  high-pover  pulses  to  the  electrical  thruster. 
The  primary  pover  source  is  defined  as  the  'raw'  pover  source  and  all  the 
components  (radiator,  generator,  converter,  etc.)  required  to  deliver 
electrical  pover  to  the  ES/PC 

The  performance  characteristics  and  mass  models  for  the  primary  pover 
sources  vere  determined  by  conducting  a  literature  survey  to  gather  the 
necessary  data,  contacting  technical  experts  to  confirm  the  results  and 
making  refinements  vhere  needed.  The  literature  contains  considerable  point 
design  data  but  very  little  perimetric  or  limit  data. 

A  vide  range  of  primary  power  options  vere  considered  including:  solar 
cell  arrays,  solar  thermophotovol  taie,  nuclear  thermionic,  nuclear  thermo¬ 
electric,  nuclear  Rankine  cycle,  nuclear  Brayton  cycle,  and  nuclear  MHD. 
Parametric  performance  models  vere  developed  for  pover  levels  ranging  from  a 
few  kilowatts  to  10  MW. 

A.1  SOLAR  CELL  ARRAYS 

Solar  cell  arrays  are  designed  in  collapsible  sections  and  are  opti¬ 
mized  to  minimize  power  transmission  losses.  Transmission  losses  are  ohaic 
losses  and  are  minimized  by  pover  transmission  at  the  highest  possible 
voltage.  Olaic  losses  are  directly  proportional  to  the  current  squared  for 
a  fixed  bus  bar  cross-sectional  area.* 

The  performance  of  the  solar  cell  arrays  is  affected  by:  (1)  plasma 
arcing,  (2)  voltage  breakdown,  (3)  voltage  transients  and  Van  Allen  Belts. 


^  Bauer,  David  P.,  Barber,  John  P.,  Swift,  Hal  lock  F.,  Vahlberg,  C  Julian, 
'Electric  Rail  Gun  Propulsion  Study  (Advanced  Electric  Technology-High 
Thrust), '  AFRPL-TR-81-02,  1981. 


Plasma  arcing  is  a  f  met  ion  of  ear  iromeat  plasma  density  (altitude) 
aad  tie  voltage  of  tie  array.  Gilmour2  listed  same  of  tie  laboratory  arciag 
data  obtaiaed  by  Steveas  at  NAS  A- Levi  s.  lie  data  sbow  that  arciag  caa  be 
expected  ia  GEO  witi  array  voltage  around  1000  V,  aad  caa  be  expected  ia  LEO 
at  around  300  V.  Voltage  breakdown  of  gases  may  occur  if  tie  operating 
voltage  is  iigier  tiaa  tie  Pa scien  minimum.  Gilmoux  presented  Pascien 
ciaracteristics  for  some  molecular  gases  (see  Figure  A-l).  Tiis  figure 
siows  tiat  for  low  pressure  (iigi  vacuum)  1000  V  may  be  feasible.  Voltage 
transients  are  experienced  immediately  after  eclipse.  Wien  tie  cold  solar 
array  is  suddenly  exposed  to  sunl  igit,  tie  voltage  (for  a  given  current)  is 
2  to  2.5  times  iigier  tiaa  tie  warm  array  voltage,  liese  voltage  transients 
can  be  overcome  if  provisions  are  made  in  tie  design.  Some  of  tie  sug¬ 
gestions  to  overcame  tiis  problem  are:  (1)  to  connect  tie  solar  array  to 
tie  power  system  only  after  it  warms  up,  or  (2)  to  clamp  tie  solar  array 
voltage  to  a  battery  voltage  and  make  provisions  to  accommodate  tie  current 
pulse  from  tie  cold  array. 


2 

Figure  A-l.  Pascien  Ciaracteristics  For  Same  Molecular  Oises. 


2 

Gilmour,  A.  S.  Jr.,  'Investigation  of  Power  Processing  Tecinology  for 
Spacecraft  Applications,'  AFWAL-TE-82-2054,  June  19  82. 
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Solar  cell  arrays  degrade  when  passing  through  the  Van  Allen  belts.  Degra¬ 
dation  of  up  to  40%  is  reported.  3  Provisions  Bust  be  made  for  this  degra¬ 
dation  if  solar  cell  arrays  axe  to  be  used  in  orbit  transfer  missions. 

Solar  cell  arrays  can  be  designed  to  deliver  any  power  level  required 
since  they  are  composed  of  solar  cells  connected  in  series  and  parallel. 
Solar  cell  arrays  capable  of  providing  power  levels  up  to  10  kff  have  been 
deployed  in  space  (used  on  Skylab).  Arrays  capable  of  multi-hundred  kilo¬ 
watt  have  been  designed  and  analyzed. 

Operating  voltages  ranging  from  28  V  to  900  V  have  been  reported. 
References  4.  5,  6,  7,  and  8  report  the  performance  of  arrays  delivering 
multi-hundred  watts  with  operating  voltage  ranging  from  28  to  42  VDC.  Ref¬ 
erences  7  and  9  report  analyses  of  solar  power  systems  ranging  from  100  to 
300  V.  References  10  and  11  reported  tradeoff  analyses  for  power  management 
and  transmission  of  multi-hundred  kilowatt  and  multi-megawatt  systems.  The 
multi-hundred  kilowatt  system  had  an  operating  voltage  of  750  for  the  DC 
transmission  and  an  operating  voltage  of  440  for  an  AC/DC  transmission. 


^Taussig,  R. ,  et  al..  'Overview  Study  of  Space  Power  Technologies  for  the 
Advanced  Energetics  Program.'  NASA  CRrl65269,  1981. 

^Brooks,  G.  R.,  'Orbiting  Solar  Observatory  (0S0-8)  Solar  Panel  Design  and 
in  Orbit  Performance.'  Intersociety  Energy  Conversion  Engineering  Conference 
(IE CEC),  789321.  1978. 

^Thronton,  'ATM  Solar  Amy  In-Flight  Performance  Analysis,'  IE  CEC,  1974. 

**Gol  dhammer ,  L.  J. ,  'Design  and  Flight  Performance  of  the  Pioneer  Venus 
Multiprobe  and  Orbiter  Solar  Arr ays.'IECBC,  809076,  1980. 

7 

Luke  ns,  F.  E.,  'Advanced  Development  of  a  Programmable  Power  Processor,' 
IE  CEC,  809076,  19  80. 

8McEinney,  H.  N. ,  Briggs,  D.  C.,  'Electrical  Power  Subsystem  for  the 
Intelsat  V  Satellite.'  IECEC,  789084,  1978. 

9Glass,  M.  C. ,  'A  Six  Kilowatt  Transformer-Coupled  Converter  for  Space 
Shuttle  Solar  Power  Systems,'  IECEC,  809146,  1980. 

10Mildice,  J.  W. ,  'Study  of  Power  Management  Technology  for  Orbital  Mul  ti- 
100  kff  Applications,'  NASA  CR159834,  July  1980. 

^Peterson,  D.  M.  and  Pleasant,  R.  L.,  'Study  of  Multimegawatt  Technology 
Needs  for  Photovoltaic  Space  Power  Systems,'  GDC- AST- 81-019,  August  1981. 
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Hie  mul  ti-megawatt  system  had  an  operating  voltage  of  900  for  AC  and  DC 
The  mul  ti-transaission  megawatt  systee  will  require  multiple  shuttle 
launches  and  construction  in  space. 

System  mass  data  gathered  from  references  1-18  are  presented  in 
Figure  A-2.  Figure  A- 2  shows  that  the  specific  mass  is  about  SO  kg/kW  for 
one  kilowatt  arrays  and  decreases  asymptotically  to  10  kg/kV  for  100  kV  and 
above.  This  data  was  curve— fitted  to  generate  a  mass  model: 

«,p  -  0.013  +  40.6 1/Ppp  (A—  1 ) 

where  mjp  is  the  specific  mass  in  kg/W  and  Ppp  i*  the  average  power  level  in 
watts.  Equation  A-l  is  true  for  power  levels  above  1000  V.  For  power 
levels  lower  than  1000  W,  the  mass  of  the  system  is  assumed  to  be  54  kg 
independent  of  power  level.  An  overall  degradation  of  20%  in  the  perfor¬ 
mance  of  the  arrays  was  assumed  for  orbit  transfer  mission. 

A.2  SOLAR  THERMO  PHOTOVOLTAIC 

The  thermophotovol  taic  (TPV)  concept  offers  a  viable  alternative  to 
conventional  solar  cell  arrays.  It  offers  a  high  efficiency  (50%),  19  thus 

12Briggs,  D.  C,  Pollard,  H.  E.,  'Intelsat  V  Solar  Array  Design  and  Develop¬ 
ment  Summary.'  IECEC.  789322,  1978. 

13Kelly,  F.  G. ,  Luf  t,  W. ,  Kurland,  R.  M. ,  'Design  Features  of  the  TSDSS 
Solar  Array,'  IECEC,  78  9323,  1978. 

1 4 

Borduas,  H.  F.,  et  al.,  'Medium  Power  Deployable  Hybrid  Solar  Arrays,' 
IECEC,  780328.  1978. 

^Peterson,  D.  G. ,  'INSAT-1  Solar  Array  Design  and  Development  Summary,' 
IECEC,  19  80. 

^Souza,  C.  J.,  'Large  Area  Flexible  Solar  Array  Design  for  Space  Shuttle 
Application,'  IECEC,  809485,  1980. 

^Woodcock,  W.  G.  Ill,  'Multi-Hundred  Kilowatt  Solar  Arrays  for  Space,' 
IECEC,  809271.  19  80. 

^Trumble,  T.  M. ,  'Space  Application  for  Gallium  Arsenide  Solar  Cells,' 
IECEC,  82  9  269,  19  82. 

1  9 

Horne,  W.  E. ,  et  al.,  'Solar  Thermophotovol  taic  Space  Power  System,' 
IECEC.  809072.  19  80. 
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Figure  A-2.  Specific  Mass  of  Solar  Cell  Array  Systems. 

a  reduction  in  photovoltaic  cell  area  and  a  possible  reduction  in  weight. 
The  photovoltaic  cells  axe  shielded  from  the  external  environment  and  do  not 
suffer  radiation  degradation  to  the  extent  that  conventional  panels  do. 

TPV  power  systems  can  be  designed  to  deliver  any  power  level  required 
and  the  same  limiting  factors  which  apply  to  solar  cell  arrays  (voltage 
levels,  breakdown,  and  transients)  apply  to  TPV  systems. 

TPV  systems  are  still  in  their  infancy  and  many  parameters  which  affect 
performance  are  still  not  well  understood.  Design  criteria  and  methods  have 
not  been  developed.  A  literature  survey  revealed  that  feasibility  and 
sensitivity  studies  have  been  conducted.1^ *20,21  Very  little  specific  mass 
data  was  located.  Reference  19  reported  seme  specific  mass  data,  which  we 
used  in  conjunction  with  a  system  energy  balance  model  (Figure  A-3), 


20Swanson,  R.  M.,  et  el.,  'Silicon  Photovoltaic  Cells  in  Thexaophotov  ol  taic 
Conversion,'  EPI-ER-478  Progress  Report,  Stanford  Electronics  Lab.,  February 
1977. 

2 1 

Horne,  W.  E.,  et  al.,  'Improved  Thermophotovol  taic  Power  Systems,'  IECEC, 
82  9018,  1982. 
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reported  in  reference  20,  to  develop  a  specific  mess  model.  The  specific 
mess  for  e  single  unit  is  presented  in  Figure  A-4.  Figure  A- 4  indicates 
that  a  panel  made  with  cl  asters  of  small  TPV  converters  would  be  lighter 
than  a  single  large  concentrator  doe  to  the  increased  performance  of  the 
heat  pipes  (in  the  converter)  as  heat  pipe  length  is  decreased. 


Concentrator  Radius  (ft) 


Specific  Mass  Versos  Concentrator  Sadias  Of  Th  eraopho  t  ov  ol 
Sy  s  t  em  s . 


Figure  A-4. 


The  assumptions  made  in  order  to  model  the  specific  mess  for  large 
arrays  were: 

(1)  72%  of  the  incident  solar  energy  is  retained  in  the  TPV  converter. 

(2)  The  efficiency  of  the  photovoltaic  cells  is  40%. 

(3)  One  concentrator  is  used  if  the  required  power  level  can  be 
achieved  with  a  concentrator  radius  less  than  3m.  (limitations  of 
the  available  data),  and  multi-3m.  radius  concentrators  are  used 
for  higher  power  levels.  i.e.,  the  specific  mass  is  constant  for 
power  levels  requiring  more  than  3m.  radius  concentrators. 

We  assumed  that  large  TPV  panels  will  be  formed  with  large  concen¬ 
trators,  eventhough  smaller  concentrators  have  a  lower  specific 
mass,  because  we  felt  that  their  advantages  outweigh  the  possible 
savings  in  mass  of  the  TPV  system.  The  advantages  of  the  large 
concentrators  are:  they  minimize  the  number  of  concentrators 
required  to  achieve  a  power  level  (leads  to  less  complex  system), 
they  reduce  manufacturing  slope  errors  of  concentrator  surfaces, 
and  they  simplify  the  sun- TPV  alignment  requirement. 

The  specific  mass  model  for  power  levels  ranges  from  1  kW  to  11  kW  is: 


»Sp  *  0.005  +  4.6xlO-7Ppp  (A-2) 

where  mJp  is  in  kg/W  and  Ppp  i*  the  average  power  level  in  watts.  For 
higher  power  levels  the  specific  mass  was  assumed  constant  and  equal  to  10.1 
kg/kw.  The  complete  model  is  illustrated  in  Figure  A-5. 


A. 3  NUCLEAR  THERMOELECTRIC 

Nuclear  primary  power  for  electromagnetic  propulsion  requires 
thermal-to-ele ctrical  converters  capable  of  operating  for  long  periods  of 
time  at  high  temperatures.  Thermoelectric  converters  are  in  a  unique  posi¬ 
tion  to  meet  these  requirements.  They  are  a  direct  energy  conversion 
system.  Static  conversion  and  modularity  are  attractive  characteristics. 
A  large  number  of  parallel  converters  provide  a  system  with  high  redundancy 
and  reliability.  Modular  design  permits  flexibility  in  design  and 
operation. 

Mission  studies  have  shown  that  a  400  kW  nuclear  thermoelectric  power 
system  can  be  deployed  in  a  single  shuttle  launch  with  the  largest  payload 
of  interest  to  this  study.  Thermoelectric  power  systems  capable  of  120  kW 
have  been  analyzed  by  JPL  and  LASL  with  the  400  kW  being  their  ultimate 
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Figure  A-5 .  Specific  Mass  Of  TPV  Systems. 

range.  Pover  systems  for  higher  power  levels  can  be  designed  bat  will 
require  mal  tiple  shuttle  launches. 

Output  voltage  levels  of  up  to  150  volts  have  been  proposed  for  120  kW 
thermoelectric  power  sy  sterns.^  *^3  He  voltage  level  is  limited  by  the 
electrical  insulating  material.  The  high  operating  temperature  degrades  the 
insulating  material.  Bauer  in  reference  1  reported  a  personal  ccmmuni  ca  tion 
with  Dr.  Ranken  (LASL)  who  projected  that  developments  in  insulating  mate¬ 
rials  will  permit  three  times  the  presently  projected  voltage  within  a 
de  ca  de . 


22 

Mondt,  J.  F.,  et  al.,  'Nuclear  Power  Source  for  Electric  Propulsion,'  14th 
International  Electric  Propulsion  Conference,  1979. 

23 

■‘•'Ranken,  W.  A.,  'Experimental  Results  for  Space  Nuclear  Power  Plant  De¬ 
sign,'  IECEC,  809142,  1980. 


A  specific  mass  node  1  was  generated  from  data  in  references  22  thru  32. 
Hie  referenced  data  was  for  different  reactors,  radiators,  and  operating 
strategies.  The  type  of  reactor  and  radiator  has  a  strong  impact  on  the 
mass  of  the  system.  The  specific  mass  model  was  generated  by  curve-f itting 
the  conservative  data.  The  specific  mass,  for  power  levels  ranging  f rcm 
1  kW  to  150  kV  (limitation  of  the  available  data)  can  be  expressed  as: 

®#p  -  0.021  +  370.3/Ppp  (A-3) 

where  mjp  is  in  kg /W  and  ppp  in  W.  Ve  will  assume  that  the  specific  mass 
model  (Equation  A-3)  is  valid  for  power  levels  above  the  range  of  the 
available  data.  Equation  A-3  and  the  referenced  data  points  are  presented 
in  Figure  A-6. 


24Kelly,  C.  E. ,  Ambrose,  G.  R. ,  'Testing  of  the  GPHS  Electrically  Heated 
Thermoelectric  Converter,'  IECEC,  82  9234,  1982. 

2  5 

Eisner,  N.  B. ,  'Development  of  An  Advanced  RIG  Using  Segmented  Thermoelec¬ 
trics,'  IECEC,  1980. 

26Koeing,  D.  R. ,  'Heat  Pipe  Reactors  for  Space  Power  Applications,'  AIAA-77- 
491,  1977. 

27 

Ewell,  R.  C.,  Stapfer,  G. ,  'Thermoelectric  Conversion  for  Space  Nuclear 
Power  Systems,'  IECEC,  82  9  238.  19  82. 

2  ft 

Hockey,  D.  E. ,  et  al.,  'Comparison  of  Evolving  Photovoltaic  and  Nuclear 
power  Systems,'  IECEC,  82  9011,  1982. 

^Mahefkey,  T.,  'Future  Space  Power — The  D.O.D.  Perspective,'  IECEC,  809016, 
'980. 

3  0 

Ranken,  V.  A.,  'Experimental  Results  for  Space  Nuclear  Power  Plant 
Design,'  IECEC.  809142,  1980. 

^Buden.  D. ,  'Space  Nuclear  Reactor  Power  Plants,'  Los  Alamos,  LA8823-MS, 
19  80. 

32Buden,  D.,  'Reactor  Technology  for  Space  Electrical  Power  (10  -  100  kWe)/ 
LA6891-MS,  1977. 
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Figure  A-6.  Specific  Muss  Of  Nuclear  Thermoelectric  Power  Systems. 

A.  4  NUCLEAR -THERMIONIC 

Thermionic  converters  are  also  good  candidates  to  meet  the  nuclear 
electromagnetic  propulsion  needs.  They  are  also  direct  energy  conversion 
systems  and  are  modular. 

Thermionic  power  systems  can  be  designed  for  the  power  level  range  of 
interest.  Power  systems  above  400  kW  will  require  multiple  shuttle  launch¬ 
es.  Thermionic  converters  offer  high  current  density  (10  A/cm^)  and  low 
module  voltage.  The  maxima  system  voltage  which  has  been  projected  is 
50  V.33  the  vol  tage  i  s  limited  by  the  strength  of  the  high  temperature 
electrical  insulation. 

Oats  from  references  26,  29,  32,  34,  35,  and  36  were  used  to  model  the 


■^Phil  1  ips,  W. ,  (JHJ,  Personal  Communication,  July  1982. 

3  4 

Dick,  R.  S.  Jr.,  'Collector  Temperature  Effects  on  the  Performance  of 
Advanced  Thermionic  Converters  and  Nuclear  Electric  Propulsion  Systems,' 
EECEC,  809359,  1980. 

^Phillips,  W.  M. ,  'Nuclear  Electric  Power  System  for  Solar  System  Explor¬ 
ation,'  AIAA  79  13374. 
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current 

average  current 

maximum  current 
axial  inductance  gradient 
resi stance 
vol tage 

voltage  required  for  acceleration 
average  voltage 
breakdown  voltage 

voltage  characteristic  of  energy  dissipated 
'electrode  fall  voltage' 
ionization  voltage 

Z  acceleration  length 

Z  effective  plasma  length 
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Figure  A-9 .  Specific  Muss  Of  Nuclear  Rankine  Cycle  Power  Systems. 

the  magnetic  field  gives  rise  to  an  electromotive  force  end  censes  e  trans— 
verse  current  to  flew  in  the  moving  ges.  The  current  is  collected  by  the 
electrodes  end  carried  to  the  external  load.  References  41  end  42  discuss 
MHO  generation  at  power  levels  (10  -  200  MW)  orders  of  magnitude  higher  than 
that  required  for  propulsion.  ‘The  only  pertinent  data  (at  the  upper  level 
of  the  power  range  of  interest)  was  found  in  reference  43.  This  data  was 
used  to  model  the  specific  mass.  The  specific  mass  for  power  levels  above 
100  kW  is: 


■Sp  -  10"4  +  1990 J./Ppp 


(A- 8) 


where  m$^  is  in  kg/W  and  Ppp  in  W. 


The  results  are  shown  in  Figure  A-10. 


41T.  Botts,  et  al.f  'A  Compact  High  Performance  Electrical  Power  Source 
Based  Upon  the  Rotating  Bed  Nuclear  Reactor,'  Brookhaven  National  Labora¬ 
tory,  January  1980. 

42 

Rosa,  R.  J. ,  'MHD  Applied  to  Aerospace  and  Deep  Space  Propulsion,'  AIAA 
82-1212,  June  1982. 

43 

Botts,  T.  E.,  et  al.,  'Nuclear  Reactors  Using  Fine  Particulate  Fuel  for 
Primary  Power  in  Space,'  Brookhaven  National  Laboratory,  August  19  82. 
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Figure  A-8.  Specific  Mass  of  Nuclear  Brayton  Cycle  Power  Systems. 

No  quantitative  specific  mass  data  of  the  nuclear  Rankine  system  was 
found  in  the  literature.  The  mass  of  a  nuclear  Rankine  system  should  be 
similar  to  the  mass  of  the  Brayton  system.  The  Rankine  cycle  operates  at 
higher  temperature  than  Brayton  cycle  but  there  are  no  other  significant 
differences.  We  assumed  that  the  lower  bound  of  the  data  for  Brayton  system 
is  the  mean  for  Rankine  system. 

The  specific  mass  is  then  expressed  as: 

■ip  *  0.013  +  178.6/Ppp  (A-7) 

and  is  illustrated  in  Figure  A-9 . 

A. 6  NUCLEAR  MHD 

An  MHD  energy  converter  consists  of  a  gas  dynamic  nozzle  (expander; 
followed  by  a  duct.  The  duct  is  surrounded  by  a  coil  which  produces  a 
magnetic  field  in  the  duct.  Electrodes  are  placed  at  the  top  and  bottom  of 
the  duct.  Ionized  gas  is  expanded  in  the  nozzle  and  converts  some  of  its 
pressure  and  internal  energy  into  directed  motion.  The  gas  motion  through 
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A. 5  NC CLEAR  BRAYTON/ RANKINE 


The  Brayton  and  Rani,  me  cycles  used  in  tnxbo-generator  s  have  l  cm  spe¬ 
cific  mass.  They  are  designed  as  modular  units  and  are  connected  m  para, 
lei  for  redundancy  and  increased  system  reliability.  Each  unit  is  usually 
designed  to  deliver  a  fraction  of  the  total  power  required  with  the  capabil¬ 
ity  of  delivering  much  longer  power  in  case  of  module  failure.  The  input 
pressure  level  to  either  cycle  is  approximately  proportional  to  the  output 
power  level  required. 

Mission  studies  requiring  power  levels  ranging  from  25  to  600  kW  have 
been  conducted  with  the  nuclear  Brayton  cy  cl  e.^  *38 ,39 ,40  He  power  was 
delivered  at  120  VDC  or  at  67  VAC,  three  phase,  1733  Hz.  Reference  37  is  a 
study  of  a  Brayton  cycle  delivering  400  kW,  500  V,  single  phase,  3000  Hz. 
The  voltage  in  all  cases  is  limited  by  the  alternator.  As  alternator 
technology  improves,  higher  voltages  and  frequencies  may  be  feasible. 

Data  for  the  nuclear  Brayton  power  system  from  references  32,  35,  37, 
38,  39,  and  40  was  used  to  model  the  specific  mass.  The  specific  mass  model 

is : 


®sp  =  0.019  +  190 .6/Ppp 

for  power  levels  above  1.5  kW  and  is  presented  in  Figure  A-8. 


(A-6) 
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specific  mass. 


The  specific  mass  is  expressed  as: 


I  “sp  =  0.018  +  514/Ppp  (A-5) 

for  power  levels  above  10  kW  (the  lower  limit  of  the  available  data).  We 
will  assume  that  this  model  is  valid  beyond  the  range  of  the  available  data. 
£  This  model  and  the  data  are  presented  in  Figure  A-7.  The  data  presented  in 

Figure  A- 7  include  system  design  and  analysis  results  for  different  reactors 
(heat  pipe,  heat  pipe  solid  core,  solid  core  fluidized  bed)  and  for  differ¬ 
ent  radiators  (conventional,  advanced,  and  open  loop).  The  type  of  reactor 
|  and  radiator  selected  has  a  considerable  effect  on  the  mass.  The  model 

presented  in  Figure  A-7  and  Equation  A-5  represents  an  average  of  all  the 
da  ta. 


Figure  A-7.  Specific  Mass  Of  Nuclear  Thermionic  Power  Systems. 
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•’^Phillips,  W.  M.,  Koening,  D.  R. ,  'Reactor  Design  and  Integration  into  a 
Nuclear  Electric  Spacecraft,'  IECEC,  789564,  1978. 


APPENDIX  B 

PERFORMANCE  MODELING  APPROACH 


B.l  INTRODUCTION 

The  thruster  input  pulse  parameters  sre  (1)  pulse  energy,  Ep,  (2)  cur¬ 
rent,  J,  (3)  voltage,  V,  end  (4)  pulse  duration,  tp.  The se  input  parameters 
are  related  to  the  thruster  output  parameters  and  many  other  thruster  char¬ 
acteristics.  'Die  thruster  output  parameters  which  affect  the  input  require¬ 
ments  are  pulse  mass  and  exhaust  velocity.  The  other  output  parameters, 
pulse  frequency  and  thruster  lifetime,  do  not  directly  affect  thruster  input 
pulse  requirements.  Other  thruster  characteristics  (such  as  size,  electrode 
geometry  and  material,  propellant  material,  etc.)  clearly  affect  thruster 
performance  and  therefore  must  be  considered.  These  other  characteristics 
cannot  be  selected  arbitrarily  since  they  do  affect  thruster  performance 
(i.e.,  thruster  input  requirements).  In  principle,  they  would  be  deter¬ 
mined  by  optimizing  the  design  of  the  thruster. 

The  efficiency  of  the  thruster  has  a  direct  effect  on  the  power  and 
energy  requirements,  and  therefore  the  size  of  the  primary  power  supply  and 
the  ES/PC  system.  Therefore,  the  optimization  criterion  for  the  thruster  is 
the  maximization  of  efficiency.  Thruster  characteristics  must  be  chosen  to 
maximize  efficiency.  Assuming  this  can  be  accomplished,  the  input  pulse 
parameters  for  optimum  thrusters  are  functions  of  only  the  pulse  mass  and 
exhaust  velocity.  (That  is,  for  each  combination  of  pulse  mass  and  exhaust 
velocity  there  is  only  one  thruster  configuration  which  has  maximum  effi¬ 
ciency  and  it  is  this  one  which  is  of  interest  here.) 

Let  us  assiHe  that  a  large  body  of  data  existed  (either  experimental  or 
analytical)  for  optimized  thrusters  and  that  this  data  contained  values  for 
the  thruster  input  pulse  parameters  over  wide  ranges  of  mp  and  ve.  A 
simple,  anil ti-variate  curve  fitting  technique  could  then  be  employed  to 
develop  the  required  relationships  for  Ep,  J,  V,  and  tp  as  functions  of  mp 
and  vQ.  This  is  the  approach  which  was  used  for  modeling  the  electric  rail 
gun.*  Table  B-l  lists  the  thruster  characteristics  which  were  chosen  opti¬ 
mally,  leaving  only  mp  *nd  ve  as  independent  thruster  output  parameters. 

*Bauer,  D.  P. ,  Barber,  J.  P. .  Vahlberg,  C.  J. ,  "Die  Electric  Rail  Gun  for 
Space  Propulsion,'  NASA-CR-165312 ,  Final  Report,  February  1981. 


TABLE  B-l .  PARAMETERS  AFFECTING  RAIL  GUN  EFFICIENCY 


Parameter 

Increasing  Effi¬ 
ciency  Results 

From  Increasing/ 
Decreasing  Hag- 
nitude  of  Para¬ 
meter 

Optimizing 
Value  or 
Physically 
Imposed 
Constraint 

Physical 

Implications 

Value  of 

Parameters 

USED  OR 

Derived 

Rail  Gun 
|  Inductance 
lpE»  Uni: 

Length 

Increasing 

Optimum 

Rail  guns  with  souare 
bore  are  optimal.* 

L'  •  0.63  mH/m 

i 

Irail  Gun 

ILEnGTh 

Increasing 

Optimum 

Practical  optimum 
length,  x,  is  a  func¬ 
tion  OF  bore  size  AND 

EXHAUST  VELOCITY  for 
GIVEN  VALUES  of  other 
parameters.** 

i  *  AS. 6  X  IO^hv2 

|kail  Tem¬ 
perature 

1 

Decreasing 

Constrained 

.Tradeoffs  required, 

TO  OBTAIN  PRACTICAL 

WASTE  HEAT  RADIATOR 

SIZE  AND  ADEQUATE 
STRENGTH  RAILS  WITH 
ACCEPTABLE  resistiv¬ 
ity. 

T  -  A50'K 

Rail  Gun 

Bore  Size 

"Increasing 

Constrained 

Specified  by  mission 
requirements. 

0.1  CM  <  H  <  1.0  CM 

Pellet  £*- 
IHAUST  VEL¬ 
OCITY 

Decreasing 

Constrained 

Specified  by  mission 
requirements.  Bounds 
were  imposed. 

5  KM/S  <  V  <  20  XM/S 

I'HELLET  Ac- 
iSE^ERAT  I  NG 

IStress 

1 

Increasing 

Constrained 

Accelerating  stress 
limited  to  maximum 

ALLOWABLE  STRESS  IN 
PELLET. *** 

•  Q  *  200  r,N/M2 

pellet  areal 
Density 

Decreasing 

Constrained 

Both  pellet  density, 

p,  AND  PELLET  LENGTH, 

D,  SHOULD  BE  MINI¬ 
MIZED.  Practical 
pellet  density  chosen. 
Pellet  length  in  di¬ 
rection  OF  TRAVEL 

MUST  BE  LONGER  THAN 
ONE-HALF  THE  bore 

SIZE,  TO  PREVENT  IN¬ 
BORE  PELLET  TUMBLING.* 

p  »  3000  kg/m^ 

D  *  h/2 

'The  combination  of  a  square  bore  and  pellet  length  equal  to  one-half  the  bore  size  allows  the  pellet  mass 
to  be  related  directly  to  bore  size  (i.e.,  n  »  rh5/2  ■  I50Dh5).  Given  the  pellet  density,  pellet  mass  and 

BORE  SIZE  MAY  BE  USED  INTERCHANGEABLY. 

"A  PRACTICAL  OPTIMUM  CHARACTERISTIC  PELLET  ACCELERATION  TIME  IS  THEREFORE  IMPLIED  AND  WAS  FOUND  TO  COR¬ 
RELATE  TO  RAIL  HEIGHT  AND  EXHAUST  VELOCITY  AS,  *  5 1  X  lO'^HV. 

HE  COMBINATION  OF  A  DEFINED  OPTIMAL  RAIL  INDUCTANCE.  I',  AND  A  DEFINED  ALLOWABLE  STRESS,  oQ,  DEFINES  THE 
-\-f  RAIL  GUN  CURRENT,  JQ,  AS  GIVEN  BY,  J  *  ( 2 Oq/L ' ) * ^ H  *  2.5?  X  lO^H. 

;  ’  R  Q  2  1 
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A  rail  gun  simulation  computer  program  (containing  these  optimum 
characteristics)  was  used  to  generate  a  large  data  base  of  thruster  input 
pulse  requirements  over  wide  ranges  of  m  and  v  .  The  data  was  then  reduced 

P  C 

using  regression  techniques  to  produce  a  set  of  parametric  equations 
describing  optimum  ERG  performance.  The  equations  are: 


n 


13.67mO-10V2" 


(B-l) 


J  «  2.2xl06m3/3  ( B-2) 

P 

t  “  4.45xlO“6mJ/3ve  (B-3) 


This  simple  data-based  approach  cannot  be  applied  to  the  other  three 
thruster  types  (PIT.  TPP.  MPD)  because  not  enough  performance  data  (either 
experimental  or  analytical)  is  available  over  a  wide  enongh  range  of  and 

ve.  In  addition,  it  is  not  always  clear  that  existing  data  represents 
optimum  thruster  performance.  For  example,  data  over  considerable  ranges  of 
mp  and  ve  is  usually  generated  experimentally  with  one  thruster  configura¬ 
tion.  Intuitively,  one  would  expect  that  the  configuration  of  an  optimtn 
thruster  would  be  different  at  these  wide  extremes.  The  task  of  generating 
a  more  complete  body  of  performance  data  (either  experimentally  or  analyt¬ 
ically)  is  considerable  and  is  beyond  the  scope  of  this  report.  Methods  of 
performance  extrapolation  based  on  a  sound  physical  understanding  of  basic 
thruster  operation  must  therefore  be  employed. 

In  the  following  paragraphs  we  will  examine  each  of  the  thruster  input 
parameters  and  develop  methods  of  extrapolating  thruster  performance  based 
on  existing  performance  data  and  physically  reasonable  assumptions  concern¬ 
ing  thruster  operation. 


B.2  PULSE  ENERGY 


The  purpose 
independent  output 
with  an  assumption 


of  this  section  is  to  relate  the  pulse  energy 
parameters  of  exhaust  velocity  and  pulse  mass, 
concerning  the  components  of  pulse  energy. 


to  the 
We  begin 
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Assumption  No.  1  -  The  pulse  energy  is  the  sum  of  3 
components:  (1)  kinetic  energy,  (2)  'ionization' 

energy  end  (3)  dissipation  energy  end  nay  be 
expressed  as  : 


E.  .  +  E.  +  E.. 

kin  ion  dis 


(B-4) 


The  first  tern  in  the  pnlse  energy  equation,  the  kinetic  energy,  is  the 
uaefnl  energy  inparted  to  the  pnlse  ness  during  the  pulse  and  is  given  by: 


»_v«/2 
p  e 


( B-5 ) 


The  second  term,  the  'ionization*  energy,  includes  all  energy  that  is 
'frozen'  into  the  propellant  and  is  characteristic  of  the  state  of  the 
exhausted  gas.  The  term  is  not  strictly  ionization  but  includes  depolymeri¬ 
zation  (where  applicable),  ionization,  and  heating  of  the  gas.  This  term 
cannot  be  easily  evaluated  without  a  further  assumption. 


Asamaption  No.  2  -  The  state  of  the  plasma  exiting 
the  thruster  is  independent  of  exhaust  velocity  and 
thruster  size  (pulse  mass). 


With  this  assumption,  the  term  can  be  written  as: 


E  •  »  me 

ion  pi 


(B-6) 


where  e^  is  a  'specific  ionization  energy'  characteristic  of  the  state  of 
the  plasma  at  the  thruster  exit  and  has  units  of  energy/mass.  Caution 
should  be  exercised  in  interpreting  this  term  physically.  The  specific 
energy,  e^(  is  just  a  parameter  and  its  value  will  be  determined  from 
experimental  performance  data. 

The  third  term  in  the  pulse  energy  equation  is  s  dissipative  energy 
loss  term.  It  includes  all  the  energy  required  to  maintain  the  plasma  at 
temperature.  For  both  the  PIT  and  TPP,  the  energy  loss  mechanism  was  char¬ 
acterized  as  a  resistance.  The  term  contains  all  dissipative  energy  losses 
in  the  thruster,  including  electrode  losses  and  radiative  losses  from  the 
plasma.  The  dissipative  loss  term  for  the  PIT  and  TPP  can  be  written  as: 


Edis  -  /RJ2dt 


(B-7) 


where  &  is  the  apparent  'resistance'  of  the  thruster  and  J  is  the  pulse 
current. 

For  the  MPD  thruster,  the  dissipative  energy  loss  term  was  broken  into 
two  components;  a  resistive  loss  term  and  an  'electrode  fall'  loss  term. 
The  complete  dissipative  loss  term  can  be  written  as: 

Edis  =  /aj2dt  +  /VfJdt  (B"8> 

where  Vj  is  the  'electrode  fall  voltage'.  The  first  component  is  identical 
to  the  energy  loss  mechanism  for  the  TPP  and  PIT.  The  second  component  is 
an  energy  loss  attributed  to  'electrode  fall.’  The  voltage  drop,  Vf,  is  the 
voltage  penalty  that  must  be  paid  to  remove  current  from  the  electrodes  and 
is  assumed  to  be  constant. 

Thruster  performance  data  strongly  indicates  that  the  predominant  re¬ 
sistive  element  in  the  circuit  is  the  plasma.  Ve  conclude  that: 

Assumption  No.  3  -  The  resistive  dissipation  in  the 
thruster  is  primarily  due  to  plasma  resistance. 

We  further  note  that  ionization  must  occur  early  in  the  acceleration 
(neutrals  are  accelerated  only  poorly  at  best).  This  implies  that  the  final 
state  of  the  plasma  is  reached  in  a  small  fraction  of  the  pulse  time  and 
that  : 


Assumption  No.  4  -  The  apparent  dissipative  resis¬ 
tance  remains  constant  during  acceleration. 

Equation  B-7  can  then  be  rewritten  as: 

Edi s  -  R/J2dt  ( B-9) 

(It  should  be  noted  that  even  if  assumption  no.  4  is  not  strictly  obeyed, 
but  the  variation  of  resistance  scales  similarly  with  time  as  thruster 
output  is  varied,  then  Equation  B-9  will  still  be  valid.) 
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Equation  B-8  can  also  be  rewritten  as: 


rr 


Edis  ■  +  Vf/jdt  ( B-10) 

Experimental  thruster  performance  measurements  indicate  that  for  ex¬ 
haust  velocities  of  interest  0  10,000  ■/ s) ,  the  thrust  due  to  electrothermal 
effects  is  very  small  and: 

Assumption  No.  5  -  Hie  acceleration  of  the  plasma 
is  primarily  magnetic. 

The  electromagnetic  impulse  can  then  be  written: 

*  (L‘/2)  /  J2dt  (B-ll) 

y  c 

where  L'  is  the  axial  inductance  gradient  for  the  device.  The  inductance 
gradient  is  a  constant  which  depends  on  device  ge  one  try  but  does  not  vary 
with  time. 

The  dissipative  energy  loss  term  for  the  TPP  and  PIT  can  be  rewritten 
by  combining  equations  B-9  and  B-ll  thus  producing: 

Edis  *  ( 2R/L' ) Bpve  (B-12) 

For  the  MPD,  the  resulting  equation  is: 

Edis  -  (28/L'J  mpve  +  Vf(2ipVe/L' )1/2tp  (B-13) 

These  equations  imply  the  very  interesting  result  that  the  dissipated  energy 
depends  on  the  impulse  bit,  ■  v  . 

r  6 

Ve  must  now  examine  how  R  and  L'  vary  with  pulse  mass  and  exhaust 
velocity.  The  plana  resistance  should  be  independent  of  exhaust  velocity 
according  to  assumptions  2  and  3  and  will  depend  on  plssma  geometry  and 
size.  The  resistance  can  be  expressed  as: 

R  ■>  IV A  (B-14) 
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where  (  is  the  plasms  resistivity,  £  is  the  effective  plasma  length,  and  A 
is  the  cross-sectional  area  of  the  plasma  perpendicular  to  the  path  of 
current  flow. 

The  pulse  mass  for  the  PIT  and  TPP  is  related  to  the  plasma  dimensions 
and  plasma  density  by: 


mp  -  plk  ( B-15) 

wh*re  £  is  the  plasma  length  and  A  is  the  plasma  cross-sectional  area. 
Obviously,  the  plasma  resistance,  R,  is  related  to  the  plasma  mass,  mp# 
through  the  geometry  of  the  plasma  and  the  material  properties  (  and  p. 

We  do  not  know,  a  priori,  very  much  about  the  plasma  geometry;  however, 
it  appears  to  be  reasonable  to  assmne  that: 

Assumption  No.  6  -  The  plasma  body  will  scale 

geometrically  with  pulse  mass  (i.e.  ,  ratios  of 
dimensions  remain  fixed  as  the  thruster  size  and 
pulse  mass  are  varied). 

Assmptions  2  and  4  imply  that  both  plasma  resistivity  and  plasma  density 
are  constant  and  independent  of  pulse  mass  and  exhaust  velocity.  From 
Equation  B-15.  we  therefore  conclude  that  plasma  (and  thruster)  dimensions, 
d,  scale  as: 


d  a  mjl/3  ( B-16) 

and  from  Equation  B-14  plasma  resistance  scales  as: 

R  a  m“1/3  ( B-17) 


which  may  be  rewritten  as: 


R  b  km 


-1/3 

P 


( B-18) 


where  k  is  a  constant  to  be  determined  from  experimental  measurements. 


Plasma  dimensions  and  resistance  scale  somewhat  differently  for  the  MPD 
which  operates  in  a  quasi— steady  rather  than  a  poised  mode.  The  MPD  channel 
is  completely  filled  with  plasma.  The  density  and  current  have  a  stable 
distribution  during  the  entire  poise.  The  current  flows  radially  in  the 
plasma.  From  asstmptions  1,  3  and  4  we  can  write  the  plasma  resistance  in 

an  MPD  device  in  terms  of  an  effective  (constant)  resistivity  as: 


R 


$ 


dr/ (2nr l) 
ri 


where  t  m  the  length  of  the  MPD  channel 
*  inner  electrode  radins 
r0  m  outer  electrode  radius 

The  integration  may  be  completed  to  yield: 


(B-19) 


R  -  51n(r0/ri)/2n.£  (B-20) 

From  assumption  no.  6  we  see  that  rQ/  remains  constant  as  the 
thruster  is  scaled  and  we  conclude  that: 

R  a  1/  l  (B-21) 

From  the  data  available  on  MPD  thrusters  we  find  that  the  mass  flow  rate 
scales  as  the  annular  area  of  the  device: 

■p  «  ^d'tdg/r^2)  ( B-22) 

From  assumption  6  we  can  therefore  infer  that: 

l  a  r0  0  rA  a  ij/2  (B-23) 

which  may  be  combined  with  Equation  (B-21)  to  yield: 

R  -  km“1/2  (B-24) 


where  k  is  a  constant  to  be  detemined  from  MPD  performance  data. 


The  scaling  relation  for  the  inductance  gradient,  L'  ,  is  the  remaining 
on  known.  He  inductance  gradient  is  a  function  of  the  geometry  of  the 
thruster  electrodes.  Parallel  plate  devices  (such  as  the  PIT)  have  a 
gradient  given  approximately  by: 

L'  *=  p0w/h  (B-25) 

where  |iQ  =  permeability  of  free  space 
w  =  separation  of  plates 
h  ■  height  of  plates 

The  gradient  for  coaxial  devices  (such  as  the  MPD)  is: 

L'  =  (n„/2n)ln(r  /r^)  (B-26) 

o  o  * 

where  r0  *  radius  of  outer  electrode 
r^  =  radius  of  inner  electrode 

From  assumption  no.  6.  L'  is  independent  of  pulse  mass  (thruster  size)  and 
exhaust  velocity. 

Combining  the  inductance  and  resistance  results  for  the  PIT  and  TPP  and 
collecting  all  of  the  constants  together  we  obtain,  (Equations  B-12,  B-18) : 

Edis  -  (B-27) 

where  k  is  a  constant  which  will  be  determined  from  performance  data.  For 
the  MPD  device  we  obtain  a  slightly  different  result,  (Equations  B-13 ,  B- 
24)  : 

Edi.  -  V^SyVVe^S  <®-28) 

where  kj  and  k2  »re  constants  to  be  determined  from  performance  data. 

The  pulse  energy  equation  can  now  be  rewritten  for  the  PIT  and  TPP  as: 

E  *  ■  v2/2  +  m  e .  +  km2^8v»  (B~29) 

p  p  e  pi  P  e 

and  the  efficiency  as: 
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q  “  1/  1 1  2e^/v^  +  2k/ (  ] 

For  the  MPD,  the  equivalent  expressions  are: 


(B-30) 


EP  "  VeV2  *  Vpai  *  V^Va  +  V  Ve)1/2tp  (B_31) 

and 

r\  -  1/  [  1  +  2ei/y2  +  2k1/(iJ/2ve)+2k2/(i2/2v^/2)]  (B-32) 

All  of  the  available  data  on  efficiency  for  varying  pulse  mass  and 
exhaust  velocity  was  collected  for  each  device.  The  values  for  e^  and  k  (k^ 
and  k2  for  the  MPD)  which  best  fit  the  data  were  then  determined  and  are 
shown  in  Table  B-2 .  Also  shown  in  Table  B-2  is  the  specific  energy  re¬ 
quired  to  reach  the  first  ionization  level  for  argon  and  Teflon  (argon  is 
the  propellant  used  in  the  PIT  and  MPD  and  Teflon  in  the  TPP)  .  The  param¬ 
eter  e^  is  just  a  parameter,  however,  if  the  physical  model  is  good,  then 
its  value  should  be  comparable  to  that  of  the  actual  ionization  energy.  The 
results  shown  in  Table  B-2  support  the  model.  The  differences  between  e^ 
and  the  ionization  potential  indicate  the  model  is  not  exact. 


TABLE  B-2.  PULSE  ENEBGY  COEFFICIENTS 


ille 


The  efficiency  data  and  model  predictions  are  illustrated  in  Figures  B- 
1.  B-2,  and  B-3 .  Inspection  of  these  figures  indicates  a  good  fit  to  the 

available  data.  Closer  examination  of  Equations  B-30  and  B-32  shows  excel¬ 
lent  qualitative  agreement  with  observed  trends: 


Figure  B-l .  Efficiency  Data  and  Model  Results  for  the  PIT. 2 


Figure  B-2  .  Efficiency  Data  and  Model  Results  for  the  TPP.3»4 
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Figure  B-3 .  Efficiency  Data  and  Model  Results  for  the  MPD. 

(1)  Efficiency  increases  with  exhaust  velocity  (but  can  never  exceed 
100*1 )  . 

(2)  The  increase  of  efficiency  with  exhaust  velocity  is  most  pro¬ 
nounced  at  low  velocities  (<10,000  m/s). 

(3)  Efficiency  increases  with  pulse  mass  (thruster  size). 

(4)  The  effect  of  pulse  mass  is  pronounced  below  a  critical  size  where 
efficiency  drops  rapidly  with  decreasing  pulse  mass. 

(3)  At  pulse  masses  above  the  critical  value,  efficiency  is  insensi¬ 
tive  to  pulse  mass. 


B.3  CURRENT/ PULSE  TIME 

An  average  current  can  be  defined  from: 


(B-33) 


’Personal  Communication 
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which  can  be  farther  reduced  by,  substitution  of  Equation  B— 11,  to: 

J«vg  -  (2Ve/L‘tp)1/2  (B-34) 

Equation  B-34  embodies  most  of  the  desired  model  characteristics.  The 
inductance  gradient,  L* ,  is  independent  of  mp  and  yq.  The  average  current 
therefore  depends  directly  on  the  impulse  bit.  However,  the  average  current 
and  the  pulse  time  are  obviously  interdependent.  Assessment  of  current 
requirements  must  be  preceded  by  an  analysis  of  pulse  time.  Each  of  the 
thrusters  (PIT,  TPP  and  MPD)  have  different  pulse  time  constraints  and  must 
be  treated  separately. 

B.3.1  The  PIT 

The  PIT  accelerates  a  flat  sheet  of  plasma  by  repulsion  from  a  flat 
coil.  As  the  separation  between  the  plasma  sheet  and  the  coil  increases, 
the  mutual  inductance  (the  PIT  equivalent  of  L')  decreases.  The  accelerat¬ 
ing  force  then  drops  and  the  plasma  becomes  'decoupled*  (i.e.,  current 
flowing  in  the  drive  coil  no  longer  accelerates  the  plasma).  The  data  in 
the  literature  indicates  relatively  good  (and  constant)  coupling  up  to  a 
separation  of  about  1/10  of  the  coil  diameter  (the  decoupling  distance 
should  scale  directly  with  the  coil  diameter).  Coil  diameter  will  scale 
with  pulse  mass  according  to  assumption  no.  6  as  expressed  in  Equation  B-16. 
Correlating  the  available  data^  we  find  that  coil  diameter,  D,  is: 

D  =  66.5*1/3  (B-3  5) 

Therefore,  the  acceleration  length,  t,  is  given  by: 

l  =  6.5ml/3  (B-36) 

The  measured  inductance  gradient  for  the  PIT  is: 

L’  =  7nH/m  (  B— 3 7 ) 


^Personal  Communication  -  Dailey,  C. L. ,  THY  Space  and  Technology  Group. 
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(4)  The  models  fit  the  available  experimental  data  extremei'  wt  1  1 
(considering  the  simplicity  of  the  models  and  the  paucity  of 
data)  . 

(5)  The  excel  lent  qualitative  and  quantitative  agreement  between  mode  1 
prediction  and  observed  performance  lend  considerable  credence  to 
the  models  and  the  underlying  assumptions. 

(6)  The  models  are  extremely  versatile,  permitting  rapid  evaluation  of 
thruster  input  requirements  over  a  wide  range  of  output  condi¬ 
tions.  The  models  can  easily  be  updated  if  and  when  new  data 
becomes  available. 

B .  6  RECOMMENDATIONS 

There  are  two  major  areas  of  uncertainty  in  the  thruster  models  de¬ 
veloped.  The  first  arises  from  the  nature  of  the  thrusters  for  which  data 
was  available.  It  is  not  clear  that  the  experimental  thruster  performance 
data  available  to  us  was  generated  with  'optimum'  thrusters  (i.e.,  that  the 
thruster  was  optimized  for  every  mp  and  ve  for  which  data  was  reported). 
The  performance  data  must  be  evaluated  more  carefully  in  order  to  answer 
this  question.  The  second  problem  arises  from  the  limited  amount  of  data 
available.  The  missions  considered  in  this  study  require  that  the  thruster 
models  apply  at  operating  conditions  (m^  and  vg)  which  could  be  orders  of 
magnitude  different  from  the  conditions  for  which  we  have  data.  We  recom¬ 
mend  that  performance  data  for  optimum  thrusters  be  generated  over  an 
expanded  range  of  mp  mp)  and  ve- 
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The  acceleration  voltage  is  derived  from  Equation  B-56  and  is: 

V,cc  -  2.07xl0-4i£/2ve3/2  (B-70) 

Hie  se  voltage  ccaponents  are  shown  in  Figure  B-16. 
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Figure  B-16.  Voltage  Components  for  the  MPD. 


B.5  CCNCLO  SIGNS 

A  number  of  specific  conclusions  can  be  drawn  from  this  investigation 
of  performance  modeling  for  electromagnetic  thrusters. 

(1)  Thruster  modeling  for  integrated  propulsion  system  optimization 
requires  that  thruster  input  requirements  (pulse  energy,  current, 
voltage,  and  time)  be  related  to  the  primary  propulsion  system 
outputs  (pulse  mass,  exhaust  velocity,  frequency,  and  total 
thrusting  time) . 

(2)  Models  which  meet  these  requirements  were  successfully  developed 
for  all  of  the  thrusters  of  interest. 

(3)  The  models  developed  for  the  thrusters  (except  for  the  EBG)  were 
based  on  an  understanding  of  the  physical  principles  of  operation 
and  a  ntmber  of  assmptions  concerning  the  performance  of  opti¬ 
mized  thrusters. 
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of  the  total  poise  tine.  The  ionization  voltage  can  then  be  written  frets 
Eqoa  tion  B-53  as  : 

Vion  =  7.5xl08»J/3/ve  ( B-64) 

The  acceleration  voltage  is  derived  from  Equation  B-55  as: 

VaCC  =  0.311aJ/3ve  (B-65) 

The  dissipation  voltage,  computed  from  Equation  B-57  is  a  constant: 

Vdis  *  87  .2  ( B-66) 

These  voltage  components  are  shown  in  Figure  B-14.  The  TPP  thruster  does 
not  require  a  'breakdown  voltage',  since  the  breakdown  across  the  electrodes 
is  initiated  by  an  external  ignition  circuit.  There  is  a  minimum  allowable 
voltage  for  operation  of  the  device.  Lack  of  experimental  data  precluded 
developing  a  relationship  for  this  parameter. 
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iigure  B-14.  Voltage  Components  for  the  TiJP. 
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Figure  B-ll.  Voltage  Components  for  the  PIT. 

'ionization  voltage'  but  reflects  the  practical  requirement  of  achieving 
adequate  electric  field  to  ionize  the  gat.  Fran  the  experimental  data. 
Equation  B-53  was  reevaluated  to  yield: 

Vbdn  -  7.5xl05«^/3  (B-62) 

which  is  plotted  in  Figure  B- 12.  Comparison  of  Equation  B-62  with  Equation 
B-59  shows  that  the  breakdown  voltage  is  higher  than  our  'ionization  volt¬ 
age'.  This  indicates  that  ionization  will  be  completed  in  less  than  10%  of 
the  pulse  time. 

B.4.2  The  TPP 

The  average  voltage  for  the  TPP  thruster  can  be  found  from  Equation 
B-52  and  is : 


V^g  -  0.311a>p/3%  +  7.5xl07»£/3/ve  +  87.2  (B-63) 

which  is  shown  plotted  in  Figure  B-13 . 
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and  is  plotted  is  Figure  B-10. 


Figaro  B-10.  Average  Voltage  for  the  PIT. 


Bach  of  the  voltage  components  was  examined  separately.  Assisting  that 
the  plasma  is  formed  daring  the  first  10%  of  the  total  poise  duration,  the 
factor  (tp/t)  is  eqnal  to  10.  The  ionization  voltage  can  be  written  from 
Eqna  ti  on  B-S3  as : 

Vion  -  4.6lxl05»J/3  ( B-59) 

The  acceleration  voltage  derived  from  Equation  B-55  is: 

Vacc  -  2.59zlO"4aJ/3v^  (B-60) 

and  the  dissipation  voltage  derived  frm  Equation  B-57  is: 

Vdi *  -  7.52zl0_2ve  (B-61) 

The  three  voltage  ccmponents  are  plotted  in  Figure  B-ll  . 

The  PIT  also  has  a  minimtm  allowable  voltage  required  for  gas  breakdown 
and  plasma  formation.  This  'breakdown  voltage'  is  identical  to  our 


The  ionization  voltage,  as  defined  here,  should  not  be  confused  with  the 
voltage  required  to  ionize  (the  breakdown  voltage).  The  factor,  (t/tp>.  is 
the  fraction  of  the  pulse  duration  in  which  the  ionization  takes  place.  For 
example,  the  PIT  and  TPP  require  a  very  high  voltage  for  plasma  formation  at 
the  beginning  of  the  pulse.  Once  the  plasma  has  formed,  this  voltage  is  no 
longer  required.  On  the  other  hand,  the  MPD  thruster  requires  voltage  for 
ionization  of  propellant  throughout  the  pulse.  For  the  MPD,  (tp/t)  «  1  and 
the  'ionization'  voltage  is  written  as: 

Vion  *  *pei/J  (B“54> 

The  voltage  required  to  deliver  the  acceleration  power  to  the  pulse 
mass  is  the  acceleration  (or  useful)  voltage.  It  is  related  to  the  kinetic 
energy,  the  average  current,  and  pulse  duration  as: 

Vacc=V«/(JV  (B~S5) 

For  the  MPD  the  acceleration  is  quasi-steady  and  the  acceleration  voltage 
can  be  written  as: 


Vacc  “  Ve/J  (B~56) 

The  voltage  that  is  associated  wi  th  the  energy  dissipated  during  the 

pulse  is  related  to  the  dissipated  energy,  the  average  current,  and  pulse 
duration.  The  dissipative  resistance  of  the  thrusters  was  assumed  to  be  a 
constant  throughout  the  pulse.  The  dissipated  voltage  is  given  by: 

vdi,  •  (B-57> 

Ve  now  examine  these  voltages  for  each  thruster. 

B.4.1  The  PIT 

Since  all  of  the  factors  in  Equation  B-52  are  known,  the  PIT  average 
voltage  can  be  written  as: 

V^g  -  2.59xl0-4«ip/3v^  +  4.61xl04m2/3  +  7.52xl0“2ve  (B-58) 
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Figure  B-9.  Current  for  the  MPD. 

The  average  voltage  is  computed  from  the  pulse  energy,  the  pulse  time 
and  the  average  current  as: 

Vtvg  “  V(tPJ»vg}  (B~52) 

The  origin  of  the  voltage  can  be  better  understood  if  it  is  written  as  the 

sum  of  voltage  components  related  to  the  components  of  pulse  energy.  These 

voltage  exponents  are  denoted  by:  Vion  ~  v°lta8e  required  to  heat  and 

ionize  the  propellant,  V&cc  -  voltage  required  to  accelerate  the  propellant, 
and  “  the  voltage  characteristic  of  the  energy  dissipated.  A  minimum 

voltage  required  to  initiate  operation  (br  eakdown  vol  tage )  has  also  been 

analyzed  for  some  devices. 

The  'ionization'  voltage  is  computed  frem  the  ionization  energy 
(Equation  (B-6))  substituted  into  Equation  (B-52)  to  yield: 

Vion  -  Vi/(JV  X  (VU  (B_53) 
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B.3.3  The  MPD 


Since  the  MPD  operates  in  a  '  qua  si- steady’  mode,  the  pulse  duration  is 
an  independent  parameter.  The  current  is  therefore  independent  of  the  pulse 
time  and  is  assned  constant  during  the  pulse.  The  electromagnetic  thrust 
for  an  MPD  thruster  can  be  written  as: 

F  -  ipVe  =  L'J2/2  (D-48) 

This  can  be  generalized  to  form  a  relationship  between  the  current  and 
the  independent  pulse  parameters: 

J  a  (»pve)1/2  (B-49) 


or  written  in  equation  form  as: 


J  =  k(mpve)1/2  ( B-50) 

where  k  is  an  experimentally  determined  constant.  From  the  available  data5 
the  constant  was  evaluated  and  Equation  B-SO  becomes: 

J  •=  2.42xl03(i  vj1/2  ( B-51) 

p  e 

which  is  illustrated  in  Figure  B-9. 

B.4  VOLTAGE 

The  thrusters  treated  are  all  electromagnetic  and  are  'current  driven' 
(i.e..  the  acceleration  force  is  related  to  current,  not  voltage).  These 
devices  do  not,  in  general,  require  a  high  voltage  to  accelerate  the  pulse 
mass  (although  some  require  a  relatively  high  voltage  for  the  formation  of 
the  current-carrying  plasma  and  to  produce  the  required  short  pulses).  The 
voltage  is,  therefore,  a  derived  parameter  in  electromagnetic  thrusters  and 
is  indicative  of  the  impedance  of  the  device.  The  primary  function  of  the 
ES/PC  system  is  to  provide  the  required  current. 
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The  poise  tine  can  be  computed  by  rewriting  Equation  B-36  as: 


‘p  -  2»pve/ (B-45) 

Unfortunately  for  the  TPP,  we  don't  know  either  L'  or  Jjvg  (°r  ewen  the 
relationship  of  J^g  to  JBJLX)  .  If  the  acceleration  poise  is  similar  in 
scaled  devices  then  the  ratio  of  J^g  to  JM1  will  be  constant  and  Equations 
B-44  and  B-4S  can  be  combined  to  yield: 


t 


P 


(B-46) 


where  k  is  a  constant  to  be  determined  from  performance  data.  Available 
datal.0*ll  was  ttsed  to  obtain: 


tp  =  150xl0_9m^/3ve 


( B-47 ) 


which  is  illustrated  in  Figure  B-6 . 


Exhaust  velocity  ik"/-, ; 

Fignre  B-8.  Pulse  Time  for  the  TPP. 


*®Vondra,  R.  J.  ,  Thcmassen,  K.  ,  Solbes,  A.,  'Analysis  of  Solid  Teflon  Pulsed 
Plasma  Thrusters,'  Journal  of  Spacecraft  and  Sockets,  Vol.  7,  No.  12. 

^Poeschel,  R.  L.  ,  'A  Comparison  of  Electric  Propulsion  Technologies,'  AIAA- 
82-1243,  Presented  at  the  18th  Joint  Propulsion  Conference  (Cleveland,  OB). 


Assmaption  No.  7  -  For  TPP  thrusters,  the  magnetic 
pressure  is  maximized  aud  is  independent  of 
thruster  size  and  exhaust  velocity. 


The  maximum  current  scales  directly  with  thruster  size  according  to  Equation 
B-42.  From  assumption  no.  6  and  Equation  B-16,  thruster  size  scales  with 
pulse  mass.  Combining  Equations  B-16  and  B-42  we  obtain: 

J«ax  «  kaj/3  (B-43) 

where  k  «  a  constant  to  be  detemined  from  data.  The  available  data^  were 
used  to  evaluate  the  constant  and  the  result  is: 

Jmax  *  1.07x107b^/3  (ft-44) 

which  is  illustrated  in  Figure  B-7 . 


Fieare  B-7.  Maximum  Current  for  the  TPP. 


^Personal  Coaununi cation  -  Dominic  Palumbo  -  Fairchild  Industries,  Inc. 
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Figaro  B-6 .  Peak  Current  for  the  PIT. 


B.3.2  The  TPP 

The  TPP  is  a  parallel  plate  electrode  device  and  the  electrodes  can, 
presumably,  be  made  as  long  as  necessary  to  reach  the  desired  velocity.  The 
poise  time  is  therefore,  not  constrained,  and  ve  must  search  ror  other 
constraints  on  the  current. 

The  TPP  operates  in  a  very  similar  manner  to  the  ERG.  The  primary 
difference  between  the  two  is  the  presence  of  a  solid  'pellet'  in  the  ERG. 
Optimum  TPP  design  should  be  similar  to  optimum  ERG  design.  The  optimi¬ 
zation  strategies  for  ERG  design  are  smaarized  in  Table  B-l .  One  of  the 


primary  optimization  strategies  is  to  maximize  the  magnetic  (accelerating) 
stress.  The  stress  is  related  to  current  and  electrode  dimensions  by: 


a  a  (J/h)2 


( B-42) 


vhere  h  *  the  electrode  height.  'Die  mailman  stress  will  be  constrained  by 
thruster  strength,  electrode  erosion  (from  the  resultant  current  density)  or 
some  other  factor.  For  scaled  TPP  thrusters  designed  for  optimum  efficiency 
ve  would  ezpect  that: 
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Figure  B-5 .  Average  Current  for  the  PIT. 


measured  pulse  times  are,  as  expected,  longer  than  Equation  B-39  predicts, 
but  only  by  about  15%.  There  are  a  few  measurements  of  peak  current 
reported  and,  as  expected,  they  correlate  with  pulse  mass  and  exhaust 
velocity  as  predicted  by  Equation  B-41 .  The  data^»8  for  peak  current 
yiel  ds  : 

Jmax  =  261®p/3ve  (B~41) 

which  is  illustrated  in  Figure  B-6 . 

Comparison  of  Equations  B-40  and  B-41  indicate  that,  for  the  particular 
ES/PC  system  used  in  the  experiments,  the  peak-to- aver  age  current  ratio  was 
approximately  1.75. 


^Dailey,  C. L. ,  'Large  Diameter  Inductive  Plasma  Thrusters,'  AIAA- 79-2 093 , 
Presented  at  the  14th  International  Electric  Propulsion  Conference  (Prince¬ 
ton,  NT )  . 

® Da i 1 ey ,  C.  L. ,  'Large  Diameter  Inductive  Plasma  Thruster,1  AIAA-79-2  093 , 
presented  at  the  14th  International  Electric  Propulsion  Conference  (Prince¬ 
ton,  MI )  . 
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SECTION  C.l 


INTSODCCTION 


Inductor*  can  be  used  either  as  energy  stores  or  as  pul se-f orming 
components  in  pulse  power  systems.  In  the  energy  storage  role  inductors  are 
attractive  because  they  offer  high  stored  energy  at  relatively  loir  mass.  In 
addition,  they  provide  a  1  ov- impedance,  constant-current  source  characteris¬ 
tic.  This  characteristic  »s  very  important  for  seme  applications  such  as 
electromagnetic  launchers  and  electromagnetic  thrusters.  The  disadvantages 
of  inductors  as  energy  stores  is  that  they  are  lossy  (i.e.,  have  a  low  Q) 
and  are  therefore  relatively  inefficient.  The  losses  can  be  eliminated  and 
the  efficiency  improved  by  employing  superconductors.  Superconductors, 
however,  are  nut  tolerant  to  pulse  discharge  and  pose  serious  thermal  design 
problems  in  high  current  applications. 

Inductors  are  also  used  as  pulse-forming  components  to  control  the 
shape  and  duration  of  pulses  delivered  from  capacitive  energy  stores.  In¬ 
ductors  can  be  used  in  either  a  normal  LC  combination  or  in  multi-staged 
pulse  forming  networks.  In  these  applications,  the  inductors  carry  current 
for  only  a  short  time  and  the  resistive  losses  inherent  in  inductors  are  not 
so  important. 

In  many  applications  of  interest  for  both  inductive  energy  storage  and 
inductive  pulse  forming,  component  mass  must  be  minimized.  The  objective  of 
this  report  is  to  develop  methods  of  minimmi  mass  selection  and  sizing  of 
inductors  for  both  energy  storage  and  pulse  forming  duty. 


SECTION  C.2 


CIRCUIT  CHARACTERISTICS 

The  circuit  performance  of  inductors  plays  an  important  role  in  induc¬ 
tor  sizing  and  selection.  The  energy  efficiency  of  the  inductor  in  the 
circuit  is  of  particular  interest.  Te  define  the  efficiency  as: 

n  “  Euae/(Edis+5W  (C“1> 

There  EdiJ  *  dissipated  energy. 

*  ueful  energy  extracted, 
lhe  dissipated  energy  is: 

Edis  -  R/I2dt  (C-2) 

where  S  «  resistance  of  the  inductor. 

The  useful  energy  is: 


Euse  -  UI2  -  if )/2  (C-3) 

where  I0  ■  maximum  current  during  period  of  interest 
If  3  minimum  current  during  period  of  interest 
Substituting  Equations  2  and  3  into  Equation  1  provides: 

H  -  1/(1  ♦  20/t)  ( C-4) 

where  0  »  /l2dt/(I2  -  I2),  the  efficiency  factor 
r  3  L/R.  the  inductor  time  constant. 

The  factor.  0,  depends  on  the  details  of  the  application  and  the  character¬ 
istics  of  the  circuit. 

In  pul  se- forming  duty,  0  can  be  readily  evaluated  from  the  circuit 
performance  characteristics  and  is  approximately  equal  to  the  pulse  dura¬ 
tion.  If  the  inductor  time  constant  is  long  with  respect  to  the  pulse 
duration  the  efficiency  will  be  high. 
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figure  C-l.  Current  Variation  in  Energy  Storage. 

In  energy  storage  duty,  0  is  related  to  the  frequency  of  poising  and 
the  depth  of  discharge.  He  basic  current-time  characteristic  for  energy 
storage  duty  is  illustrated  in  Figure  C-l.  The  discharge  time  is  typically 
orders  of  magnitude  shorter  than  the  charging  time.  He  efficiency  of  the 
inductor  is  therefore  primarily  determined  by  the  charging  process.  Two 
different  charging  variations  are  illustrated,  a  linear  current  rise  and  a 
square  root  current  rise  (which  more  nearly  approximates  the  desired  con¬ 
stant  power  charging).  He  factor  0  can  be  evaluated  for  both  cases  and  is 
related  to  the  frequency  by: 


0  -  a/2f  (C-5) 

where  a  =  an  efficiency  factor.  He  efficiency  factor  is  a  function  of  the 
current  ratio,  If/Io,  and  is  illustrated  in  Figure  C-2.  For  reasonable 
current  ratios  (<0.5)  the  efficiency  factor  is  independent  of  current  ratio 
and  is  approximately  equal  to  one. 

He  efficiency  can  be  evaluated  from  Equation  C-4.  He  results  for 
inductors  with  various  time  constants  are  illustrated  in  Figure  C-3.  He 
inductor  time  constant  and  its  dependence  on  inductor  mass  is  of  major 
importance.  I 
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Figaro  C-3.  Inductor  Efficiency. 


SECTION  C.3 


INDUCTOR  SIZING 


C.3 .1  THE  INDUCTANCE 

Fran  die  point  of  view  of  inductor  sizing  and  design,  die  inductance  is 
not  a  design  paraaeter.  The  designer  is  not  free  to  choose  the  inductance 
in  order  to  minimize  mass.  Hie  inductance  is  usually  fixed  by  the  appli¬ 
cation  and  the  designer's  job  is  to  size  and  design  an  inductor  with  the 
required  inductance  and  ainiami  aass.  Froa  the  designer's  point  of  view, 
therefore,  the  inductance  is  an  independent  paraaeter. 

The  inductance  of  an  inductor  is  dependent  upon  the  configuration  and 
the  size  of  the  inductor.  The  inductance  asy  be  conveniently  expressed  as: 

L  *  kNz  (C-6) 


where  L 
k 

N 

z 


■  inductance 
“  inductance  factor 
*  umber  of  turns 
>  length  of  conductor. 


The  geometrical  factor,  k,  depends  on  the  choice  of  configuration  (toroid, 
solenoid,  etc.)  and  on  the  geometry  of  the  inductor  ( length-to-radius  ratio, 
etc.).  The  value  of  the  inductance  factor  is  independent  of  coil  size. 
Values  for  the  inductance  factor  can  readily  be  calculated  for  all  common 
configurations  and  geometries.  In  this  report  we  concentrate  on  soleniods 
as  they  are  simple  and  display  the  same  basic  scaling  characteristics  as 
other  multi-turn  concepts  such  as  toroids. 


The  geometry  of  a  solenoid  is  illustrated  in  Figure  C~4.  The  induc¬ 
tance  tactor  for  solenoids  can  readily  be  evaluated  from  the  formulas  of 
Grover^  and  Wheeler.^  The  inductance  factor  is  a  tunction  of  the  geometri¬ 
cal  ratios  c/a  and  b/a  and  is  illustrated  in  Figure  C-5.  Ve  note  that  for 


^Grover,  F.  W.,  'Inductance  Calculations,'  D.  Van  Nostrand  Co.,  New  York, 
1946. 

^Wheeler,  H.  A.,  'Inductance  Formulas  for  Circular  and  Square  Coils,'  IEEE 
Proc.,  70,  pp  144  9-14  5  0,  19  82. 
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long  coils  (b/a  >  1)  tho  inductance  factor  is  not  a  strong  function  of  the 
thickness  ratio  (c/a).  The  parameters  for  a  Brooks  coil  are  also  identified 
in  Figure  C-S.  The  Brooks  coil  is  important  as  it  is  the  solenoid  config¬ 
uration  which  has  the  highest  inductance  for  a  given  length  of  conductor. 
It  is  sometimes  thought  of  as  a  minimum  mass  configuration.  However,  it  is 
only  a  minimum  length  configuration.  As  we  shall  see  in  the  following 
sections,  inductor  mass  depends  on  other  factors  as  well  as  conductor 
length. 

The  inductance  factor  for  solenoids  can  be  erpressed  approximately  as: 

k  3  k  -  ki  (C-7) 

o 

where  kQ  =  2xl0~7 [2 . 78/ ( (b/ a ) +1)  +  ln ( 1+0  .3 9(b/a ) ) ] 
kj^  *  6  .2  8xl0~^(»/b)  (c/a)k' 

k'  =  (0. 331-0. 011/ (b/a) )oxp[-<0. 361+0. 088/ (b/a) ) (c/a )  ]  . 

The  approximation  is  very  good  over  the  range  of  geometries  of  interest. 
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Figaro  C-5 .  Iaductan.ce  Factor  for  Solenoids. 


C.3 .2  IN DOCTOR  MASS 

The  mass  of  a  solenoid  can  be  expressed  as: 

■  *  pAz  ( C-8) 


where  ■  =*  inductor  mass 

p  =■  mass  density  of  the  conductor 
A  •  cross-sectional  area. 

Inductor  mass  is  minimized  by  minimizing  the  product  of  z  and  A.  Inspection 
of  Equation  C-6  shows  that  z  will  be  minimized  if  N  is  maximized  (for  a 
fixed  inductance)  and  that  inductor  mass  will  be  directly  related  to  the 
inductance.  Equations  C-6  and  C-8  provide  no  information  on  the  selection 
of  conductor  cross  section.  A.  For  a  minimum  mass  inductor,  A  must  be  as 
smal  1  as  possibl  e. 


C.3. 3  CONDUCTOR  CROSS-SECTIONAL  AREA 

The  conductor  cross-sectional  area  can  in  principle  be  reduced  until 
the  stress  and/or  the  temperature  (due  to  resistive  dissipation)  becomes 
unacceptable.  Further  reduction  in  cross  section  will  result  in  coil 
failure. 


C.3. 3.1  Thermal  Constraints 

The  thermal  design  of  an  inductor  can  proceed  in  either  of  two  direc¬ 
tions.  In  the  first  approach  the  inductor  is  not  actively  cooled,  but  is 
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Figure  C-ll.  Specific  Energy  vs.  Current. 


C.4J  INDUCTANCE 

The  effect  of  inductance  on  specific  energy  is  more  clearly  illustrated 
in  Figure  C-12.  The  specific  energy  increases  with  Figure  C-ll  inductance 
at  low  current  levels  and  is  independent  of  inductance  at  high  current. 


Figure  C-12.  Specific  Energy  vs.  Inductance. 
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small  c/a  (<0.2).  The  effect  of  thickne  ss-to-radins  ratio  is  il  lastrated 
more  clearly  in  Figure  C-9. 

C.4.3  CONDO CTOR  TEMPERATURE 

The  conductor  temperature  primarily  affects  the  resistivity  bat  also 
has  an  effect  on  heat  transfer.  The  variation  of  specific  energy  wi  th 
temperature  is  illustrated  in  Figure  C-10.  The  specific  energy  of  the 
inductor  is  almost  independent  of  temperature.  The  total  mass  of  the  system 
must  include  the  cryogenic  system  and  will  not  necessarily  show  the  same 
effects  as  the  inductor  only. 


Figure  C-10.  Specific  Energy  vs.  Conductor  Temperature  for  Al-1100-0. 

C.4.4  CURRENT 

The  effect  of  current  on  inductor  specific  energy  is  illustrated  in 
Figure  C-ll.  At  low  current,  the  specific  energy  is  proportional  to  the 
current  and  has  a  slight  dependence  on  the  inductance.  At  high  current,  the 
energy  density  is  independent  of  current  and  inductance. 
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C.4.1.3  The  Geometry 


Tie  geometry  was  varied  to  determine  inductor  mass  sensitivity  to  these 
factors.  The  ranges  examined  were 


0.03  <  c/a 

<  1 

(C-32) 

0.125  <  b/a 

<  16 

(C-33) 

s/A  =*  2500 

m~l 

(C-34) 

C.4.1.4  The  Ontpnts 

The  program  calculates  a  number  of  inductor  parameters.  For  the  pur¬ 
poses  of  comparison  the  specific  energy  (stored  energy  per  unit  mass)  was 
chosen  as  a  figure  of  merit.  The  higher  the  specific  energy  the  'better' 
the  inductor  (i.e.,  a  lower  mass  for  a  given  stored  energy).  This  parameter 
also  permits  rapid  comparison  of  inductors  with  other  types  of  energy  stor¬ 
age  for  which  specific  energy  is  usually  known.  The  specific  energy  was 
examined  over  the  entire  range  of  parameters  described  in  the  preceding 
paragraphs. 

C.4.2  IN  DOCTOR  GEOMETRY 

The  1  ength—  to-radius  and  the  thiclmess-to-radius  ratio  both  have  an 
effect  on  the  specific  energy  of  inductors.  A  typical  electromagnetic 
thruster  case  was  chosen  to  examine  the  effect.  The  inductance  was  lOnfl, 
the  current  was  lOOkA.  the  conduction  time  was  set  at  10^s  (an  energy 
storage  case)  and  a  temperature  of  300°S  was  assumed.  The  resnl  ts  are  shown 
in  Figures  C-8  and  C-9. 

Figure  C-8  shows  the  effect  of  length  to  radius  ratio  at  various  thick¬ 
ness  to  radius  ratios.  For  these  conditions  the  specific  energy  is  a  few 
hundred  joules  per  kilogram,  an  order  of  magnitude  higher  than  dielectric 
capacitors.  The  specific  energy  for  very  long  solenoids  is  approximately 
1/2  that  of  very  short  ones.  Increasing  the  th  i  c  kne  s  s- t  o- r  adi  us  ratio 
improves  the  specific  energy,  however  the  results  for  very  short  (b/a  <  1) 
solenoids  are  not  accurate.  The  formula  currently  used  to  compute  the 
inductance  factor,  k,  is  only  accurate  for  either  large  b/a  (>1)  and/or 
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SECTION  C.4 
RESULTS 


The  computer  program  was  exercised  over  a  wide  range 

eters.  The  purpose  was  to  examine  the  ability  of  the 

realistic  minimum  mass  inductors  and  to  determine  the 

such  inductors. 

of  inductor  pa ram- 

program  to  select 

characteristics  of 

C.4 .1  THE  PARAMETERS 

C.4  .1 .1  The  Independent  Parameters 

The  independent  parameters  are  the  inductance,  L,  the  current,  I,  and 

the  conduction  time,  t.  The  ranges  investigated  were  derived  from  electro¬ 
magnetic  thruster  considerations  and  were: 

O.lpH  <  L  <  25|iH 

(C-29) 

103A  <  I  <  106A 

(C-3  0) 

10*6s  <  t  <  lOyr 

(C-31) 

C.4. 1.2  The  Material  Properties 

The  properties  of  Al-1100— 0  aluminum  were  nsed.  The  performance  at 
three  temperatures,  21°X.  77°K  and  300°K  was  examined.  The  properties  used 
are  summarized  in  TABLE  1. 

TABLE  1.  MATERIAL  PROPERTIES,  AL-1100-0 


Temperature 

(°K) 

Densi  ty 

(kg/m2) 

Stress 

(MN/m3) 

Re  sistivi ty 

(0-m) 

Temperature 

Rise 

(°K> 

Action 

Constant 

(A2-s/m4) 

Hea  t 

Tr ansf  er 

( w/m3  ) 

21.00 

2700.00 

70.0 

0.110E-8 

50.0 

0.130E17 

0.100E5 

77.00 

2700.00 

70.0 

0.400E-8 

0.110E17 

0.500E5 

300.00 

2700.00 

70.0 

0.255E-7 

50.0 

0.300E16 

0.100E7 
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C.3.5  COMPUTATIONAL  TECHNIQUE 


Tie  minima  mass  configuration  is  fonnd  by  simultaneous  solution  of  the 
equations  developed  in  the  previous  sections.  A  closed  fom  solution  is  not 
possible  and  a  nsaerical  approach  aust  be  used.  A  computer  program  vas 
developed  to  find  the  solutions.  The  inductance,  current,  and  conduction 
tiae  vera  treated  as  the  independent  variables.  The  aaterial  properties  of 
resistivity,  density,  action  constant,  and  maximum  stress  aust  be  chosen. 
The  thermal  transfer  rate  and  perimeter— to-area  ratio  must  also  be  selected. 
Finally,  the  geometrical  factors  c/a  and  b/a  aust  be  chosen.  The  program 
then  finds  the  minimum  aass  indnctor  and  determines  the  mass,  resistance, 
number  of  turns,  stress,  and  dimensions  of  the  inductor.  The  choice  between 
active  and  passive  cooling  is  made  based  on  minimizing  mass.  Reduced  resis¬ 
tance  configurations  are  found  by  incrementing  the  conductor  cross  sectional 
area  and  repeating  the  computation.  The  stress  and  thermal  requirements  are 
reduced  in  reduced-resistance  inductors. 
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inductance  during  rapid  discharge.  During  charging,  current  and  magnetic 
field  diffuse  into  the  conductors.  Yhen  the  inductor  is  subsequently  dis¬ 
charged,  the  magnetic  field  diffused  into  the  conductors  (and  the  energy 
associated  vith  it)  cannot  be  recovered.  Only  the  field  and  energy  stored 
in  the  'open'  area  of  the  inductor  and  in  the  skin  depth  of  the  conductors 
can  be  recovered.  The  result  is  an  apparent  drop  of  inductance  during 
discharge.  The  effect  is  most  pronounced  in  thick  coils,  such  as  the  Brooks 
configuration,  vhere  reductions  of  20%  or  more  can  be  expected.  In  addi¬ 
tion,  for  multi-layer  coils,  internal  circulating  currents  are  induced  which 
result  in  significantly  increased  resistive  losses. 

These  high  frequency  effects  can  be  controlled  by  minimizing  inductor 
thickness,  c,  and  by  maximizing  inductor  thickness-to-radius  ratio,  c/a.  If 
c/a  is  constrained  to  be  less  than  about  0.1,  the  high  frequency  effects  are 
modest. 

C.3.4  INDUCTOR  RESISTANCE 
CJ  .4  J.  Minimum  Mass  Inductors 

In  the  preceding  paragraphs  the  conductor  geometry  for  minimum  mass 
inductors  was  developed.  The  resistance  of  these  inductors  can  be  readily 
cal  cul  ated  from 


R»5z/A  (C-28) 

The  cross  sectional  area.  A,  is  detexmined  from  thermal  constraints  and  the 
conductor  length,  z,  is  determined  from  stress  considerations.  The  minimum 
mass  constraints  result  in  maxima  resistance.  For  same  applications  this 
resistance  may  be  unacceptably  high.  The  resistance  can  always  be  reduced 
but  there  is  a  mass  penalty. 

C.l.4.2  Resistance-Mass  Scaling 

In  order  to  reduce  the  resistance  the  size  of  the  inductor  must  be 
increased.  The  process  begins  by  increasing  the  cross  sectional  area.  The 
size  and  other  characteristics  of  the  inductor  can  then  be  recalculated 
using  the  previously  developed  formulas.  There  is  no  simple  relationship 
between  inductor  mass  and  resistance. 
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AL  (1100-0) 


Figure  C-7 .  Tho  Electrical  Skin  Depth  for  Aluminum. 

C>3  J  .4 .1  Resistance 

For  poises  in  which  the  skin  depth  is  less  than  the  condnctor  thick¬ 
ness,  the  resistance  of  the  condnctor  increases  above  the  DC  resistance 
inverse:  7  with  the  skin  depth.  'Dierefore,  an  inductor  can  have  a  low  DC 
resistance  bat  yet  have  a  high  discharge  resistance. 

in  energy  storage  '  appl  ica  ti ons  the  inductor  must  respond  in  two  quite 
different  time  scales.  Charging  is  a  relatively  slow  process  related  to  the 
frequency  at  which  pulses  are  required.  The  discharge  time,  or  pulse  dura¬ 
tion,  is  usually  much  shorter  (typically  by  orders  of  magnitude)  than  the 
charging  time.  (Otherwise  an  energy  store  would  not  be  required.)  The  skin 
effect  resistance  is  not  usually  a  problem  during  charging,  but  can  be 
important  during  discharge.  The  thinner  the  conductors  are,  the  less  depen¬ 
dent  the  resistance  is  on  pulse  time.  Eie  high  frequency  resistance  of  an 
inductor  can  usually  be  designed  to  acceptably  low  values  without  seriously 
affecting  inductor  sizing  or  configuration. 

C.3.3.4.2  Inductance 

other  effect  of  diffusion  limited  skin  depth  is  a  reduction  of 
inductance.  This  effect  can  result  in  a  substantial  reduction  of  apparent 
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forces  the  coil  thickness,  c,  to  increase  (in  order  to  retain  the  required 
conductor  cross-sectional  area).  Under  these  circumstances,  the  required 
inductor  thickness,  c,  can  be  considerably  greater  than  the  mean  inductor 
radius,  a.  This  of  course  is  an  impossible  configuration.  The  stress  (and 
therefore,  the  number  of  turns)  must  be  reduced  in  order  to  achieve  a 
geometrically  possible  inductor.  The  implication  is  that  for  certain  appli¬ 
cation  regimes  a  highly  stressed  inductor  cannot  physically  be  constructed. 
The  geometrical  constraint  may  be  expressed  as: 

c/a  <  1  ( C-26 ) 


and  must  always  be  applied. 

C.3.3.4  High  Frequency  Effects 

Both  energy- st  or  age  and  pulse-forming  inductors  are  subjected  to  high 
frequencies  during  pulsing.  The  effects  of  high  frequencies  on  the  resis¬ 
tance  and  the  inductance  of  the  coil  can  be  important,  particularly  for  high 
current  designs  where  conductor  size  is  large.  When  current  is  changed 
rapidly  in  a  conductor,  current  diffusion  effects  confine  the  changing 
current  to  a  thin  surface  layer  or  skin.  The  depth  of  this  layer  is  related 
to  the  pulse  duration  by: 


cs  =  (n5tp/n0)1/2 


(C-27 ) 


where  cs  3  shin  depth 
tp  =*  pulse  time. 

The  skin  depth  for  aluminum  at  various  temperatures  is  illustrated  in  Figure 
C-7.  The  skin  depth  for  millisecond  pulses  (those  of  interest  here)  is 
less  than  1  cm.  If  the  skin  depth  is  less  than  the  conductor  thickness  then 
both  the  resistance  and  the  inductance  of  the  coil  are  affected. 
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ff  *  (i0(NI)‘ta/2?cb-‘ 


(C-22) 


We  aote  chat: 


NAg  -  cb  (C-23) 

and  Equation  C-23  can  be  reduced  to: 

•  “  <d0I2/2>(N/Ag)U/b)/p  (C-24) 

Equation  C-24  can  be  rearranged  to  provide: 

N  -  3<2oAg/tloi2)  (b/a)  (C-2S) 

The  total  number  of  turns  in  the  solenoid  is  therefore  determined  by  the 
allowable  stress  in  the  conductors. 

Inspection  of  Equations  06  and  C-8  shows  that  coil  mass  is  minimized 
by  maximizing  the  ntmber  of  turns  in  the  inductor.  From  Equation  C-25  we 

see  that  the  number  of  turns  is  maximized  by  maximizing  the  allowable  stress 

and  by  decreasing  the  allowable  current.  The  effect  of  inductor  aspect 
ratio  (b/a)  on  the  number  of  turns  is  not  so  clear  because  the  aspect  ratio 
|  also  has  a  direct  effect  on  the  inductance  factor,  k.  Increasing  the  geo¬ 

metrical  cross  section.  Ag,  increases  the  number  of  turns,  but  also  directly 
increases  inductor  mass. 

In  summary,  the  conductor  cross  section  is  sized  by  thermal  constraints 
while  the  number  of  turns  is  determined  by  allowable  stress.  Finally,  there 
are  two  other  constraints  which  also  must  be  considered  in  inductor  sizing. 
These  constraints  are  geometry  and  high  frequency  affects. 

C.3.3.3  Geometrical  Constraints 

Simultaneous  application  of  the  stress  and  thermal  constraints  de¬ 
scribed  in  the  previous  paragraphs  can  lead  to  impossible  physical  inductor 
dimensions.  The  problem  is  particularly  severe  for  low  current  and/or  long 
pulse  duration  inductors.  In  order  to  meet  the  stress  design  constraints 
(Eqn.  C-19)  the  number  of  turns  required  becomes  very  large.  This  in  turn 


The  presence  of  fins  and/or  cooling  channels  increases  the  geometrical 
cross  section  of  the  conductors.  The  geometrical  cross  section  and  the 
physical  cross  section  are  then  related  by: 

A  -  pAg  (C-18) 

where  Ag  =  geometrical  cross  section 
P  =  'solidity'  ratio  <_  1.0. 

The  solidity  ratio  coaid  be  as  low  as  0.5  in  order  to  achieve  maximum 
cooling  rates. 

C.3.3.2  Stress  Constraints 

A  solenoid  acts  as  a  magnetic  pressure  vessel.  The  magnetic  field 
exerts  a  bursting  pressure  on  the  solenoid  and  this  in  turn  induces  stress 
in  the  windings.  The  stress  in  the  conductors  of  a  solenoid  is  related  to 
the  magnetic  pressure  by: 


a  “  PBa/cP 


(C-19) 


where  a  =  conductor  stress 

Pm  =  magnotic  pressure. 

This  equation  is  strictly  true  only  if  the  stress  is  uniform  throughout  the 
conductors.  The  magnetic  pressure  is  related  to  the  ma^aetic  flux  density 
by: 


P*  =  B2/!^ 


(C-20) 


where  B  =  magnetic  flux  density 

L*0  =  permeability  of  free  space. 

The  magnetic  flux  density  is  related  to  inductor  geometry  and  current 
through  the  Biot-Savart  law  by: 

B  =  HoNI/b  (C-21) 

Combining  Equations  C-19,  C-20,  and  C-21  we  obtain: 
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by  increasing  the  beat  transfer  rate  (i.e.,  forced  convection/nuci  ear  boil¬ 
ing),  and  by  increasing  the  perime ter-to-area  ratio  (i.e.,  finning  or  chan¬ 
neling  of  the  conductors). 

Heat  transfer  rates  vary  fro*  a  few  kl/m?  for  cryogenic  cooling  to  a 
few  MV/ ar^  for  boiling  water  cooling.  The  maximum  achievable  transfer  rates 
depend  strongly  on  the  detailed  design  of  the  coolant  channels  and  the 
velocity,  pressure  and  Reynolds  Nnaber  of  the  coolant  flow.  Ve  shall  assume 
the  following  heat  transfer  values  for  this  study: 

4H  -  104w/a2  (C-13) 


for  liquid  hydrogen. 


4^  *  5xl04w/m? 


(C-14) 


for  liquid  nitrogen,  and 


qw  -  10** w/ a?  (C-15) 

for  water. 

The  periaeter-to-area  ratio,  s/A,  is  maximized  by  finning  and/or  by 
providing  many  saal  1  channels  for  coolant  flow  in  the  conductor.  In  addi¬ 
tion  the  surface  of  the  fins  or  channels  should  be  finely  grooved  to  further 
increase  the  available  cooling  area.  If  this  is  done,  the  ratio  is  related 
only  to  the  mean  thickness  of  the  conductor  fins  or  webs  by: 

s/A  =  4/d  (C-16) 

where  d  »  mean  thickness  of  the  fins  or  webs.  The  perimeter  to  area  ratio 
is  maximized  by  minimizing  the  fin  or  web  thickness.  The  minimum  thickness 
is  in  practice  probably  related  to  the  mean  conductor  thickness  (i.e.,  c  in 
Figure  C-4).  Ve  shall  assume  that  the  minimum  is  given  by: 


I 


d  *  0.1c 


(C-17) 


Figure  C-6.  Action  Constant  for  Aluminum. 


C.3.3.1.2  Active  Cooling 

Active  cooling  of  a  conductor  can  be  characterized  by  equating  the 
resistive  dissipated  power  to  the  thermal  transfer  from  the  cooled  surfaces. 
This  may  be  written  as: 


I2$z/A  -  q*z 


(C-ll) 


where  {  ■  resistivity 

q  =  heat  transfer  rate 

s  3  'wetted'  perimeter  available  for  cooling. 

Equation  C-ll  can  be  rewritten  as: 

A  =  I(5/4(s/A))1/2  (C-12) 

The  perimeter-to-area  ratio,  s/A,  is  a  function  of  the  geometry  of  the 
conductors.  The  cr  os  s- se  ct  i  ona  1  area  (and  therefore,  inductor  mass)  is 
reduced  by  reducing  the  resistivity  (i.e.. 


reducing  conductor  temperature). 


designed  as  a  thermal  inertia.  H e  conductor  temperature  is  allowed  to  rise 
and  the  cross  section  is  selected  to  keep  the  temperature  rise  acceptable 
for  the  required  current  level  and  the  required  conduction  time.  In  the 
second  approach,  the  conductors  are  actively  cooled  and  the  cross  section  is 
sized  to  permit  steady-state  transfer  of  the  dissipated  energy  into  the 
coolant.  The  approach  which  wil  1  provide  the  minimum-mass  coil  depends  very 
strongly  on  the  application.  The  tradeoffs  and  proper  selection  are  made 
during  computation. 

C.3JJ.J.  Thermal  Inertia 

The  required  cross-sectional  area  for  a  thermal  inertia  design  can  be 
readily  determined  from  the  action  integral: 

A2dt  -  gxA2  (C-9) 


where  I  3  current 

t  =»  time  of  conduction 

g^  *  action  constant  (dependent  on  temperature  rise) 

A  ■  cross-sectional  area. 

The  action  constant,  g^,  is  a  function  of  the  conductor  material,  the  ini¬ 
tial  temperature,  and  the  desired  temperature  rise.  The  values  for  aluminum 
are  illustrated  in  Figure  C-6. 

Assuming  a  constant  current.  Equation  C-9  can  be  rewritten  as: 

A-Kt/g^1'2  (C-10) 

The  cross-sectional  area,  and  therefore,  the  inductor  mass,  can  be  reduced 
by  reducing  the  current  or  the  pulse  time  or  by  increasing  the  action 
constant.  The  action  constant  can  be  increased  by  material  selection,  by 
increasing  the  allowable  temperature  rise,  or  by  decreasing  the  initial 
temperature.  Copper  has  a  higher  action  constant  than  aluminum,  however, 
that  is  more  than  compensated  by  the  higher  density  of  copper.  Aluminum  is 
usually  the  best  overall  material  choice  for  minimum  mass  inductors. 
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C.4.6  CONDO CT ION  TIME 

The  specific  energy  is  si  so  s  function  of  the  conduction  time  end  is 
shown  in  Figure  C-13.  In  pul  se-f  oraing  duty,  tine  corresponds  to  the  poise 
tine.  In  energy  storage  application,  conduction  is  usually  continuous  or 
nearly  so.  For  very  short  pulse  tiaes  the  inductor  is  sized  as  a  thermal 
inertia.  The  ainiaua  conductor  cross  section  is  then  directly  related  to 
the  pulse  tiae.  He  specific  energy  of  such  inductors  can  be  very  high.  As 
pulse  tiae  is  increased  the  inductor  auat  be  actively  cooled  and  the  conduc¬ 
tor  cross  section  is  independent  of  conduction  tiae.  The  specific  energy  is 
usually  lower  for  long— pulse  actively— cool  ed  inductors  than  for  short-pulse 
thermal  inertia  designs. 
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Figure  C-13.  Specific  Energy  vs.  Conduction  la.me. 


SECTION  C.5 


CONCLUSIONS 


Sev «r al  important  conclusions  can  be  drawn  from  tbis  study  of  induc¬ 
tors.  The  conclusions  are: 

(1)  Tie  specific  euergjr  of  normal  inductors  is  very  attractive  when  com¬ 
pared  to  that  of  dielectric  capacitors  and  is  comparable  to  inertial 
stores. 

(2)  Cryogenic  operation  has  little  effect  on  the  specific  energy  of  induc¬ 
tive  stores.  The  mass  of  the  cryogenic  system  must  also  be  included  to 
arrive  at  a  system  level  conclusion.  Cryogenic  operation  does,  how¬ 
ever,  significantly  reduce  the  inductor  resistance. 

(3)  The  mass  of  an  inductive  store  is  a  very  complex  function  of  material, 
operational,  and  application  parameters.  It  is  not  possible  to  arrive 
at  a  single  parameter  (such  as  specific  energy)  which  will  completely 
characterize  inductive  energy  stores  for  any  application.  The  mass 
must  be  evaluated  specifically  for  each  application.  The  computational 
technique  developed  in  this  study  meets  this  need. 
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PARAMETRIC  MODEL  OF  SUPERCONDUCTING  ENERGY  STORE 
FOR  SPACE  PROPULSION 


D.l  OVERVIEW 

Superconducting  inductive  energy  storage  is  an  alternative  method  for 
long-term  energy  storage  with  the  inherent  advantage  that  the  energy  is 
immediately  available  in  electrical  form  from  a  low  impedance  source.  Com¬ 
pared  to  storage  inductors  made  frcm  normally  conducting  metals,  supercon¬ 
ducting  inductors  have  no  joule  heating  and  an  essentially  infinite  time 
constant.  However,  there  are  losses  that  vary  with  the  rate  of  field  change 
in  the  superconductor  and  there  are  refrigerator  power  demands  that  depend 
on  these  losses,  on  the  surface  area  of  the  dewar,  and  on  the  currents  that 
must  flow  from  ambient  to  cryogenic  temperatures.  Maximum  current  densities 
and  optimal  field  strengths  are  determined  by  the  critical  parameters  of  the 
superconductor.  For  large  high-field  systems,  the  mass  and  size  of  the 
force  containment  structure  become  the  limiting  factors.  Temperature  may  be 
dictated  by  operational  considerations  or  by  refrigeration  costs. 

The  objective  of  this  note  is  to  develop  a  simple  parametric  model  of  a 
superconducting  energy  store  for  space  propulsion  systems  so  as  to  permit  a 
comparison  with  other  storage  methods  in  terms  of  system  mass  and  energy 
efficiency.  Presently  envisaged  systems  require  pulse  energies  from  frac¬ 
tions  of  a  kilojoule  to  half  a  megajoule,  currents  between  15  and  1200  kA, 
voltages  generally  below  1  kV,  pulse  lengths  of  1  )is  to  6  ms,  and  pulse 
repetition  frequencies  between  0.4  and  19  Hz. 

This  comparison  is  somewhat  complicated  by  the  fact  that  two  different 
types  of  energy  efficiencies  can  be  defined  for  a  cryogenic  system: 

(1)  An  overall  efficiency,  q2  ,  based  on  total  power  consumption, 
including  that  of  the  refrigerator,  and 

(2)  A  conversion  efficiency,  n2 ,  which  only  looks  at  the  'high 
quality*  pulse  power  inputs  and  outputs  and  disregards  the  power 
consumption  of  the  refrigerator  or  the  energy  content  of  any 
'used'  ciryogenic  fluid. 

The  efficiency  q2  will  be  more  or  less  constant  for  a  given  operation, 
but  will  decrease  as  the  repetition  rate  decreases  and  when  long  periods 
of  stand-by  status  are  encountered. 
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D.1.1  Components  of  Snperconductina  Energy  Store 


The  main  components  of  a  sapor  conducting  energy  storage  system  are: 

1.  Hie  superconducting  inductor, 

2.  Hie  current  leada  into  tie  cryogenic  regioa, 

3.  Hie  devar  systoa,  and 

4.  Hie  refrigerator. 

D»2  1HB  SDTBtaXdnCTDB  INDUCTOR 
D.2.1  Opera tint  Coaaideratioaa 

To  fully  specify  a  superconducting  storage  systea  that  will  deliver  a 
specified  energy  pulse  at  a  specified  repetition  rate,  requires  a  number  of 
iterations  on  tie  following  variaoles  or  paraaeters: 

°  Optiaua  field  streagti  for  ainian  size  and  suss, 

0  Field  siape  and  fringing  as  dictated  by  operational  considerations, 

°  Ratio  of  stared  to  delivered  energy, 

°  Type  of  superconductor  and  its  operating  teaperature, 

0  Fora  and  aaouat  of  stabilizer, 

0  Internal  design  of  tie  conductor  to  minimize  hysteresis  and  eddy 
current  lossea,  and 

0  Coolant  flow,  plus  location  and  sizing  of  tie  cooling  channels. 

Clearly,  this  is  somewhat  beyond  tie  scope  of  the  present  analysis. 
However,  it  is  always  true  that  tie  highest  possible  field  will  give  the 
most  compact  storage.  Nevertheless,  whenever  high  currents  and  high  repeti¬ 
tion  rates  are  required  in  a  self-sufficient  env  iroaaent,  refrigeration 
oconcmics,  as  detexainad  by  tie  Carnot  factor,  will  push  the  systea  to  tie 
highest  possible  teaperature. 

Most  superconducting  magnets  built  to  date  operate  at  a  teaperature  of 
about  4.2  E  in  liquid  ieliua  boiling  at  ataospheric  pressure.  Superconduc¬ 
ting  systems  with  dedicated  closed  cycle  refrigerators  permit  operation  at 
any  teaperature  between  2  and  20  E.  For  tie  present  task,  an  operating 
teaperature  of  13  E  is  suggested,  in  conjunction  with  a  refrigerator  output 
at  12  E.  This  is  not  a  conventional  choice  for  either  sagnets  or  storage 
inductors,  but  it  is  feasible  as  Stevenson  (1973)1  has  reported  a  5  T  gas 
cooled  solenoid  operating  at  13  E. 
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D.2.2  Mass  of  Inductor 

In  the  following  discussion  the  usual  notation  for  solenoidal  coils 
will  be  used.  'Die  ratio  of  outer  to  inner  radius  is 

«  ■  *2  /  ax  .  (2-1) 

and  the  ratio  of  half-length  to  inner  radius 

P  *  b  /  .  (2-2) 

Very  large  superconducting  energy  storage  systems  have  been  proposed 
over  the  past  15  years  by  Roger  Boas  and  his  associates  at  the  University  of 
Wisconsin.  The  latest  of  these  reported  by  Eyssa,  Boom,  et  a  1 .2  is  a 
100  hlh  (360  MJ)  study  for  NASA.  It  uses  an  earlier  parametric  analysis  by 
Moses3  and  suggests  10  T,  2  K  coils  at  current  densities  of  250  to 
500  A/mm?.  It  also  accepts  the  difficult  task  of  nonruniform  distribution 
of  superconductor,  as  dictated  by  local  field  strengths,  and  proposes  novel 
pure  tension  structures  for  force  support.  Since  the  design  is  not  con¬ 
cerned  about  stray  magnetic  fields,  it  favors  very  short  solenoids 
(0  *  0J).  where  most  of  the  stored  energy  is  outside  the  winding. 

For  smaller  systems,  one  can  afford  tj  be  somewhat  less  optimistic  and 
consider  longer  coils  with  uniform  current  distributions.  The  mass  of  the 
inductor  can  be  divided  into  three  components:  the  superconductor,  the 

stabilizer,  and  the  structure  providing  force  containment,  support,  and  flow 
channels  for  the  cryogen. 

Realistic  assumptions  would  be  that  the  superconductor  proper,  with  a 
density  of  8  gm/ cm3  will  permit  a  current  density,  JQ,  of  10^  A/ m2  (or  1000 
A/mm2),  »nd  that  the  stabilizer  and  channel  structure  of  density,  dg,  will 
be  present  at  an  'equivalent  current  density1,  J,  of  1  to  5'108  A/m2  (100  to 
500  A/ mm2).  The  value  of  dg  would  be  9000  kg/m3  for  copper  and  stainless 
steel,  and  2700  kg/sP  for  alminum  and  its  alloys. 

A  still  simpler  approach  can  be  based  on  an  analysis  of  Label  1,4  who 
shows  that  for  solenoids  with  uniform  current  density  the  specific  mass, 
M/E$,  is  proportional  to  the  average  density,  d4Vg,  (usually  close  to  dg) 
divided  by  the  square  of  the  maximum  field  on  the  conductor,  times  some 
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(graphed)  coefficients  that  depend  on  a  and  jL  For  small  magnets  (E s<  oj 
MJ),  the  mass  of  the  winding  and  coil  mandrel  predominate s,  and  at  dg  *  9000 
kg/m3 

My  /  (B^B2  )  -  45,000  [  kg'JtJ^T-2  ]  ,  or  (2-3) 

My  /  (Ej-B2*^)  -  5  [  n3-MJ"1*f2  ]  .  (2-4) 

For  large  magnets  (&s>  3  MJ) ,  the  mass  of  the  force  containsent  structure. 
Mg  exceeds  that  of  the  conductor  winding.  My,  and 

(My  +  Mg)  /  (Ej'B2)  -  8,500  [  kg*MJ~1*T“2  ]  ,  or  (2-5) 

(My  +  Mg)  /  (E j'B^'dg  )  -  1  [  m3*MJ~1*T~2  ]  .  (2-6) 

For  a  mostly  aloninnm  structure,  dg  »  3000,  and  a  field  of  5  T  these  equa¬ 
tions  predict  specific  masses  of  600  and  120  kg/MJ  respectively. 

It  is  worth  noting  that  the  minimal  mass  of  the  structure  according  to 
Levy's3  interpretation  of  the  virial  the  ora  comes  to  11.5  kg/MJ  for  Axmco 
21-6-9  stressed  to  200,000  psi.  Label  1  also  graphs  the  specific  masses  for 
a  n amber  of  coils  between  0.05  and  20  MJ  and  shows  values  between  600  and 
200  kg/MJ  with  a  clear  downtrend  in  the  transition  region.  Fig.  1. 

For  space  rated  0J  to  1  MJ  5  T  alminum  alloy  solenoidal  coils  at  13 
E,  I  would  thus  suggest  a  conservative  design  value  of 

(My  +  Mg)  /  E,  *  300  kg/MJ  .  (2-7) 

D.2.3  Surface  Area  of  Inductor 

The  tol  loving  is  not  an  indnctor  design  by  any  means,  but  merely  some 
rough  assumptions  to  arrive  at  a  plausible  geometry  and  a  conservative  upper 
limit  tor  the  surface  area.  They  are: 

1.  The  inductor  is  a  straight  solenoid  of  rectangular  cross-se ct . on, 
with  an  inner  radius  a^< 

2.  The  length-to-dismeter  ratio  £  is  3.  This  approximates  the 
energy  density  inside  an  infinite  solenoid  within  15%.  The 
Length  is  thus  6**^,  and  the  enclosed  field  volme  is 


D— 4 


V 


(2-8) 


=  6  -tt  . 

All  tlie  field  energy  is  contained  within  that  volume  so  that  the 
stored  energy  is 

Es  *  VB2  /  2*hq  ,  and  (2-9) 

*13  =  Mo'Es  1  (3-jpB2)  .  (2-10) 

The  average  current  density,  J,  in  the  winding  volume  is  100 
A/mm2,  or  10®  A/m2,  so  that  the  effective  cross-sectional  area  of 
the  conductor  is 


Ae=I  /  J  . 

(2-11) 

The  nmber  of  ampere  turns  per  meter,  n*l,  are 
long  solenoid  formula  as 

obtained  from  the 

n-I  =>  B  /  n0  . 

(2-12) 

The  winding  thickness  is 

dj  =  n'I/J  =  B  /  (p0*J)  ,  so  that 

(2-13) 

*2  ~  ai+dl  =  al+  B  /  (|i0*J)  and 

(2-14) 

a  =  1  +  B  /  ((ig-J'aj)  . 

(2-15) 

There  is  a  50-mm  space  around  the  coil,  so  that 

Dewar  Radius  =  +  50  nun.  (2-16) 

End  supports  and  plumbing  take  up  300  mm,  so  that 


Dewar  Length 


|5'a1+  300  mm. 


(2-17) 


D.2.4  LOSSES  IN  THE  INDUCTOR 


Changing  fields  cas  se  energy  dissipation  in  s  super  condnct  or.  The 
magnitude  of  these  losses  is  strongly  dependent  on  the  detailed  internal 
1  ay  oat  of  the  coadactor  wire  and  cable.  Fig.  2  shows  the  'sized  matrix' 
structure  for  a  composite  superconductor  specifically  designed  for  milli¬ 
second  pulse  applications.  The  super  conducting  filaments  are  made  as  fine 
as  possible  to  reduce  hysteresis  losses  in  the  superconductor  proper.  Copper 
stabilizer  is  in  close  contact  with  each  superconducting  filament,  but 
neighboring  stabilizer  elements  are  decoupled  from  each  other  by  'shields' 
of  high  resistivity  Cu-Ni  alloy.  This  reduces  eddy  currents  in  the  copper 
stabilizer  and,  in  conjunction  with  a  wire  twist,  limits  similar  induced 
'coupling*  currents  that  are  partly  in  the  filaments  and  partly  in  the 
ma  triz. 

The  approximate  relations  from  the  previous  section  can  also  be  used  to 
estimate  the  total  conductor  vol ume  as 

Vx  >  *  (2‘*1  +  di)  *  (2-18) 

and  that  of  the  contained  superconductor  as 

VQ  3  V1*J  /  JQ  .  (2-19) 

Pinning  centers  are  believed  to  enable  superconductors  to  carry  bulk 
supercurrents  up  to  a  critical  current  density  J c.  For  a  wide  range  of 
fields  the  pinning  force  is  constant,  so  that: 

Pinning  Force  =*  Jc  z  B  .  (2-20) 

The  hysteresis  losses  per  unit  volume  of  superconductor  represent  the 
work  dona  in  pushing  flux  across  pinning  centers  (Bean^),  md  are  for  each 
half-cycle  (charge  or  discharge): 

Q0/Vq  -  •iB-Jc,d0  /  4  [  J  /  ffl3  J  .  (2-21) 

where  d0  is  the  thickness  of  superconductor  perpendicular  to  the  changing 
field  or  the  diameter  of  the  superconducting  filaments. 


These  losses  are  essentially  independent  of  the  rate  of  flax  motion, 
and  are  only  a  function  of  the  magnitude  of  the  field  change,  AB,  the 
critical  current  density,  Jc  =  1.2-J0  (by  design),  and  the  filament 

di ame ter. 

Most  commercial  superconducting  materials  are  readily  available  with 
filament  diameters  down  to  25  pm,  and  5  pm  is  usually  considered  a  practical 
limit.  However,  recent  work  by  Larbal esti er^  promises  1  pm  filaments  with 
improved  current  densities  as  well. 

The  rate  dependent  coupling  losses  in  mixed  matrix  conductor  that  has 
been  given  the  maximum  practical  twist  pitch  of  eight  wire  diameters 
(Laquer8)  are. 


Qr/  ( VQ)  =  (3n/8)**2,AB-Bdot/(p3-D)  ,  (2-22) 

where : 

a  =  Wire  diameter,  typically  0.5  mm 
Bdot  =  dB/dt  ,  Rato  of  field  change  [  T/s  ] 

P3  =  p  of  high  resistivity  web,  typically  10-^Q-m 
D  =  geometry  dependent  coupling  parameter  (0.148  max.). 

In  applying  these  equations  we  shall  assume  that  the  entire  volume  of 
superconductor  is  subject  to  the  full  AB  and  then  divide  by  two  in  Eq  (2-21) 
and  by  four  in  Eq  (2-22).  A  more  detailed  analysis,  again  requiring  a  more 
complete  coil  design,  would  include  integration  of  the  losses  over  the 
actual  field  distribution. 

D.2.5  Limitations  on  Fast  Pulses 

The  above  loss  equations  are  probably  valid  for  field  change  rates  of 
up  to  several  thousand  Tesla  per  second.  For  faster  rates  and  shorter 
pulses,  however,  one  has  to  consider  the  distributed  capacitance  of  the  coil 
winding  and  the  fact  that  the  flux  pinning  mechanism,  or  its  description 
breaks  down  for  very  high  Bdot  values.  If  the  repetition  rates  are  not  too 
high  it  might  be  possible  to  extract  longer  pulses  and  then  shape  them  in  a 
non-superconducting  device. 
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D.3  anUUNT  LEADS 


D.3.1  General  Design  Principles 

The  Inaction  of  cryogenic  laada  is  to  conduct  a  current  f ran  a  power 
source  at  aabient  temperature  to  an  electrical  dewice  at  cryogenic  tempera¬ 
ture.  In  doing  so  they  also  conduct  beat  between  these  temperatur es.  and 
generate  I2*C  losses.  Since  for  a  natal,  thermal  and  electrical  conductivi¬ 
ties  are  fundamental  ly  related,  as  shown  by  the  near  constancy  of  the  Lor o nr 
ratio,  L“p  k/T,  there  is  no  way  to  have  one  withont  the  other. 

The  accepted  procedure  is  to  fabricate  leads  as  'counterflow1  heat 
exchangers  and  size  them  so  as  to  minimize  the  heat  flow  into  the  cryogenic 
enwironment  at  the  rated  current.  If  the  cryoelectrical  dewice  operates  in 
boiling  liquid,  all  energy  influx  is  tahen  up  in  the  latent  heat  of  vapori¬ 
zation  of  the  cxyogen.  All  of  the  boil-off  can  be  channeled  through  the 
leads,  and,  provided  the  leads  are  thermally  anchored  in  the  cryogen,  lead 
boil-off  and  flow  can  be  self-regulating.  No  heat  will  be  conducted 
directly  to  the  device. 

Minimum  energy  input  from  well-designed  counterflow  leads  into  boiling 
helium  amounts  to  1.2  watt/kA/ lead,  as  shown  empirically  and  frem  numerous 
analyses,  e.g.,  Ceil  in  and  Cl  imenko^  and  Efferson.1® 

D.3 .2  Lead  Refrigeration 

In  an  open  system  with  boiling  liquid,  no  use  can  be  made  of  the 
sensible  heat  of  the  cold  escaping  gas.  Most  of  the  lead  refrigeration  is 
'free'.  However,  in  a  closed  system  the  sensible  heat  taken  up  by  the  lead 
flow  gas  in  returning  to  abient  temperature  also  has  to  be  produced  by  the 
refrigerator.  The  more  gas  is  'bypassed*  into  the  current  leads,  from  the 
noxmal  return  flow  through  the  heat  exchanger,  the  mare  difficult  it  becomes 
to  maintain  refrigerator  efficiency. 

Gerhold22  has  discussed  the  problem  in  general  terms.  Using  his  anal¬ 
ysis  together  with  seme  unpublished  Los  Alamos  work,  I  suggest  a  provisional 
formula  tor  the  compressor  power,  P^,  needed  to  refrigerate  a  pair  of  leads 
at  a  steady  current 

PL1  *  20,1  /  T0  I  W  1  (3-1) 
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where  TQ  is  the  minium  gee  temper  store  in  the  system.  It  should  be  noted 
that  this  formula  assumes,  in  effect,  22%  of  Carnot  efficiency.  It  is 
probably  optimistic  and  needs  to  be  examined  in  more  detail.  Neverthel  ess, 
it  serves  to  point  oat  the  tremendous  refrigeration  penalty  imposed  by  high- 
current  leads.  A  pair  of  steady  state  SO  kA  leads  at  12  I  would  require 
83  kW  of  compressor  power. 

For  quickly  varying  and  pulse  currents,  it  should  be  possible  to  time- 
average  the  current  by  relying  on  the  heat  capacity  of  the  lead.  However, 
depending  on  pulse  length  and  repetition  rate,  such  a  scheme  could  introduce 
flow  and  pressure  instabilities  and  degrade  the  heat  exchange  efficiency  of 
the  lead. 

D.3 .3  Lead  Resistance 

Current  leads  optimized  for  minimis  refrigeration  have  a  resistance 
that  produces  a  fixed  voltage  drop  of  about  45  mV  at  the  design  current. 
This  optimum  is  quite  independent  of  the  material,  its  resistivity,  the 
temperature  coefficient  of  resistivity,  or  the  length  of  the  lead.  Nor  does 
it  matter  whether  the  lead  is  of  constant  cross-section  or  tapered  toward 
the  cold  end.  We  thus  have  for  the  optimum  resistance  of  a  pair  of  leads: 

%,  , opt  “  0.09/1  [  MKS  ]  ,  (3-2) 

and  for  the  corresponding  power  dissipation: 

Pu  “  0.09*1  [  MIS  ]  .  (3-3) 

A  SO  kA  current  thus  will  dissipate  4.5  kW  in  a  pair  of  leads. 

D.3. 4  Lead  Mass 

The  length  of  each  lead  is  usually  set  in  the  range  of  0.5  to  1  m  by 
space,  mechanical,  and  insulation  considerations.  Since  R*p*X  /A,  we  get 
from  eq'n  (3-2)  for  each  lead 


(X/A) 

rtnf  3  0.045  /  (I*P.*.  )  [  MKS  ]  . 


Clearly,  &  material  with  tho  lowest  average  resistivity,  Psyg,  over  the 
relevant  temperature  range,  i.e.,  a  high  parity  metal  such  as  A1  or  Ca, 
would  then  permit  the  lightest  leads.  Unfortunately,  such  leads  are  less 
stable  and  less  tolerant  of  excess  currents  (acre  likely  to  burn  out)  than 
those  aade  treat  less  pure  aaterials  with  higher  resistivities  and  saaller 
tea  per  a  tore  coefficients.  The  increased  aass  and  heat  capacity  of  alloy 
leads  also  iaproves  their  theraal  inertia  for  pulse  service,  albeit  at  the 
expense  of  increased  power  dissipation  when  the  current  is  increased  above 
the  optimum. 

An  extreae  choice  would  be  brass  with  a  resistivity  of  8*i0~^:*a  at 
300  C  and  half  that  value  at  10  L  If  the  lead  is  tapered  for  a  constant 
p/A  ratio  along  its  length,  1,  we  get  an  average  cross-sectional  area 

Aavg  *  10‘p1Tg-X*I  -  6-10_7-I  [  MKS  ]  .  (3-5) 

For  0.5  b,  50  kA  brass  leads  the  average  crosssection  would  be  300  ern^  and 
the  aass  of  two  leads  would  be  a  massive  261  kg. 

A  better  material  would  be  Deoxidized  High  Phosphorus  Copper  (DHP)  with 

—  o 

a  resistance  ratio  of  7  and  an  average  resistivity  of  1‘10  °:*m.  A  pair  of 
DHP  50  kA  leads  would  weigh  27  kg  if  straight  with  a  30  cm2  cross-section, 
and  32  kg  with  a  linear  5:1  taper  froa  12  to  60  ca^. 

Information  for  aluainua  alloy  leads  has  not  yet  been  collected,  but 
similar  resistivities  are  available.  I  therefore  suggest  an  interim  formula 
of 

Ml  -  O-2'I  [  kg  /  kA  ]  .  (3-6) 

D.4  DEfAHS 

Dewar  vessels  are  passive  structures  whose  only  purpose  is  to  reduce 
the  heat  flow  from  ambient  temperature  to  a  cryogenic  device,  experiment,  or 
fluid.  In  the  vacuum  of  space  two  nodes  of  heat  transfer  that  are  important 
in  an  atmospheric  envirooaent  are  absent.  These  are  gas  convection  and  gas 
conduction.  Only  radiation  heat  transfer  between  surfaces  and  conduction 
through  solids  remain. 

The  rate  at  which  a  surface  of  area.  A,  emits  thermal  radiation  is 
given  by  the  Stef an^Bol  tzaann  relation 


P  “  WfA'T4 


(4-1) 


where:  s  is  the  eaissivity  of  the  surfaoe,  sad 

o  ■  5.67*10“®  [  f ]  the  rsdistion  coastsat. 

The  aet  exchsage  of  rsdisat  energy  between  two  snrfsces  of  equal  (low) 
emissiv ity  sad  equal  sres  st  Tx  sad  TQ  tskes  place  st 

pj-  •*wA«(T14  -  T04)  •  (4-2) 

Clesrly.  radiative  best  trsasfer  is  ooatrolled  slaost  eatirely  by  the 
high  temperature  surface.  A  300  C  si mian  surface  with  e-0.02  will  radiate 
9.2  W/m^.  The  retura  rsdistioa  free  a  surface  at  10  or  20  K  is  lower  by 
four  or  five  orders  of  Magnitude,  sad  therefore  is  insignificant  in  reducing 
the  heat  leak. 

Conduction  through  solids,  on  the  other  head,  is  expressed  by  the 
Fourier  equation,  which  is  linear  in  the  temperature  gradient: 

P2-  k*A*(T^  -  T0)  /  X  .  (4-3) 

Optimally  designed  dewars  usually  have  equal  heat  losses  through  radia¬ 
tion  and  conduction.  There  is  no  sense  in  trying  to  make  one  heat  leak 
orders  of  magnitude  lower  than  the  other.  As  a  matter  of  fact,  any  signifi¬ 
cant  additional  energy  losses  in  the  systea  should  also  be  considered  when 
designing  the  dewar  insulation. 

D.4.1  Multilayer  Insulation 

The  ideal  insul  ation  woul  d  consist  of  a  nnnber,  N,  of  concentric,  mass¬ 
less,  low-emissivity,  isothermal,  floating  radiation  shields  in  a  perfect 
vacuum.  Such  an  arrangement  would  reduce  the  heat  leak  given  by  equation 
(4-2)  by  a  factor  of  N+l. 

Commercial  multilayer  insulation,  known  as  'super-insul ation',  tries  to 
approximate  the  ideal  by  interleaving  thin  aliainum  foils  with  fiberglass 
paper  spacers  or  by  using  multiple  wraps  of  crinkled  altminized  mylar.  The 
latter  is  the  lighter  option  and  also  can  be  altminized  with  strips  or  small 
squares  so  as  to  minimize  eddy  currents  induced  by  changing  magnetic  fields. 
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Ho  temperature  profile  in  those  insulations,  ho* over,  is  not  as  one  sight 
expect  tor  a  series  of  pare  radiation  shields,  but  becoaes  linear  below  70 
L  His  indicates  that  conduction  through  contacts  is  predominant  ia  the 
low  reaper  a  turo  region,  effectively  shorting  out  the  insulation  that  would 
be  pro* i dad  by  the  radiation  shields.  On  the  other  hand,  vith  the  absence  of 
gravity  in  a  space  dewar,  it  is  conceivable  that  contact  forces  and  the 
resulting  heat  floe  could  be  reduced  through  proper  design. 

Altaiaized  aylar  insulation  is  available  in  bulk  fora.  Application 
data,  published  by  National  Research  Corporation12,  suggest  that  the 
aaterial  is  aost  effective  when  applied  at  16  layers/cm  (40  layers/ inch)  and 
then  will  have  a  heat  flux  between  300  E  and  20  E  of 

Pj  -  18.9  /  N  [  (Ta”2*  layers-1  ]  ,  (4-4) 

provided  the  naaber  of  l^ers,  N.  is  between  30  and  130. 

The  mass  of  the  insulation  is 

Mj  -  8J4*10~3*N  [kg*a-2*  layers  ]  .  (4-5) 

As  long  as  the  dewar  will  only  be  operated  in  space,  its  inner  al  uainum 
al  1  oy  shel  1  need  be  no  aore  than  5  aa  thick  and  its  outer  si  el  1  a  mere  1  mm. 
This  would  give  a  aass  of 

Md  -  16  JL  [  kg’ a- 2  ]  .  (4-6) 

For  short  pulses,  eddy  currents  in  the  dewar  shell  may  be  undesirable  and  a 
fiberglass  epoxy  dewar  of  siail  ar  aass  may  have  to  be  used. 

D.5  KEWSEXATOM 

Information  on  spacerated  refrigerators  is  rather  scarce  in  the  open 
literature,  probably  for  classification  and  proprietary  reasons.  However, 
Strobridge1^  has  done  a  thorough  survey  of  cryogenic  refrigerators.  After 
taking  out  the  inevitable  Carnot  factor 

”C  “  T0  ’  (Ta  “  V  '  (5-1} 
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the  cabined  thermal  end  mechanical  efficiency,  does  not  depend  on  the 

temperature,  TQ,  between  1.8  end  90  L.  Fig.  3  taken  frca  his  paper  shows 
that  it  ia  only  a  function  of  the  refrigeration  capacity  (f rigiwatts),  PQ. 
For  saall  refrigerators  with  large  surf  ace- to-  vol  me  ratios  is  saal  1.  It 
reaches  15  %  for  kilowatt  aachines  and  levels  off  at  about  33  %  for  aegswatt 
uni ts. 

Strobridge  also  gives  a  relation  for  costs,  Cgf  (which  have  to  be 
adjusted  for  inflation  since  1969)  as  a  function  of  the  ccapressor  power  at 
aabient  teapersture,  P 


*  PQ  /  (nc’ng  )  tnd 
CR  *  48*Pa0*7  [1969  i,  watt  ]  . 


Finally,  we  can  extract  the  following  relations  frca  his  plots  of  aass  and 
volute  against  frigiwatts  for  refrigerators  in  the  10  to  30  K  range: 


*  60‘Pg*7  [  kg.  watt  ]  ,  and 


VR  *  0.055*Po0,85  [  a3,  watt  ]  . 


D.5 .1  Refrigeration  Alternatives 

Stored  cryogens  present  an  alternative  to  aechanical  refrigeration,  but 
only  for  relatively  short  aissions  and  low  refrigeration  powers.  The  latent 
heat  of  vaporization  can  be  used  effectively  above  the  triple  point  tempera¬ 
ture,  Tp.  Various  cryogens  have  been  used  in  this  manner  for  cooling  either 
equipment  or  radiation  shields.  For  liquid  hydrogen  the  cooling  power 
aaounts  to  about  450  kJ/kg.  The  triple  point  is  at  13.8  E,  mandating  a 
slight  increase  in  the  operating  temperature  of  the  inductor  or  a  short  span 
(magnetic  ?)  refrigerator.  Detailed  analysis  of  these  options  is  beyond  the 
scope  of  this  note. 


D.6  EXAMPLE 

We  suaaarize  the  use  of  the  parametric  model  with  the  example  of  a  500 
kj,  50  kA,  5  T  superconducting  storage  inductor  at  13  E  and  refrigeration  at 
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12  K,  recognizing  that  most  of  the  numbers  are  intended  to  be  conserve tiv e. 


i.e.  ,  on  the  high  side. 


Coil  length: 


i  ge  t 

the  mass 

M*  +  Mg  »  150 

Inner 

radius : 

tl  -  0.139  m 

l  : 

2*b  »  6^  -  i 

Vind. 

thictaus : 

di  ■  0.040  m 

Outer 

radius : 

*2  -  0.139  +  1 

Dewar 

radius : 

ij-*-  o.o5  ■  o.: 

Dewar 

length : 

2-b  +  0.3  -  1 

ice : 

A„  -  1.96  m3 

Cond. 

Volume: 

VA  -  3  .33  *10-: 

Dewar  surface: 


This  does  not  include  any  structure  and  corresponds  to  a  mass  of  about 
100  kg  compared  to  the  above  conservative  value  of  120  kg. 


From  2-19  Supercond  Volume 


VQ  -  3.33-10"3  m3 


From  2-21  Hysteresis  Loss  for  AB  ■  5  T 

do  *  5  |im 

Q0  *  21  J  /  half  cycle 

if  the  entire  conductor  sees  the  full  field  change,  otherwise  divide  by  2 
and  get  10.3  J. 

At  a  rep  rate  of  1  Hz  this  dissipates 

PQ  »  2’10.5  =*  21  frigiwatts. 

(Same  loss  for  charge  and  discharge). 

From  2-22  Coupling  Loss  for  the  given  typical  values  and 

AB  -  5  T 

Bdot  »  2000  T/s  (1  ms  pulse) 

Qg  *  1660  J  /  discharge  cycle 

if  the  entire  conductor  sees  the  full  field  change,  otherwise  divide  by  4 
and  ge  t  413  J. 

At  a  rep  rate  of  1  Hz  this  dissipates 

■  413  frigiwatts. 

(Assuming  that  Bdot  for  recharge  is  negligible). 
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From  3-1  Lead  Refrigeration  for  50  kA 

Pli  “  83  kW  Compressor  Power 
Divide  by  2  for  average  current  of  25  kA. 

PL1  "  42  kW 

From  3-3  Lead  Dissipation  for  50  kA 

?L2  “4.5  kW  of  Pulse  Power 
Divide  by  4  for  average  power  dissipation. 

From  3-6  Lead  Mass  for  50  kA 

Ml  -  10  kg 

Stas  of  average  low  temperature  power  dissipation  from  hysteresis  and 
coupling  (neglecting  'ordinary1  eddy  currents  and  deferring  the  lead  prob¬ 
lem) 

Pq+Pl*  436  frigiwatts 

If  these  losses  are  present  all  the  time,  dewar  losses  would  need  be  no 
less  than  44  w. 

From  4-4  Layers  of  Superinsol ation  to  achieve  44  V 

N  *  18.9'Ad/44  =  0.84  , 

or  a  single  radiation  shield. 

A  simple  40  layer  insulation  system  would  have  a  heat  leak  of 

Pj=  18 .9'Ajj/  40  =0.93  frigiwatts. 

This  would  minimize  standby  power  requirements. 

From  4-5  and  4-6  Dewar  Mass 

Mj+Mjj  -  0.7+31.8  =  32.5  kg. 

From  5-1  Carnot  Factor  at  12  I 

»)c  =  0.042  or  4 .2% 

From  Fig.  3  Refrigerator  Efficiency  (iterated) 

=  0.13  or  13% 

From  5-2  Compressor  Power  (neglecting  leads) 

Pa*  437/(0.042-0.13)=  80  kW 

Total  Compressor  Power 

PL1  +  pa  “  122  KW 

From  5-2  and  5-4  combined  Refrigerator  Mass 
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Total  Mass 


Mg”  5680  kg 


1 50+1 0+3 2+56 80  =■  5872  kg 

To  calculate  efficiencies  sad  t\2  mOT 8  inforaatioa  is  needed  oa  operatioa 
aad  duty  cycles  aad  oa  tlxe  priae  power  supply. 

Clearly,  the  biggest  power  aad  mass  requirements  are  associated  with 
the  size  of  the  refrigerator,  which,  ia  turn,  is  determined  by  the  curreats 
ia  the  leads  aad  the  coupliag  (eddy  current)  losses  ia  the  composite  super¬ 
conductor. 
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Fig.  1.  Specific  Mass  of  Superconducting  Coils,  from  Lubell  (1973) 

Fig.  2.  Mixed  Matrix  Superconductor  Designed  for  Millisecond  Pulse  Applica¬ 
tions,  from  Laquer  (1976) 

Fig.  3.  Efficiency  of  Low  Temperature  Refrigerators  as  a  Function  of 
Refrigeration  Capacity  (Frigiwatts),  from  Strobridge  (1969) 
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